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The Design of Active Distributed Filters 


By ROBERT E. PARKIN 


(Manuscript received July 13, 1967) 


We develop a synthesis method which leads to simple realizations of low- 
pass and high-pass active filters using distributed resistance-capacitance 
networks. In the low-pass filter there is one distributed network terminated, 
through a negative impedance converter, in a load consisting of a lumped 
network or another distributed network. Thus we can obtain a Chebyshev- 
like low-pass filter which has a faster cutoff than a four pole Chebyshev 
realization with the same ripple. In the high-pass filter design an addi- 
tional feedback element is involved. We can obtain a complex pole pair 
to be placed anywhere in the complex frequency plane except the real fre- 
quency axis. Design charts are given. 


I. INTRODUCTION 


In a recent book, P. M. Chirlian covers the main filter synthesis 
techniques using distributed resistance-capacitance networks.’ This 
paper takes a more practical design-oriented approach. Consistent 
with the original terminology of Wyndrum,” we call the uniform dis- 
tributed RC network URC. The aim here is to design active distributed 
filters with minimal circuit element complexity. We consider the two 
basic types of URC shown in Figure 1, the single URC, and the common 
ground high-pass double URC. The low-pass filters are synthesized 
with a single URC and a Linvill realization® and the high-pass filters 
by a common ground high-pass double URC and a Yanagisawa 
realization.* 
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(a) (b) 
Fig. 1— The circuit symbols for (a) the single D.C.R. and (b) the cross-con- 
nected double D.C.R. 


Il. THE LOW-PASS FILTER DERIVATION 


2.1 The Poles of the Transfer Function 


The transfer function for the circuit shown in Figure 2(a) can be 
shown to be° 


yo () 


MS cosh 6 — Yvr se + 2ys(@sinh 6 — yyz cosh @) 





where 
VJ procol’, 


Yv_ = TolLyt , 


6 


ens — 2s/ToL, 


19 is the resistance per unit length, cy is the capacitance per unit length, 
Lis the length of the network, and :p is the complex frequency variable. 
If zg = 0, the poles are given by 


A a 
Uwe = tonh 6 2) 


ToL. 


a 


(a) . (b) (c) 


Fig. 2— The basic low-pass filter circuit. 
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Notice that since p = 6’, the negative real axis in the p plane is 
the imaginary axis in the @ plane. 


2.2 Lumped Load 


If the load y, is the RC combination shown in Figure 2(b), and 
if 6 = x + jy, equation (2) becomes 


1 ‘ } 
Uwe = Gt (e+ jy W = eee (3) 
where 
_ ki =(Cr 
v= <i and W = a 


An infinity of roots of equation (8) lie on the j axis of the @ plane 
for any U or W, and these poles are free to move’ a maximum distance 
of r. When W = 0 (the case of resistive loading) the dominant pole 
pair in the 6-plane can lie at +j7 when U = 0, can move to the origin 
when U = 1, or move along the real axis towards +o as U > 0. 

If U = oa, 


1 


WO en ale 6 


(4) 
and a pole pair caused by the capacitance loading, the loading pole 
pair, lies on the real axis in the 6-plane for W > 0. As W — 0, the loading 
poles tend to + and in the limit W = 0 do not cause instability. 

It is possible to choose values for U and W so that the unbounded 
loading and dominant poles can combine and produce roots of equa- 
tion (3) for 2 and y nonzero. The graph of U and W for values of x 
and y nonzero is given in Figure 3. 

The stability plot of U and W is shown in Figure 4; the plane of 
reference is the p plane. 

Feeding the URC by a current source, the poles of V,/Z,; can be 
shown to be” 





Yv~ = 6 tanh @. (5) 


Again, there is a unique solution of equation (5) for x and y nonzero. 
For the load considered previously, the solution of equation (5) is 
given in Figure 5. 
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Fig. 3.— The values of U and W for complex poles at 9 = x + jy when zs = 0. 


2.3 Distributed Load 
If load y, in Figure 2(a) is the shorted URC shown in Figure 2(c), 


for 25 = 0 the transfer function is given by 
a Sa i pe a a 
Bia sab & + days] sinh 6— VELA ooh 6 | 
(6) 


where Arp and Bcp are the resistance and capacitance per unit length 
of the load. 


cosh 6 — 
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When zs = 0, the poles are given by 


tanh 6 = +~/A/B-tanh (~/ AB: 8). (7) 
It is possible by suitable choice of A and B to place a pole anywhere 
in the 6 plane. In addition, depending on the values taken by A and 
B, more than one complex pair of roots can occur for x and y nonzero, 
but no matter what values A and B assume, positioning one complex 
pole pair fixes all the poles of the system. When a root is positioned 
near the x = y line in the @ plane with all other roots in the stable 
region, to obtain an equiripple filter, the other complex roots (if they 
occur) are close to the x = 0 line. 
The solution of equation (7) for the section of the 6-plane near the 
x = 6 line is shown in Figure 6. 


III. LOW-PASS FILTER DESIGN CHARTS 


3.1 Lumped Load 


From equation (1) and the load of Figure 2b, it can be seen that 
the position of the complex pole is determined by U, W, and zyg. By 
choosing these three factors, the ripple in the passband is fixed. When 


U 
NEGATIVE REAL UNSTABLE COMPLEX OR 
‘ AXIS POLE POSITIVE REAL AXIS POLES 
STABLE 
es 
POSITIVE REAL 
AXIS POLE 
5 STABLE COMPLEX 
NS OR NEGATIVE 
UNSTABLE REAL AXIS POLES 
N 
70.5 BN, 0.5 1 
—-W Me EON 
NEGATIVE REAL POSITIVE REAL 
AXIS POLE AXIS POLE 
STABLE UNSTABLE 
S-U 


Fig. 4— Stability plot of the low-pass lumped load filter when zs = 0. 
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the ripple is small the complex pole pair is near the origin. When the 
ripple is large the complex pole pair is a long way from the origin and 
close to the x = y line, the jw axis in the p plane. The design chart of 
Figure 7 was obtained by setting 6 = x + j(x+8x) for some value of 
Zyg, Where 8x is small, and determining U and W for this condition. 
With W held constant, U was decreased until equiripple conditions 
were achieved. This was repeated for a complete range of values x and 
zyg. In addition, Figure 7 shows the conditions for constant p,, where 
p- is the 3 dB cutoff frequency. 


Example 


Suppose it were required to synthesize a low-pass filter with a fast 
cutoff having a 3 dB cutoff frequency of 75 kHz, source resistance of 
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Tig. 5 — The values of U and W for complex poles at ® = x + jy when zs = ©. 
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Fig. 6 — The values of A and B, the distributed load factors, for complex poles 
at ® = x + jy when 2s = 0. 


5K Q and working into another circuit with a driving point resistance 
of 10K Q, and the maximum ripple tolerated in the passband is to be 
1 dB. 


If p, is the normalized cutoff frequency, 





ToLeoL = aa 10° ; 
and 
_ 5X 10° 
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Fig. 7— The lumped load low-pass filter design chart. 


One factor can now be chosen from Figure 7, zys, ToL, or col. Let- 
ting roL = 5 X 10+, then 


Zvng = 0.1 
po. = 30 Hz 
U = 0.2335 
W = 0.150 


Thus coL = 8260 pF. The value U is the total normalized resistor 
load comprising R; and the 10K © driving point impedance of the 
next stage. The realization of the circuit is shown in Figure 8. The 
theoretical and measured frequency responses of this circuit are 
shown in Figure 9; notice the close correlation between the two. The 
negative impedance converter used in making the measurements is 
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YL = 50kn Col = 8260pPF 


1239 PF 


Fig. 8 — A low-pass filter circuit having lumped loading. 


the NIC of Drew and Gorski-Popiel.? The rate of fall of the response 
at cutoff is 20.75 dB per octave, and the rate at the 10 dB point is 
26.4 dB per octave. 


3.2 Distributed Load 


The design procedure proceeds as before, by placing the complex 
pole on the x = y line and increasing A until equiripple conditions 
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Fig. 9 — The response of the filter circuit of Fig. 8. 
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Fig. 10 — The distributed load low-pass filter design chart. 
are achieved. The poles of the transfer function are given by 


VBfA__ _ (cosh 0+ O-sinh O-2ys)_ 
tanh (VAB-6) inh 6+ O-cosh O-2ys) TB ©) 


If (AB)%-6 = x + ja’, the left side of equation (8) can be split 
into real and imaginary parts, so that 


a _ _ sinh 2’ cosh 2’ 
6  sinz’ coszx’ (9) 

Thus (AB)*-6 is calculated, and substituting in equation (8), A 
and B are known explicitly. 

Figure 10 shows the design chart for a low-pass filter with dis- 
tributed load. It was found that p, was independent of A for any value 
of B, so that cutoff frequency lines were not shown in Figure 10 but 
were given in Figure 11. 
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Fig. 11— The relation between p. and B. 


Example 

Using the low-pass filter with distributed load to synthesize a filter 
having the same requirements as in the previous example, but this 
time neglecting the effect of loading of the next stage, then 


A = 0.323, 
B = 0.446, 
Zvg = 0.38, 
rol = 1.514 X 10%, 


Col = 0.0274 uf, 


for p, = 31 Hz. The circuit is given for this filter in Figure 12. The 
measured and theoretical frequency responses are given in Figure 13. 
Again, notice the close correlation between the theoretical and: meas- 
ured frequency response. The slope of the response at cutoff is 21 dB 
per octave, and at the 10 dB point it is 26.75 dB per octave. 





YoL = 15.15 kn Col = 27400pF 


Fig. 12 — A low-pass filter circuit having distributed loading. 
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Fig. 13 — The response of the filter circuit of Fig. 12. 


3.3 Elimination of Inband Loss 


The inband loss obtained in the examples of Sections 3.1 and 
3.2 can be eliminated or set to any value required by adjustment 
of the active device. For instance, if —yz is realized by yr, and a 
current inversion negative impedance converter described by 


vin = Kut 
try = Kteat ? 
the equivalent form of equation (1) for z, = 0 is given by 


V2 1 


(cosh 9 — deen ) 


k is the factor controlling the inband loss. 
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IV. THE HIGH-PASS FILTER 


4.1 Transfer Function Singularities 


The parameters of the common ground high-pass URC have been 
developed by Hager.*® Imbedding this element with the current negative 
impedance converter (current NIC) and y4 — yz networks as shown 
in Figure 14, the transfer function can be found as 











Vv 
on tanh Vv 1~2yva 
a v ¥_\__ 
2Yyvat2yvn—1 tanh zteve| uve +2 ya tanh 3) ~ tanh ¥ 
(10) 
where 


w= V2proL’. 


Notice that replacement of the current NIC by a voltage NIC does 
not yield the same result.® It can be seen from Equation (10) that 
the loss at high frequencies is zero if and only if zg ~ Oasp— o. 

When 2g = 0, the poles of equation (10) are given by 


1 Vv 
Yves = 2 —— 9 tanh — Una (11) 


and the zeros by 





1 Vv 
Yva ~ 9 (ts 1) : (12) 


where 


YNA = ToLYs and YUNA = ToLYz . 





Fig. 14— The basic high-pass filter circuit. 
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4.2 Lumped Load and Feedback 


If ys and yz have the form shown in Figure 14, for a second order 
zero at the origin, equation (10) yields 


isa) 


yp 
(Ee) 


where Wp = ya/(pcol). 
Letting ¥ — 0 in equation (13), Wr = 14 Substituting in equation 
(11) 


—-1-WW, =0 (13) 


1 Vv a 
We 5 Sean 7 3. (14) 


There is a unique solution of equation (14) for x and y nonzero, when 
etjy=y. 


4.3 Distributed Load and Feedback 


The common ground high-pass URC is still used as the central filtering 
element in this high-pass filter, but the lumped elements are replaced 
by URGs as shown in Figure 15. Substituting for yy, and yyz in equa- 
tion (10), and for zys = 0, 











Dis lA al | eae = Se 
tan 3 1 4 cosh — 1, 
Vo _ tanh WY © cosh ¥ + 1 
Ve 2(v sinh ¥+2 cosh ¥+2) ‘pen. 4y sinh p (15) 
sinh W( cosh W +1) tanhd —a(vsinh»y-+2(cosh y +1) 
A( cosh err rs Ga Vil 

where 

v= Qpr Col 

Y= V/2pArBol? = ~/AB-3 


V/2par Berl? = +/aB-d. 


—~ 
I 





Fig. 15 — The completely distributed high-pass filter circuit. 
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For second order zeroes to occur at the origin, B = 14. 
Notice the similarity to the previous result. 


V. HIGH-PASS FILTER DESIGN 


A high-pass filter can be realized from the circuit of Figure 14 by 
replacing yz with a resistor Rz and capacitor Cy, in series, and y4 by 
a capacitor. Unfortunately, to place the complex poles where re- 
quired could require Rz; > 0, C, < 0. To make Cz > 0, connect an- 
other resistor R, on the other port of the NIC. 

To assess the best position for the complex poles, three factors can 
be used which, when minimized simultaneously, produce desirable 
properties in the filter. 


(t) Make the slope at cutoff as near as possible to 12 dB per octave, 
the low-frequency rate of change of attenuation. 

(it) Make the difference between the frequency at which the response 
crosses the constant attenuation line passing through the high-fre- 
quency attenuation, and the frequency at which a 12 dB per octave 
line drawn at the cut-off point would cross the high-frequency line, 
a minimum. 

(iti) Make the overshoot a minimum. 


These three factors are shown in Figure 16, where G, = 12, the 
slope of the response line at cutoff in dB per octave, G. = the fre- 
quency at which the response attenuation is equal to the high-fre- 
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Fig. 16 — Three factors governing the quality of the high-pass filter. 
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Fig. 17 — A lumped load and feedback high-pass filter. 


quency attenuation —1.25 x (cutoff frequency), and Gs; = peak 
overshoot. Combining these three factors together as 


G=GQ4+@4+6, 


where 


os 40G, 
> 1.25 X (cutoff frequency) 





and minimizing G', the best filter under these circumstances if zs = 0 
and the overshoot is limited to 1.25 dB, is shown in Figure 17. The 
response of this filter is shown in Figure 18. Notice that G, is slightly 
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Fig. 18 — The response of the circuit of Fig. 16. 
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modified to G{ to normalize for the cutoff frequency, and to introduce 
a loading factor. 


VI. CONCLUDING REMARKS 


It is evident that for best filtering characteristics, the URC must 
be connected to use its transmission line properties. The low pass 
filters use a URC as a transmission line, and the negative-valued 
loads serve merely to boost the high-frequency responses and extend 
the cutoff frequencies in a way not unlike that employed in video 
amplifiers, where an inductor is often used to compensate the high- 
frequency capacitive drop in the response.? 

The completely distributed low-pass filter has a faster cutoff than 
the partially distributed filter. The difference in response is caused 
by the slightly different noncomplex pole positions near the origin; 
the further away the first negative real axis pole in the p plane is 
from the origin, and so closer to the other noncomplex poles, the 
further away the complex poles must be from the origin and nearer 
the limit of stability line to obtain an equiripple filter. Thus, the 
cutoff frequency is higher, the effect of the nondominant poles lying 
on the negative real axis of the p plane on the transfer function at 
frequencies about cutoff is more pronounced, and a faster cutoff re- 
sults. 

For frequencies up to the 20 dB point, to obtain an equivalent re- 
sponse to that of the distributed low-pass filter in lumped circuitry, 
would require at least five resistors and five capacitors. 

To improve the response of the high-pass filters would require a 
higher order zero at the origin than the second order zero now ob- 
tained. Until this is achieved with a minimum of active devices, the 
distributed high-pass filters of the type considered here will show 
less promise than the low-pass filters. 

Sensitivity has not been considered in this paper; to do jacties: to 
the subject would require a depth of study sufficient to warrent an- 
other paper. However, for the low-pass filters, the information con- 
tained in Figures 7 and 10 could provide a basis for future sensitivity 
studies. 
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Design of Monopulse Antenna Difference 
Patterns with Low Sidelobes* 


By E. T. BAYLISS 


(Manuscript received August 9, 1967) 


The flexibility of modern monopulse radar antenna systems makes 
possible the independent optimization of sum and difference patterns. 
The two parameter difference pattern, developed here for the circular 
aperature antenna, is designed to have nearly equal sidelobes similar to 
those of the Taylor sum pattern. The difference pattern is asymptotic to 
a model difference pattern which has the greatest slope (angle sensitivity) 
for a given sidelobe level. The model function is unrealizable because it 
has uniform sidelobes which are infinite in extent. The two parameter 
difference pattern ts realizable and is expressed in a Fourier—Bessel series 
of N terms in a manner similar to Taylor’s treatment of the sum pattern. 
The other parameter, A, controls sidelobe level. 

Comprehensive tables of the Fourier—Bessel coefficients are given for 
both the circular aperture series and the difference pattern series. Directivity 
and angle sensitivity are investigated and found to have maximum values 
that decrease as sidelobe level decreases. The monopulse system performance 
using the asymptotic difference pattern and the Taylor sum pattern compares 
favorably with a maximum likelihood angle estimation system. Development 
of a line source difference pattern is presented in the appendix. 


I. INTRODUCTION 


Monopulse radar systems have, in recent years, achieved a high 
degree of flexibility in their antenna patterns. This flexibility is the 
result of development of multihorn feed structures and array anten- 
nas. It now appears feasible to optimize independently the sum and 
difference patterns of the monopulse system. A good deal of work has 
been done to improve the sidelobe performance for sum beam pat- 

* This work was supported by the Army Materiel Command under Contract 


DA-30-069-AMC-333(Y) through the Nike-X Project Office, Redstone Arsenal, 
Alabama. 
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terns. However, less attention has been paid to improving the angle 
sensitivity and sidelobe performance of the circular aperture dif- 
ference pattern. 

The objective of this investigation was to develop a difference pat- 
tern which possesses characteristics that are compatible with those of 
the Taylor sum pattern. Specifically, the goal was to obtain the maxi- 
mum angle sensitivity commensurate with a given sidelobe level. Low 
sidelobes are desirable in both sum and difference patterns for the 
suppression of near-target clutter, ground clutter, and jammers. Re- 
quiring large angle sensitivity and low sidelobes for the difference 
pattern is analogous to. requiring a narrow beamwidth and low 
sidelobes in the case of the sum beam. The pattern that meets these 
requirements must be produced by a reasonably well behaved aper- 
ture illumination. 

This paper solves the problem of generating difference patterns by 
using a technique that parallels Taylor’s approach to the sum pattern 
design. The problem is first detailed in terms of making good angle 
estimates. It is attacked by applying a general synthesis technique 
to the approximation of a model difference pattern. Finally, the re- 
sulting asymptotic difference pattern is examined for angle sensi- 
tivity. Pertinent design information is presented. 


II]. THE PROBLEM—-ANGULAR ESTIMATION IN NOISE USING 
ANTENNA DIFFERENCE PATTERNS 


A monopulse system can make a maximum likelihood angle esti- 
mate? under one of the following assumptions: 


(t) The primary source of noise is spatially and temporally uncor- 
related noise from the radiation field 

(it) The primary source of noise is thermal noise which is independ- 
ent in the sum and difference channels. 


The angle estimate is made by correlating the sum and difference 
channel outputs. The sum and difference patterns required to make 
the maximum likelihood estimate unfortunately have quite high side- 
lobes (for example, the difference pattern first sidelobe is only 14.5 
dB below the sum beam maximum). For applications where clutter 
or active noise sources are a problem, such high sidelobes are clearly 
unacceptable. Techniques used to suppress sidelobes can be expected 
to reduce the angle sensitivity. (This tradeoff is examined in Section 
IV.) Thus, the problem is to suppress antenna sidelobes in both the 
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difference and sum patterns without severely degrading the angle 
sensitivity of the system. 

The objective of this paper is to provide a unified approach to dif- 
ference pattern design for the monopulse system. Specifically, the 
objectives are to: 


(t) Design a difference pattern with equal sidelobes, similar to side- 
lobes of Taylor sum pattern. 

(%1) Solve the practical problem of generating the aperture function 
for the circular aperture. 

(12) Balance sidelobe level and angle sensitivity to get the best 
angle sensitivity for a given sidelobe level. 

(tv) Compute and tabulate such design parameters as angle sensi- 
tivity, directivity, sidelobe level and, of course, the aperture generat- 
ing functions. 


Ill. THE SOLUTION—-AN ASYMPTOTIC PATTERN FUNCTION WITH 
UNIFORM SIDELOBE LEVEL 


3.1 Aperture and Pattern Functions 


The aperture screen concept is very useful in the analysis of planar 
antennas. It is well known that from the tangential E and H fields on 
a surface enclosing a source one can, in principle, calculate the radi- 
ated fields. If the enclosing surface is an infinite plane, then some 
useful simplifications can be made. The equivalence principle and 
image theory® can be used to determine the field in the source-free 
half-space (2 > 0) from its tangential components on the plane (2 = 
0 in Fig. 1). For this special case the E-field can be determined solely 
from the E-field on the plane.* 

Plonsey* has pointed out that it is possible to describe the field in 
the Fraunhofer region for all 6 < «/2 where the usual Kirchhoff scalar 
diffraction approximation’ is good only for small 6. 

The aperture and pattern expression can be put in standard form 
by making the following change of variables: 


2a. v3 
w= > sin 8 and P=7e (1) 


where wu is the angle variable measured in standard beamwidths, p 
is the normalized aperture variable and a is the radius of the antenna, 
(Fig. 1). It is of interest to consider the case where the vector field 
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Fig. 1— Aperture geometry. 


in the aperture is factorable: 


Ep z | = i +a,)00,¢); pS . 
: 0; elsewhere 


where y is a complex constant, a; is a unit vector and g(p, ¢) is the 
complex scalar aperture function. Notice for linear polarization, y = 0, 
and for circular polarization, y = +7. Under the above assumptions 
and with the usual Fraunhofer approximations (that is, p< R,1<« kR) 


E(R, 6, 6) = =e {ar X [a, X (a. + va)}} (22) Fe, $); 0< : 
(3) 


where 


Fu, ®) =~ [ae | of, &) exp Lup cos (® — ep dp. (4) 


The pattern function, F(u, &), which is the Fourier-Hankel transform 
of the aperture function, g(p, 9), will be the starting point for the 
next section, 
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3.2 Synthesis Method for Circular Apertures 


The procedure for synthesizing an arbitrary pattern function with 
a circular aperture is an application of Woodward’s technique to the 
circular aperture synthesis problem.°® 

The aperture function can be represented as a Fourier-Bessel series 


0 


(Pp, ¢) = 2 4 ye)? D4 (uerp)} ee (5) 


k=-0 1 


This series has orthogonality and completeness properties for appro- 
priately chosen y,: . It converges in the mean to an arbitrary aperture 
function, j(p, ¢), when the complex coefficients, a,; , are found by taking 
the inner product of the function with each orthogonal function. 

The Sturm—Liouville boundary conditions give rise to the eigenvalues 
for the angular and radial eigenfunctions: 


(2) The angular eigenfunctions must satisfy periodic boundary condi- 
tions. Thus their eigenvalues are the integers k. 
(7) The radial eigenfunctions must satisfy one of two cases: 


1. J i,(uii™) = 03 wir >0,k 20,1 = 0, 1, 2; Pn 
2. Cd (Mer) + Mart { (ua) = 0;c 2 0, c ae k > 0, wer > O for 
L = 0,1, 2, --- except when c = k = 0 then poy = 0. ; 


The selection of case 1 or 2 and the constant c determine the set of 
positive distinct eigenvalues to be used in the radial eigenfunctions 
of (5). 

The aperture function expansion that allows nonzero boundary 
values (case 2) is to be preferred in the expansion of equal sidelobe 
patterns because: 


(i) As Taylor’ has shown, the slowest far sidelobe decay rate results 
from nonzero boundary values. 

(it) The “central region” (| u | < uzx) zeros are confined to a smaller 
region. This means that better sidelobe control can be obtained in 
the central region for a given main-beam performance and a fixed 
number of terms, JV. 


Thus, choosing case 2, let c = 0. The synthesis of a difference pattern 
allows simplification of (5) because of the symmetry required. The 
difference pattern requires an antiphase aperture [one sign reversal, 
F(u, ®) = —F(u, ® + 7)], thus only first order Fourier components 
are required (k = +1). Thus choosing k = 1, the boundary conditions 
reduce to* 


* Primes denote differentiation with respect to the argument. 
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J'(u7r) = 0 b= '0,.1,.4>+ (6) 
where the k subscript on p,; has been dropped. See Table I for values 
of Mi. 

To further simplify the choice of coefficients, let the difference axis 
be g = @ = 0. This requires setting a,, = —a_,, = B,/2 and gives - 


N-1 
i(p, ¢) = cose 2, BiFs(uip); pr. (7) 


The N-term approximation to the required function @ or F is designated 
by a tilde g or F. The pattern function expansion corresponding to 
g(p, ¢) in equation (7), that is, transformed by (4), is 


ud) eo 
— 4" 





F(u, 6) = j cos ® Ds B,J (ur) 7 (8) 


The terms of (8) 


ud{ pire, 


Ji(usm) 7 eae 





are shown in Figure 2. 

The pattern function series (8) has the important “sampling” 
property which many approximation techniques use. If (8) is set equal 
to the required pattern function, F(u, 6), at the sample point (zu, , 0) 
then the mth coefficient, B,, , can be evaluated because all other terms 
of (8) are zero there. Equation (8) is evaluated at (u,, , 0) by applying 
L’Hoéspital’s rule, substituting Bessel’s equation and the boundary 
value (6) to obtain 


Flim 50) = Bub 5 Fiwue)(L — Grin) *). @) 


TasLe [—Bersseu FunctTIon Zuros, J,/(ui7) = 


l M1 l Mi l br L M1 
0 0.5860670 5 5.7345205 10 | 10.7417435 15 15.7443679 
1 1.6970509 6 6 .7368281 11 11 .7424475 16 16 .7447044 
2 2.7171939 7 7.7385356 12 | 12.7480408 17 17 .7450030 © 
3 3.7261370 8 8.7398505 13 | 138.7435477 18 18 .7452697 
4 47312271 9 9.7408945 14 | 14.7439856 19 19 .7455093 
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Tig. 2— Bessel series pattern function terms. 


Now setting (9) equal to P(u, &) at the sample points, the coefficients 
of (8) are evaluated 


A Zo ee eee 
Di = Mn 5] 0) je? Jy "(u,.0)(1 _ (TH) “] " 


m=0,1,2,---,N-—1. (10) 


In Section 3.4 (10) is used to evaluate the coefficients from sample 
points on the model function, F,,, and thus determine the aperture 
function (7) and pattern function (8). 


3.3 Construction of a Model Pattern Function 


The previous section provides a method of expanding any realizable 
pattern function. In this section a model difference pattern function is 
developed. It will serve as a model for the construction of the asymp- 
totic difference pattern in the next section. Thus its role is analogous 
to Taylor’s “ideal function.”* 

The main characteristic desired in a model function is maximum 
angular sensitivity for a given sidelobe level. It has been shown’ for 
an array that the pattern function that produces the greatest angular 
sensitivity with a given sidelobe level must also have sidelobes of 
equal height. This theorem is analogous to that proved by Dolph for 
the sum pattern.® 

Two important differences occur when the above result for an 
array is extended to the case of the continuous aperture: 
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(t) The pattern function has maximum slope for a given sidelobe 
level (not maximum sensitivity). 
(71) The pattern function is not a realizable antenna pattern. 


These two characteristics result from requiring the difference pat- 
tern to have sidelobes of equal height. These equal sidelobes now are 
infinite in extent (for all vw). Any pattern function whose sidelobes do 
not decay at least as w°/? is not realizable because it requires an un- 
bounded aperture function. Because the function is unrealizable its 
angular sensitivity is undefined. Nevertheless it does have the maxi- 
mum slope for a given sidelobe level and will serve as a model func- 
tion for synthesis of a realizable asymptotic difference pattern in Sec- 
tion 3.4. 
Thus, the main characteristics of the model function are: 


(1) Maximum slope at the origin for a given peak-to-sidelobe ratio. 
(iz) Sidelobes of uniform magnitude and infinite in extent. 


Price and Hyneman pointed out that an exact functional form for 
such a function is not known. What follows is a method for construct- 
ing a very close approximation to the equal sidelobe difference pat- 
tern function. 

The ideal sum pattern used by Taylor’ is of the form 


F.(u) = cosrV/u? — A? (11) 


where A is a parameter that determines the sidelobe level. This func- 
tion has equal sidelobes that are infinite in extent. Taking the deriva- 
tive of the ideal sum pattern a difference pattern is obtained. 


Fy(u) = rw — A®)'?yu sin r(w? — A?)*™, (12) 


The first few sidelobes of this function are not of equal height as 
can be seen from the dashed curves in Figure 3. However, it is pos- 
sible to modify the first few sidelobes of the above function so that 
they are equal to the asymptotic sidelobe level. The zeros of the 
above function are given by 


= ‘ ee ; (13) 
+(A?>+n)!; n=1,2,3,--- 


A model difference function can be constructed by moving the 
first T zeros on either side of the origin so that the sidelobes are made 
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Fig. 3— Model difference pattern compared with unmodified function for 
sidelobe levels of 25, 35, and 45 dB. 


equal. The model difference function is given by 


2 


Py(u) = Fw) 3 (14) 


rere ee 
where the zeros, z,, are moved to &,. The location of the new zero &, 
which make the sidelobes equal can be found by an iterative proce- 
dure.® Very good results were obtained by moving only the first four 
zeros on either side of the origin. The constructed model functions 
are shown in Figure 3 as solid curves. 

To expedite future computations the sidelobe parameter A and the 
zeros were fitted as polynomial functions of the sidelobe level in dB. 
Fourth degree polynomials were fitted to values of A, é,, and the loca- 
tion of the difference peak, po, that were obtained from the iterative 
computation. Curves of the parameters along with their polynomial 
coefficients can be found in Figure 4. 

The use of fitted polynomials inevitably leads to errors in the re- 
sults. The effect of these errors is more pronounced in the low side- 
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Fig. 4— Model difference pattern parameters. 


Polynomial Coefficients 





Polynomial 
Name 
Co C1 Ce Cs C4 
A 0.380387530 | —0.05042922! —0.00027989)| —0.00000343] —0 .00000002 
Er 0.98583020 |—0.03338850} 0.00014064) 0.00000190} 0.00000001 
£2 2.00337487 |—0.01141548! 0.00041590| 0.00000373) 0.00000001 
&3 3.00636321 |—0.00683394| 0.00029281} 0.00000161) 0.00000000 
4 4.00518423 |—0.00501795;} 0.00021735; 0.00000088; 0.00000000 
Po 0.47972120 |—0.01456692| —0.00018739)] —0.00000218} —0.00000001 


lobe pattern. For example, notice the depressed first sidelobe in the 
—45 dB pattern of Figure 3 (and later in Figure 8). If good fidelity is 
desired for low sidelobe patterns it is advisable to use the calculated 
zeros directly and avoid use of the polynomials. 


3.4 Synthesis of Asymptotic Difference Patterns 


Now the problem of approximating the model function using the 
synthesis method described in Section 3.2 is considered. In this sec- 
tion an asymptotic form of the model function is developed. The 
asymptotic form is expressed in terms of the Fourier-Bessel series. 

In Section 3.3 the model function for the ideal difference pattern 
Fy, 1s generated in equations (12) and (14). The outstanding feature 
of this model function is its equal sidelobes. Fy; is an entire function 
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of uw and is completely described by its zeros, 


0; n= 0 
ZL, = \+6,;3 n=1,---,T7. (15) 
hte +n); n=T+1,-+ 


For large | wu | the zeros approach +n asymptotically. The model 
function is not realizable because it does not have decaying sidelobes 
(nonfinite energy). - 

The steps in constructing the asymptotic difference pattern are 
illustrated in Figure 5. The behavior of the model function in the far 
sidelobe region is changed so that it conforms to the asymptotic be- 
havior of a realizable pattern function with decaying sidelobes. To 
do this, the zeros, Z,, of Fy, for alln = N are moved to u,. These 
new far sidelobe zeros are zeros of the truncated pattern function 
series (8). This series must represent a realizable pattern function. 

These zeros also satisfy (6). By changing only the far sidelobe zeros 
we have retained the essential model function characteristics in the 
central region. This central region behavior is embodied in the central 
zeros, Z, for n < N. The Nth zero pair of the model function, +Zy , 
and the Nth zero pair of the realizable function, --uy , do not coincide. 


|F(uw)| 











CENTRAL «<———-—— | ——————> FAR SIDELOBE 
REGION REGION 


U<pny My<u 


-3/ 
u>/*asyMPTOTE 
REALIZABLE 
FUNCTION 


/ 
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Mn+e 





Ky+3 


TP, o2 oZ aZ oZ 


Fig. 5— Construction of asymptotic difference pattern. 
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Usually, the distance between Zy_, and py is somewhat greater than 
one. Such a situation will give rise to a rather large sidelobe at the 
transition between central region and far sidelobe region. 

To eliminate this high transition sidelobe the zeros of the central 
region will be dilated by a factor « which makes Zy and py coincide 
(16) 


c= EX. 

Zy 

For sufficiently large N, o will always be greater than one. 
The steps illustrated by Figure 5 are: 


(1) A realizable pattern function with decaying sidelobes replaces the 
model function in the far sidelobe region. 

(zt) The model function (dashed curve) is dilated so that the zero 
(¢Zy) coincides with the first zero of the far sidelobe region at u = py. 


The asymptotic form of the model function is essentially complete. 
It can be expressed as a cononical product of its zeros 


F,(u, ®) = Cur cos ® u E — (2) | is E — (4) | (17) 


where C' is a constant. The second product term must be reduced to 
closed functional form for this expression to have any practical value. 
The zeros in the second product of (17) satisfy J{(u,r) = 0. Thus 
the second product term in (17) in closed form is equal to the first 
term of (8) with all zeros removed for |u| < py. 


ff -()}-enwefflr-@)). os 


Combining (17) and (18) 
N-1 au 2 
- (Gy 


i [s- (Cy 


F(u, ®) is a realizable asymptotic difference pattern. It is nearly equal 
to the model function in the central region (that is, | wu | < uy). In fact, 
as the size of the central region increases without bound, the asymptotic 
function approaches the model function 


Lim F,(u, ®) = Fy,(u, ®). (20) 
N20 








Fu, &) = C (cos &)2ruJ {(ur) (19) 


ANTENNA PATTERNS 635 


The limit can be deduced from (17) by noting that asymptotic behavior 
of uy (6) gives 
3 

— ___#w at 4 m N. 21 

(+N Nn Peres 5 1<T< (21) 
Now as N is increased « — 1. Thus as N goes to infinity (and « — 1) 
the second product term (17) disappears and only the first remains. 
The first product term (with « = 1) is exactly equal to the model 
function. Thus F’,(u, &) is truly asymptotic to Fy,(u, ®). 

The asymptotic function (19) is expressed as a Fourier—Bessel series 
using equation (10) to evaluate the coefficients. First evaluate F.(u, ®) 


at the sample points (n,, , 0) 
N-1 | u 2 
- (4) | 
I LJ Jy. Ji{(ur) 


F(tn ) 0) eae F(un ’ 0) = C2tpm N-1 2 U . 
i= (s) [eG 
1=0 Ki ais 


l4~m 


o 








m=0,1,---,N—1. (22) 


Evaluating the limit and substituting Bessel’s equation yields 








Fy(um 0) = C(rpn) TL — rpm) “VI a(unr) 5 f - (12) 


m=0,1,---,N — 1. (23) 


Notice also that because of the construction of the asymptotic function 
F (un, 0) = 0 for m = N, Finally the coefficients from (10) are: 


lee on 
|-cj2, E ~ (=) | 
tian) 2 ’ NS OL oe YN el 


1=0 by 
l4~m 


0 : m=N,N +1,-::: 








where y,, is given by (6), Z, by (15) and o by (16). The constant C 
is evaluated so that the peak of the asymptotic difference pattern is 
unity, that is, F.(opo , 0) = 1. 

The pattern function expansion (8) is identically equal to the asymp- 
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totic difference pattern, F,(u, ®), because only the first N coefficients 
(24) are nonzero, thus the series terminates. Furthermore, the expan- 
sion (7) is equal to the aperture function of the asymptotic difference 
pattern. Thus equations (24), (8), and (7) define the solution to the 
problem of finding a difference pattern compatible with the Taylor 
sum pattern. (A similar function is developed for the line source in 
the appendix.) This realizable two-parameter difference pattern is 
asymptotic to the model function. The parameter A controls the 
sidelobe level of the pattern. The parameter N controls the size of the 
central region; hence, it controls the degree of approximation to the 
model function. 


IV. THE RESULT—DESIGN DATA AND PATTERN CHARACTERISTICS 


4.1 Pattern and Aperture Functions 


The usefulness of a pattern function is determined largely by the 
availability of design data about the pattern. Coefficients B,; were 
calculated in equation (24) for sidelobe levels from —17.5 dB to —45 dB 
at 2.5 dB intervals and for N = 3 to 30. Coefficients for selected values 
of N appear in Tables II and III. The coefficients are normalized to 
make the difference pattern peak equal to unity. 

The coefficients are used, in equation (7), to generate the aperture 
function with ¢ = 0 (Figures 6 and 7). Some observations can be 
made regarding the effects of the parameters on the aperture function: 


TABLE I] — CoBFFICIENTS TOR bax PATTERNS 


N= 5 7 8 11 14 17 21 27 30 30 30 30 
COEFFICIENTS SLL= -17.5 ~20.0 ~22.5 25.0 —27.5 —30.0 ~32.5 35.0 —37.5 -40.0 42.5 -45.0 
BO 0.79738 0.79422 0.78677 0.78215 + 0.77598 0.76878 0.76183 0.75508 0.74706 0.73851 0.73009 0.72200 
Bl 0.14654 0.20764 0.30017 0.36076 0.42971 0.49981 0.56485 0.62432 0.68963 0.75578 081714 0.87429 
B2 0.13474 0.11727 0.06774 0.04336 0.01691 —0.00460 —0.01848 ~0.02572 ~0,02877 -0.02588 0.01699 —0.00303 
B3 —0.15632 —0.18338 —0.13874 -0.12031 -—0.09011 --0.06174 0.04006 —0.02421 -0.01019 0.00022 0.00651 0.00973 
B4 0.09650 0.18907 0.16174 0.16100 0.13248 0.10080 0.07503 0.05502 0.03628 0.02125 0.01079 0.00375 
BS —0.15105 -0.15245 —0.18152 —0,16116 —0.12935 -0.10140 ~—0.07855 -—0.05604 -—0.03749 —0.02412 ~0.01464 
BG 0.08042 0.11498 0.18236 0.17722 0.14914 0.12134 0.09701 0.07139 0.04984 0.03410 0.02272 
B7 0.05983 —0.16641 -0.18258 —0.16223 -0.13703 --0.11266 —0.08443 —0,06016 -—0.04233 -0,02932 
B8 0.13589 0.17757 0.16881 0.14884 0.12604 0.09575 0.06900 0.04929 0.03485 


ies] 
o 


~0.09434 —0.16273 ~—0.16888 -0,15679 ~0,13730 -0,10558 —0.07661 --0.05523 --0.03952 
0.04678 0.13912 0.16253 0.16080 0.14645 0.11398 0.08307 0.06024 0.04344 
- 0.10843 ~—0,15009 —0.16080 -—0,15342 —0,12094 -0.08841 --0.06437 --0.04667 
0.07294 0.13214 0,15681 0.15811 0.12642 0.09261 0.06762 0.04921 
~0,03554 —0,10956 —0.14893 -0.16046 ~0.13036 —0.09565 —0.06999 -0,05107 

0.08348 0.13742 0.16040 0.13271 0.09750 0.07146 0.05226 

--0.05531 -0.12261 ~0.15794 —0,13343 -—0,09813 —0.07202 -—0.05276 

0.02676 0.10500 0.15309 0.13250 0.09753 0.07167 0.05258 

— 0.08518 —0.14594 -0.12993 —0,09572 -0.07040 - 0.05172 

0.06387 0.13661 0.12575 0.09270 0.06825 0.05019 

~—0.04188 —0,12528 —0.12000 —0.08853 —0.06524 —0.04803 

0.02018 0.11216 0.11275 0.08324 0.06140 0.04525 

~ 0.09749 ~—0,10409 —0.07691 -0.05677 —0.04188 

0.08144 -0,09395 0.06947 0.05133 0.03790 

—0.07256 —0.08964 -—0,06646 ~—0.04925 -0.03649 


MONDDMOARBDIDWVWAWAVD 
NNNNSoSSSSa Boo 
WNHeEDOOMRNDHEWNH=AS 


B24 0.04915 0.07314 0,05417 0.04010 0.02967 
B25 ~0.03185 -0.06044 -0,04481 —0.03320 —0.02460 
B26 0.01530 + 0.04762 0.03534 0.02621 0.01944 
B27 : "0.03472 -0,02579 -0,01915 -0.01422 
B28 0.02199 0.01635 0.01215 0.00903 


B29 -0.00976 -0,00725 -0.00539 -0.00400 
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TABLE [J] — CoEFFICIENTS I’oR DIFFERENCE PATTERNS 


COEFFICIENTS \ SLL= -17.5 
BO 0.79738 
Bl 0.14654 
B82 0,13474 
B3 ~ 0.15632 
Ba 0.09650 

COEFFICIENTS \ stL= 17S 
BO 0.81073 
B1 0.02204 
B2 0.26549 
B3 —0.37699 
B4 0.43812 
BS 0.45439 
B6 0.42252 
B7 —0.34848 
BB 0.24219 
B9 0.11903 


20.0 
0.78751 
0.26831 
0.06408 
~0.09756 
0.06431 


-20.0 
0.80063 
0.15201 
0.16480 
—0.25968 
0.30936 
—0.32444 
0.30349 
—0.25158 
0.17579 
—0.08696 


22.5 
0.77806 
0.37445 
0.01520 
—0.05612 
0.04125 


-22.5 
0.79056 
0.26817 
0.09021 
--0,17324 
0.21485 
—0.22918 
0.21625 
—0.18048 
0.12694 
—0,06326 


25.0 
0.76893 
0.46772 
~0.01692 
0.02743 
0.02488 


-25.0 
0.78064 
0.37274 
0.03653 
—0.11027 
0.14597 
0.15972 
0.15259 
—0,12849 
0.09111 
~0.04580 


Ne=5 


27.5 
0.76033 
0.55036 
~0.03613 
~0,00808 
0,01343 


-30.0 
0,75219 
0.62400 
—0.04535 
0.00445 
0.00556 


N=10 


-275 
0.77095 
0.46739 
0.00030 
—0.06513 
0.09622 
—0.10940 
0.10637 
0.09065 
0.06493 
—0,03297 


-30.0 
0.76154 
0.55343 
—0.02345 
0.03349 
0.06071 
—0.07324 
0.07303 
0.06323 
0.04586 
~0,02356 


—32.5 
0.74458 
0.69001 
—0.04690 
0.01209 
0,00029 


-32.5 
0.75248 
0.63195 
~0.03549 
—0,01200 
0.03576 
—0.04751 
0.04915 
0.04349 
0.03205 
-0.01670 


—35.0 
0.73722 
0.74921 
—0.04252 
0.01623 
~—0.00311 


—35.0 
0.74367 
0.70372 
—0.03846 
0.00191 
0.01858 
0.02942 
0.03219 
—0.02935 
0.02208 
~0,01171 


—37.5 
0.73042 
0.80278 
~—0,03364 
0.01796 
—0.00517 


-37.5 
0.73523 
0.76955 
—0,03404 
0.01023 
0.00708 
0.01690 
0.02027 
~0.01932 
0.01495 
0.00810 


—40.0 
0,72392 
0.85116 
—0.02137 
0.01804 
—0,00631 


~40.0 
0.72715 
0.83004 
—0,02362 
0.01451 
—0.00031 
0.00841 
0.01201 
—0.01227 
0.00987 
~0,00550 


~42.5 
0.714792 
0.89522 
—0,00658 
0.01706 
--0.00680 


42.5 
0.71935 
0.88566 
0.00834 
0.01596 
—0.00478 
—0,00281 
0.00638 
~-0.00737 
0.00629 
—0.00364 
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—45.0 
0.71213 
0.93516 
0.01005 
0.01544 
0.00686 


-45.0 
0.71189 
0.93694 
0,01085 
0.01549 
-0.00721 
0.00074 
0,00262 
~0.00401 
0.00380 
—0,00232 


first, when JN is fixed (Figure 6), one sees that increasing the sidelobe 
level decreases the peak but increases the edge illumination. Con- 
versely, when the sidelobe level is held constant (Figure 7), the effect 
of increasing N reduces the peak only slightly but increases the edge 
illumination tremendously. One of the prominent characteristics of this 
type of pattern synthesis is that the edge illumination increases as V 


increases. 
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Fig. 6 — Circular aperture functions for N = 10. 
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Fig. 7 — Circular aperture functions for —30 dB sidelobe level. 


The pattern functions exhibit a more straightforward dependence 
on the two parameters. The sidelobe envelopes are plotted in Figure 8 
and 9 that correspond to the aperture functions in Figures 6 and 7. 
Only the envelope is plotted with a O to mark the sidelobe peak. 
Obviously the sidelobe parameter controls the maximum sidelobe level. 
The suppression of the first few sidelobes for very low sidelobe levels 
(Figure 8) is caused by using too small an N. When nondecreasing 
sidelobes are encountered the pattern is in a nonasymptotic region 
indicating that there are not enough terms in the series. 

The effect of changing N can best be seen by keeping the sidelobe 
level fixed (Figure 9). As N is increased the sidelobes become more 
nearly equal as well as increasing the region affected. If we examine 
the sidelobe more closely we see that an inflection occurs at about 
u = N. The inflection results from the transition between the central 
region where the zeros have been modifled and the far sidelobe region 
where the natural decay envelope is u~*”” (Figure 10). 

The asymptotic behavior of a given pattern function for large u can 
be examined in a manner similar to Taylor’s’, giving the asymptotic 
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form of the pattern function as 


(n)<" Te + 1) 
rv/2 yo? 


h(, ®) exp iu = ; = 2) aaa, exp | — ie — ; — s) 
) 2 = 
for |w|— oo and Reu> 0. (25) 


F(u, &) ~ 


The above notation follows Taylor’s except that the aperture function 
is defined over the entire aperture not just on a line. The function 





























SIDELOBE LEVEL IN DECIBELS 
































8 
u, STANDARD BEAMWIDTHS 


Fig. 8 —Sidelobe envelope of pattern functions for N = 10. 
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h(p, ®) is the analytic part of the aperture function, that is, g(p, ®) = 
h(p, ~)(r” — p*)*. For nonzero boundary values a = 0 and A(p, ¢) = 
g(p, ¢). Evaluating (25) for the difference pattern gives 


Feu, ©). ~ V2 ofr, &) EA. (26) 


As one might expect when dealing with aperture-limited functions 
the asymptotic form of the pattern function depends only on the 
aperture function behavior at the boundary of the aperture. The 
asymptote in Figure 10 was computed from equation (26). 

To pursue the topic of asymptotic behavior further leads to con- 
sidering the superdirectivity ratio and ‘“Q” of the antenna. A rough 
qualitative definition would describe the superdirectivity ratio as the 
ratio of total power flow through the aperture to the power that is 
actually radiated. The definition of Q is the ratio of energy stored in 
the near field (evanescent waves) to that radiated per cycle. Increasing 
the superdirectivity increases the reactive component of power flow 
(also Q) and thus increases the ohmic losses in the antenna. 




















SIDELOBE LEVEL IN DECIBELS 














8 
u, STANDARD BEAMWIDTHS 


Fig. 9—Sidelobe envelope of pattern functions for —30 db sidelobe level. 
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Fig. 10 — Asymptotic sidelobe behavior for —25dB sidelobe level and N = 11. 


Large superdirectivity ratios are to be avoided (expect in very small 
antennas) for a number of practical reasons: high ohmic losses, high 
error sensitivity, and narrow bandwidth. How then is superdirectivity 
affected by the choice of difference pattern? Superdirectivity increases 
as the integrated power of the pattern function in the region u > 2a/h 
increases. This means that for a given sidelobe level superdirectivity 
increases with NV. However, for a given pattern function, superdirectivity 
decreases with increasing antenna diameter. Because of the undesirable 
effect of the reactive component, the parameter N should be chosen 
so that N < 2a/). 

The asymptotic behavior of the pattern function as N is increased 
is best illustrated by the dilation factor, o (Figure 11). The dilation 
factor decreases to one as N is increased toward infinity. In the non- 
asymptotic region (Figure 11) the dilation actually increases with NV. 
Because of severe sidelobe anomalies it is best to avoid using patterns 
with N in the nonasymptotic region. For N = 30 the pattern dilation 
is less than 3 per cent for all sidelobe levels. 
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Fig. 11 — Pattern dilation factor, o. 


4.2 Pattern Characteristics 


A figure of merit by which the performance of a pattern can be 
measured against other patterns is a great help in selecting a specific 
pattern. Two such figures are directivity and angular sensitivity. 

The directivity function is defined so that the aperture function 
can be used to find the total power 


4 
2.2 = | Fu, ®) ki 
Gu, @) = 2 ——* —_____ 
| dy i | g(p, ¢) |’ pdp 


This definition is based on Silver’s work,!° and makes the following 
assumptions: 





(27) 


(t) The effect of the obliquity factor is neglected. 
(it) The antenna is large in wavelengths. 
(tz) The effects of ohmic losses in the antenna are neglected. 


The maximum directivity of a uniform-phase circular aperture is 


4 22 
G =r (28) 
The maximum directivity of a difference pattern naturally occurs at 
the peak of the difference beam (uw, 61) rather than at boresight. 


The directivity at the pattern maximum normalized by the maxi- 
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mum aperture directivity is defined as the relative directivity, n: 


4 2 
_ Gu, &) _ ~ | F@ , &) | 


(Fo 





re - (29) 
| de | | g(p, ®) | p dp 


where (ww; , &;) is the pattern maximum. 

This quantity is similar to aperture taper efficiency. Computation 
of relative directivity can now be carried out for the asymptotic dif- 
ference pattern by noting that the coefficients computed for Tables II 
and III were normalized to make F(u,, 6;) = 1 and that the aperture 
function series is orthogonal. 


1 = SYS Bt Mal - cueny) (30) 


Relative directivity in dB (10 log 7) was plotted as a function of NV 
and sidelobe level (Figure 12). It can be seen that there is one N that 
gives a maximum relative directivity. The maximum is obtained when, 
as N increases, the power removed from the main beam by narrowing 
is just offset by the power added to the sidelobes by raising more 
sidelobes. The maximum relative directivity that can be achieved in 
any difference pattern (antiphase aperture) is —2.47 dB.” 

The most important characteristic of a difference pattern is its angle 
sensitivity. The angle sensitivity function is defined” using equation (27) 


+ 
kw = ot volt /volt/aid. BW (31) 





nm IN DECIBELS 
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Fig. 12— Difference pattern relative directivity, 7. 
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where ®, is the difference pattern axis. The angle sensitivity of a 
pattern is K(0) which is the maximum value at boresight. The maxi- 
mum angle sensitivity for any aperture is produced by a linear odd 
aperture function’? (that is, g = jv). For a circular aperture its pat- 
tern function is 


Fu, 0) =-7 ne (32) 


This pattern has —11.6 dB sidelobes and angle sensitivity 


Ky = 5 G3 volt/volt/std. BW. (33) 


Again it is convenient to normalize to this maximum angular sensi- 
tivity and thus define relative angle sensitivity to be 


KO) 224 


Ko T 


where F (u, &) is normalized to unity. 
By evaluating the partial derivative of the difference pattern func- 
tion at w = 0 the relative angle sensitivity is calculated as 


oF (u, ®,) 
Ou 








(34) 


u=0 


ree J. 
fg eB: Furr) | (35) 
1=0 Mi 


Relative angle sensitivity in dB (20 log 5) was plotted as a function 
of N (Figure 13). The curves are very similar in behavior to those for 
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Fig. 18 — Difference pattern relative angular sensitivity, 6. 
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relative directivity as might be expected from equation (35). Certainly 
it is not unreasonable for the derivative of the directivity function 
at boresight to be strongly influenced by the magnitude of the peaks 
on. either side. The dependence on sidelobe level is greater than for 
relative directivity. The dashed curve in Figure 13 is the maximum 
angle sensitivity (Sn.x) for each sidelobe level. The pattern coefficients 
are given in Table II for 12 sidelobe levels with N chosen to maximize 6. 

The decrease in angle sensitivity as the target moves off boresight 
is of interest in determining the linearity and useful angular sector 
of operation. The normalized angular sensitivity function, K(u)/Ko is 
plotted for six 6... patterns in Figure 14. Although the boresight 
sensitivity is higher for high sidelobes than for low sidelobes, it also 
‘decreases more rapidly as a function of angle. 


4.3 System Performance 


The sum and difference signals derived from the antenna just 
discussed can easily be put in a form that allows processing by any 
of the three types of angle detection systems.1* The system generally 
used makes the angle estimate by forming a ratio of the difference to 
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Fig. 14— Normalized angular sensitivity function, K(u)/Ko. 
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sum signal. The statistic of this ratio was recently treated’* with 
regard to maximum likelihood estimates of angle. For the high signal- 
to-noise case the variance of the maximum likelihood angle esti- 
mate? 1* is shown to be: 





2 f/omr\-t 
a, 2th (22) (36) 


eS Ben 


where A and B correspond to G3 and K in our notation,* Z is the signal 
energy and N, is the noise power per cycle of receiver bandwidth. 

Manasse had shown earlier’ that the minimum standard deviation 
for the maximum likelihood angle estimate made using a circular 
aperture is 

(ox)or = (oi)aee(=2) = 2 (222) " sta. BW (37) 
mn r \N, 
where 2/,/N, is the output peak signal-to-noise power ratio in a per- 
fectly matched receiver with optimum antenna system (that is, uni- 
form aperture function). 

The effects of antenna illumination on system performance can be 
expressed noting that realized signal energy F is also affected by antenna 
sum directivity. Thus H = HE yn:L where L accounts for losses such as 
filter mismatch loss, integration loss, and transmitter and propagation 
losses. Then (36) expressed in standard beamwidths is 





_ Gu) (2#ansl) 
o, = K(u) N, std. BW. (38) 
Now by noting that 
KO) _ G3(0) __ : Ku) — Gs(u) 
ke = 6 and G, = 77>; K(0) = 1 and Gs(0) = 1 
for |u| <1, and for a circular aperture Ky = (1/2)G?, then (38) becomes 
ie ean 
% = TSN, std. BW. (39) 
Thus at high signal-to-noise ratios the antenna dependence is 
pd: 
Oy we 5 


This dependence shows that the difference pattern angle sensitivity is 
of prime importance in determining the over-all system angle accuracy. 


* The subscript = refers to sum pattern parameters. 
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The optimum system from the standpoint of angle estimation is one 
with 6 and ns equal to unity. As we pointed out before such a system 
would have unacceptably high sidelobes. The price of lower sidelobes 
is reduced sensitivity (Figure 15). The top curve shows the maximum 
relative directivity, mmx(2), obtainable from a Taylor sum pattern.* 
Similarly, the NV can be selected to maximize relative angle sensitivity, 
Snax » for the asymptotic difference pattern (Figure 15). 
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Fig. 15 — Relative directivity and angular sensitivity for maximum sensitivity 
pattern. 


The relative directivity (7) of the 6,,,, pattern is shown at the bottom 
of Figure 15. Notice that it is 3 dB below the sum pattern relative 
directivity. Thus, to equalize the sidelobe clutter returns from the 
highest sidelobes in the sum and difference patterns, the difference 
pattern sidelobe level should be 3 dB greater than the sum patterns 
sidelobe level. For example, if a sum pattern with —30 dB sidelobe 
level is chosen (NV = 8) the difference pattern sidelobe level should 
be —27 dB, and 6,,,, is —1 dB when N = 14. 

One of the major reasons for suppression of sidelobes is the presence 
of clutter. Although clutter can be a problem throughout the pattern 
it is often more prevalent at low angles. If clutter does not occur near 
the main beam it might be advantageous to choose a small value of NV 
so that the sidelobes decay rapidly (Figure 9). If clutter is very bad 
it may be necessary to taper the transmitter illumination, too. 


* Where N corresponds to n in Taylor’s notation, 
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A second reason for suppression of sidelobes is the presence of active 
noise sources. These are the essence of the radio astronomy problem. 
For the radar problem they are jammers (electronic countermeasures). 
They generally are considered to be point sources which may appear 
anywhere in the pattern. Hence the requirement to minimize the 
maximum sidelobe level. 

In selecting a difference pattern a major practical consideration is 
that the aperture function be realizable. Remember that the aperture 
function represents the resulting field on the face of the aperture (that is, 
after all ohmic losses and mutual couplings have been accounted for). 
Certainly one could not choose an aperture function with a lip or 
spike at the edge for use with a horn fed antenna or array. Even in 
an array of independent elements it usually is not desirable to let the 
edge excitation exceed the peak. 

These considerations along with constraints of the system problem 
such as required accuracy, clutter, noise, and so on, influence the 
selection of sidelobe level and NV. 


V. SUMMARY 


Full theoretical and design documentation has been given to a new 
two-parameter difference pattern. The salient features of the pattern 
are asymptotically equal sidelobes and near optimum angle sensitivity. 
The two parameters, A and N, give complete control of the sidelobe 
level and decay behavior. The difference pattern was designed to be 
compatible with monopulse systems using the Taylor sum beam 
illumination. 
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APPENDIX 


Asymptotic Difference Pattern for a Line Source 


The construction of an asymptotic difference pattern for a line 
source can be carried out in a manner similar to that described for 
the circular aperture. The equations for the line source difference 
pattern, which parallel those in the body of this report, are outlined 
here. 

The line-source-to-pattern-function transform is a finite Fourier 
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Foy = i g(a)e*"* dx (40) 
where wu = 2a/) sin 6, 6 is measured from the normal, 2a is the length 
of the line source and x is normalized to z+/a. Expand g(x) in a sine 
series because an odd function is required 

N- 


1 
Brsm m2; —T 
g(x) = = l Mit, 


IA 


x 


IIA 


‘: (41) 
0; 


elsewhere. 
A nonzero boundary value is again required for the truncated series. 


The eigenvalues must be u, = 1 + 1/2. Now transforming the series 
by (40) and setting wu, = 1+ 1/2 yields 


N-1 U 
PN ces (—1)'u cos u, 
Pe) = 25 oe 


(42) 
The desired model function, F3;(w) (see equations 12 and 14) is again 
described by its zeros. 








The dilation factor is 
_ ay ~N+3. 
co 7 (43) 
The asymptotic function is given by the canonical product 
N-1 2 ee) 2 
Fu) = C a - (4y} a - (4) 44 
) i I oZ, sg Bi ( ) 
where C is a constant. 


The asymptotic difference pattern in closed form is: 


N~-1 u 2 
pe (2 | 
| N-1 U 2 f (45) 

=) 

I Mi 
The function may be expressed in series form by evaluating (45) 
at U = p» and setting it equal to (42) to find B,,. Evaluating the 
series (42) at w= p» gives 





F(u) = Cu cos (ru) 





F(un) = Lim F,(u) = jrB,, . 


Urn 


(46) 
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Substituting p, = n + 1/2, the coefficients are: 


lee m + a) | 
I] E = ( oZLn ; 


ee Ce ges Se 
By = Lt \r¥4 


l=0 
lm 


m=0,1,2,---,N-—1 
0; m=N,N+1,-::: 


(47) 


As N is allowed to increase without bound it may be noticed that o— 1 
from above and the zeros of F,(u) approach those of /'s,(u). Thus 


Lim F(u) = F4,(u). 


N02 


For small N the values of B,, can be evaluated by hand by using the 
model function parameters from Figure 4. 
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Error-Controlled High Power Linear 
Amplifiers at VHF 


By H. SEIDEL, H. R. BEURRIER, and A. N. FRIEDMAN 
(Manuscript received July 21, 1967) 


Two amplifiers, each covering a different band in the VHF range, were 
constructed for test set applications. These amplifiers were required to meet 
stringent specifications beyond the general scope of previous art. These 
spectfications related to the simultaneous availability of relatively large 
output power, tight linearity, minimal delay distortion across a 40 per cent 
band, high input and output return loss, large dynamic range, and tume 
stability. 

The performance problems were solved by using: 


(t) Quadrature couplers in a corporate structure array to provide a 
high multiplicity of “paralleled” transistors, yielding both the 
power and return loss capabilities. 

(it) The emitter follower configuration as the basic amplifier element, 
providing broad band, high level performance owing to the large 
self degeneration of that configuration. 

(iit) A “feed-forward,” instead of a feedback, system for error control 
and noise cancellation. 


I. INTRODUCTION 


There are times when problems arise in device implementation that 
appear to be the product of the sheer perversity of the system de- 
signer. Worst of all, he can justify his need—as in the need to de- 
velop a solid-state 10-watt amplifier between 60 and 90 Me, with 
low intermodulation distortion, time delay distortion, and gain rip- 
ple, with a high return loss, and with excellent time stability. Sub- 
sequently, a similar need was expressed over a band of 25 to 35 Me. 

A format of amplifier design was devised meeting these require- 
ments and subsequent constructions justified the design philosophy. 
While some of the techniques used are known, we believe that the 
nature of their use here emphasizes latent potentialities which had 
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not been generally apparent. In particular, a “feed-forward” error 
compensation system is used which has parallels to McMillan’s and 
Van Zelst’s generalization of Black’s error injection methods.1~* There 
is a very specific recognition here, however, of the physical role 
played by time, which dominates the form of the error network and 
which provides, among other advantages, a means of arbitrary error 
control of large time delay amplifiers which are intractable to feed- 
back methods.* 

We shall describe further the use of “corporate” amplifier binary 
arrays using quadrature couplers and emitter-follower amplifier stages. 
The word “corporate” is intended to suggest the process of a cor- 
porate coordination of many sources to deliver their powers in some 
desired collective mode, in much the same fashion as in a corporate 
antenna. With respect to the emitter-follower amplifier, much to the 
contrary usage of that device only as a unit gain isolation stage, 
there are many advantages highly germane to its use as an amplifier 
with certain basic advantages over the common emitter and common 
base modes of transistor operation. 

The name “feed-forward” error control scheme, coined to distin- 
guish it from feedback, is fully descriptive. Error is detected, ampli- 
fied, and injected after proper time delay in the forward time stream 
of the amplifier, and error cancellation is accomplished over the en- 
tire band of interest. The error, it should be emphasized, is com- 
posed of ripples in the amplifier transfer function, distortions through 
nonlinearity, and amplifier noise. To some extent the low level error 
amplifier has noise; nevertheless, there is advantage gained in noise 
performance by exchanging that noise for the much greater noise of 
the higher power amplifier. Further, the feed-forward system makes 
demands on the active elements of the system only within the band 
of interest, and great advantage is gained in using more powerful, 
lower quality transistors without the concerns of high frequency dis- 
persion instabilities attendant upon conventional feedback systems. 


II. CORPORATE POWER STAGES 


2.1 Coupler Arrays 


High-frequency transistors generally have a smaller dynamic range 
than the lower frequency and power transistors. One may, however, 


* During the late stages of preparing this manuscript, we learned of a patent 
issued by the German Federal Republic to Heinrich Schmidt-Briicken in which he 
recognizes time as a factor in error compensation. In Appendix 3 we compare his 
implementations with present implementations, and we believe the differences are 
significant. 
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resort to transistor arrays to multiply the power capability inherent 
in a single stage while retaining the high-frequency performance. The 
notion of an array of elements to provide power is well within the 
framework of conventional electronics as observed in the early ring 
oscillator, an antecedent of the existing distributed amplifiers and 
oscillators. 

In the ring oscillator, as in the magnetron, there is a coordination 
of the various elements into a collective mode. These modes comport 
with the rotational symmetry of the ring and so are characterized by 
eigenvalues which are the various rational roots of unity. Often the 
modes are virtually degenerate and mode strapping and mode killing 
techniques are a necessity to isolate the desired mode of operation. 
This isolation is frequently marginal and mode-hopping is a well- 
known consequence of improper excitation or loading of the system. 

The difficulty in mode suppression stems from the requirement that, 
for an n degree of freedom system, n — 1 modes must be suppressed 
without hampering the transmission qualities of the nth. Further, 
reactive strapping may operate, not so much to suppress, but simply 
to alter the dispersion characteristics of the network permitting the 
unwanted modes to pop out elsewhere and, often, with deleterious 
effects. 

Seidel has examined the problems of active network synchroniza- 
tion,* and the resolution of the problems of stability would appear to 
exist through the use of hybrid or conjugate couplers with absorbing 
suppressors. Basically, through the use of couplers, the terminating 
loads on the ideally unexcited ports provide exactly the n — 1 inde- 
pendent absorbers needed; furthermore, these absorbers are precisely 
where they should be to provide the necessary orthogonality against 
absorption of the desired mode. 

As opposed to the ring, the mode structure of the hybrid coupled 
amplifier is not characterized by a simple eigenvalue description. 
Further, the forms of building the coupler array may be quite vari- 
able and the choices reflect on considerations of some moment. There 
are two particular considerations which are of major consequence. 


(1) Out-of-band stability of the harnessed system 
(12) Low level of multiple reflection in interconnection of amplifier 
stages. 


While these two considerations are phrased independently, a small 
interdependence can exist. 
With respect to out-of-band stability, assume the possible use of a 
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transmission line hybrid such as, for example, a rat race or 3 dB in- 
phase coupler. While mecting proper in-band specifications on match 
or power, what are the consequences of performance of this harnessing 
structure upon the element whose activity persists well beyond this 
point? Further, with respect to the second consideration, what are the 
proper harnessing formats to permit arbitrary cascading of elements 
whose intrinsic isolations may, at best, be modest? 


2.2 Quadrature Couplers 


A solution which contains many practical advantages is the use of 
arrays of lumped quadrature couplers. A quadrature coupler, as an 
elementary device, is a matched reactive reciprocal four-port having 
two orthogonal axes of symmetry. This characterization suffices to 
endow the coupler with the following properties: 


(t) Corresponding to an excitation port, labelled 1, there exists a 
conjugate, or hybrid, port labelled 4. 

(w) The entire incident power emerges in quadrature phase from 
the remaining ports 2 and 3. 


An evident symmetry of these coupling characteristics obtains with 
respect to rotations about the two symmetry axes. The about results 
are demonstrated in Appendix 1. 

When we apply the device symmetries to the above labelling 
scheme we find the following scattering matrix: 


0 Si Sis 0 
Sex 0 0 Sis 
Sis 0 O Sipe 

0 Sis Sp 0 


S = (1) 


Application of the unitary properties of the scattering matrix, namely 
the relationship SS = 1, yields 


Si2S% + 8355 = 0 (2) 
S283 + 8,38,3 = 1. (3) 
From (2) S% = —S.s(8/S12) which we substitute in (8) to provide 
2 2 Sis 


2 — Sis = oy (4) 


12 
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Equations (2) and (4) may be put in a more useful form. Equation (2) 
provides that Re S,,.S* = 0 so that S,.S% is imaginary implying that 
S;. and S,; are in phase quadrature. Equation (4) leads to a result 
which we shall employ crucially: 


[St _ Sis| ad: (5) 


Equation (5) is of major importance, since it states that it is pos- 
sible to use quadrature couplers in an all-pass arrangement, as shown 
in Fig. 1, without regard to the frequency sensitivity of Sy. and S43. 
Since this result is independent of the frequency sensitivities of the 
matrix elements of S, it is possible to use the structural elements of 
Fig. 1 to implement far more complex power divider arrays without 
sharpening the frequency characteristics of transmission. Figure 2 
shows a four-way all-pass division and the extension of a 2” division 
is clearly evident. Figure 3 shows the harnessing of two 4-way dividers 
to form the 8-way divider. Shown, to, for the first time, is the place- 
ment of the elementary amplifiers. 

There is a key shown in Fig. 2 which is the algorithm for patterning 
the phase shifter arrays which introduce the relative 0° and 180° 
phase additions. Given any 2” array of phase shifts, the pattern for 
2"*t is given by forming the complement array and adjoining the lat- 
ter to the former. For example, if the 8-way array from top to bottom 
is (+——+—++-—), where + and — have substituted for 0° and 180°, 
respectively, the complement pattern is (—++—+——+) so that the 
16-way pattern is (+——+—++—; —++—+-—-—+), the 32-way is 
(+--+-4+4+-—-+4-4--4+; —-+4+-4+--4++4+--4+-++4+-), 
and so forth. 

Why quadrature coupler arrays instead of more conventional in- 
phase divider arrays? Why introduce the complexity in phasing the 
array? The essential answer is inherent in the simplest of the struc- 
tures shown in Fig. 1. The coupler has a 3 dB split with energy 
incident upon each of the amplifiers having 90° difference. Under 





Fig. 1— Quadrature couplers in all-pass arrangement. 
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reflection from identical amplifier inputs the phase difference mounts 
to 180°,.as-returned to the source, and a high degree of the reflection 
is cancelled. On the other hand, the reflected energies rcturn in phase 
to the port conjugate to the input and are dissipated in a dummy load. 

This technique, using 3-dB quadrature couplers to provide high 
input return loss to arrays of mismatched elements, was first used 
for short-slot couplers in radar transmit-receiver switches to prevent 


ALL~- PASS 





ALL- PASS WITH RELATIVE 
180° PHASE REVERSAL 


Fig. 2— All-pass four way quadrature power divider. 


mismatch back to the magnetron. While this technique has been ex- 
tended to tunnel diode and transistor amplifiers, nevertheless the 
specific arrays shown here differ in two basic respects: 


(t) The arrays can have arbitrarily large multiplicity. 

(iz) All-pass transmission results from the use of 180° phase dif- 
ference sections.* 

The arrays shown in Figs. 1 through 3 were developed through 
a graphical induction process. Their properties are, however, more 
clearly phrased by returning to an analytic description. In terms of 
the scattering coefficients, the transmission 7’ and input reflection K 
of the array are 


—_ T = (Sto = iat (6) 


* The use of 180° phase difference sections and their application to high multi- 
plicity dividers was independently conceived. Original recitation of their use in 
two section dividers was made by EK. A. Marcatilli, and D. H. Ring, IEEE 
Brien Groups Microwave Theory & Techniques MTT-10 #4 (July 1962) 
pp. 251-8. 
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Fig. 3 — Eight way quadrature divider with amplifiers. 


A 


and 
K = (Sis + Sis)"k (7) 


where ¢ and k are the transmission and reflection, respectively, of each 
amplifier element and where 2” is the number of amplifier elements. 


We observe from (5) and (6) that | 7 | = |¢|, confirming the all-pass 
construction. 
At the center frequency of each coupler, | Sy2| = | Sis| = 1/(2) 


and, via (2), there always exists phase quadrature. Therefore, Si, + S75 
vanishes at centerband and, whatever the multiplicity of that zero, (7) 
shows the reflection zero to be raised to nth degree. Simple coupler 
synthesis techniques permit an arbitrary multiplicity of the degree of 
the zero and no major difficulty exists in obtaining octave performance 
and much greater. 


2.3 Coupler Economics 

It is instructive to consider the simplest couplers meeting the ob- 
jectives of the amplifier to be constructed. The coupler enters into 
the technical design of the amplifier with respect to the return loss 
of the coupled amplifier array over the band together with a deteriora- 
tion of the dynamic range relative to a single amplifier element, On 
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the other hand, the choice of coupler enters into the economics as a 
greater cost for a more complex coupler. 

We shall find that the return loss is not the essential problem of 
the array but that it is the dynamic range which is limiting. While 
dynamic range deterioration does dominate the coupler choice, we 
shall see that the degree of deterioration is not major and that the 
feed-forward compensation may, in many instances, make the prob- 
lem negligible. Let us now investigate both the return loss and the 
dynamic range problems. 


2.3.1 Return Loss 


We shall consider two coupler types; one and two stage couplers. 
Corresponding to the number of stages is the multiplicity of zeros in 
the reflection factor, k, of (7) over the band. We intend to discuss 
these couplers in detail in another publication® and we shall consider 
it sufficient here to state simply that the couplers are compact and 
are composed primarily of lumped elements. They meet the charac- 
terizations of Section 2 to a high degree over extremely wide band- 
widths and all the properties described are met substantially in 
practice. 

We shall be concerned here only with the magnitude of k. Equation 
(7) yields 


| K | = || Sie? — | Sis |? [| I, (8) 

which produces the return loss L in db; 
L = 20n logio | | Sie es — | Sis | | + 20 logio | k . (9) 
Recognizing that |S;2|? + |Sis]* = 1, from (5), equation (9) yields 


Let us choose a normalized frequency variable, 9, such that 
Q = 1 at center band. We then assert for one and two stage couplers, 
respectively, that 


| Si. |? = 1/1 +? (single stage) (11) 
and 
| Sis |? = 1/1 + [1.5452 cos (20.6°2) — 1.1950? sin (20.6°0)}? 


(double stage) (12) 
where the coefficients in the two stage unit are chosen such that the 
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differ by the same amount at centerband and at the edges of a 1.5:1 
bandwidth. The behavior of (12) is shown in Fig. 4. 

For a specific choice of a 1.5:1 bandwidth, as was required for the 
60-90 Mc amplifier, we find 


2| Si, | — 1 S 0.22 (single stage) 
S 0.038 (double stage). 
If we assume that our manner of use of the amplifier elements as 
emitter followers produces |k| generally close to unity, we may 
neglect its inclusion in (10) and 
L = —13.16n dB (single stage) (18) 
—28.4n dB (double stage). (14) 


The major need for high return loss resides in the interaction 
between amplifiers. The dominant error produced by multiple reflec- 
tions is the “third time around” incidence which is down from the 
initially incident signal by the sum of the return losses of the output 
and input of the first and second amplifiers, respectively. Since the 
gain per stage exceeds 6 dB for emitter follower amplifier elements, 
each stage is capable of feeding a next stage composed of four times 
as many amplifier elements. Therefore, if the first stage contains 
2” elements the second contains 2"+?. If Z,12 = LZ, + Le in dB. Then, 
from (13) and (14) 


D ise = —26.32(n + 1) dB (single stage couplers) (15) 
= —56.8(n + 1) dB (double stage couplers). (16) 
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Fig. 4 — Two stage coupler designed for 1.5:1 bandwidth. 
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The fine structure in transmission ofr n = 1 or greater is very small 
and represents extremely small power variation. This presents a 
negligible taxation of the power handling capabilities of the feed- 
forward compensation system which will wash it out. Under these 
circumstances, the choice between single and double stage couplers 
does not rest upon return loss considerations. 


2.3.2. Saturation Power 


We may find from (6) a means for characterizing the energy paths 
through the coupler array. In much the same way as that in which 
the Pascal triangle builds up the binomial expansion in terms of the 
number of paths connecting two points, we may use the expansion 
of (6) to enumerate the number of paths from input to output. As 
an example, consider n = 3 which corresponds to an eight element 
amplifier stage. We have 


T= (Sts _ 38i3Sis se SSiSis = Sis} . (17) 
As may confirmed using Fig. 3 as a guide, we have 


(1) One path making six transitions of an S;. form 

(it) Three paths making four transitions of an S12. form and two 
of an S;3 form, together with a net 180° phase shift 

(ii) Three paths making two transitions of an S;. form and four 
of an S,3 form 

(tv) One path making six transitions of an S,;3; form together with 
a net 180° phase shift. 


As is evident from this enumeration, different paths contain some- 
what different amplitude content because of the frequency sensitivities 
of the scattering parameters Sj. and S,3. This suggests that the am- 
plifiers are somewhat unevenly excited which must have deleterious 
consequences to the net dynamic range of the composite amplifier. 

Equation (17) contains too much information in that it does not 
distinguish between paths prior to amplifier excitation, which would 
bear on saturation, and paths leading from the amplifiers which have 
only to do with collection, which do not bear on saturation. This 
dichotomy into exciting and collecting paths may be accomplished 
simply by factoring (6) into 


T = (Sie a I Si3)"(Si2 aa S43)". 
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We then define two structure factors; an excitation function F, and a 
collection function Ff. We have 


i = (Sie + Sy3)” (18) 
FP, = (Siz — Sis)”. (19) 


If in Fig. 3, the 180° phase structures had been placed prior to the 
individual amplifiers, F, and Ff, would have reversed in formulation. 

From the point of view of saturation we are only concerned with 
the form of the structure factor F,. Again considering an eight am- 
plifier array which corresponds to n = 3, we have 


PF, = - = 882513 — SSiais = Sis (20) 
which provides the following excitation paths 


(t) One threefold transition Sj. 

(71) Three twofold S,. transitions with one S,3 transition 
(wz) Three single So transitions with a twofold S,3 transition 
(iv) One threefold transition S,3. 


To form relative magnitudes, take the relative phase of S,. to be 0° 
and that of S;, to be 80°. We have the following values at the lower 
band-edge in a 1.5:1 frequency range. 


Single-Stage Douple-Stage 
Coupler Coupler 
Sie 0.78 0.720 
Sis z 0.625 z 0.694 


In Table I are tabulated the terms of (20) where the quantity 


1 

ah 0.354 

would represent the equal term values if the coupler were an ideal 
3 dB power splitter. It is evident that the amplitude variations 
shown in the F, column for the single stage coupler may be quite 
significant and, while they are more modest for the double stage 
couplers, they are, nevertheless, still present. Let us propose a harsh 
model of saturation in which any one amplifier saturates completely 
at the 0.854 level. In the SAT column of Table I are listed the relative 
amplifier outputs after the application of saturation with the assump- 
tion of no substantial phase shift effects. 
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TABLE I—JTxclTaTION AND Output TERM VALUES 
FOR SATURATION MopEL 











Value for Single Stage Value for Double Stace 
Term Multiplicity Coupler oupler 

F, SAT F, SAT 
S13 1 0.474 0.354 0.386 0.354 
S12?S13 3 20.381 10.354 20.3585 10.354 
S12813? 3 0.3045 0.3045 0.347 0.347 
S13 1 —210.244 —10.244 —210.334 —10.33 

1 

san 0.354 0.354 








If we recognize that the individual amplifiers are symmetrically 
disposed between input and output, as shown typically in Fig. 3, then 
the term value corresponding to the transmission from the input ter- 
minal to the input of any one amplifier corresponds, within a minus 
sign, to the term value for the transmission from that amplifier to the 
output terminal. Therefore, to form the transmission magnitude under 
this model of saturation, we sum the magnitudes, respectively, of the 
product of multiplicity, the term value, and the value listed in the 
saturation column, corresponding to each of the term values. Thus, 
using the single couplers we have 


T = t{0.474 X 0.354 + 3 X 0.381 X 0.354 +3 X (0.3045)” + (0.244)"} 
0.9095¢ (single stage coupler, 8 element amplifier), (21) 


while far the double stage couplers we have 
T = {0.386 X 0.354 + 3 X 0.358 X 0.354 + 3 X (0.347)? + (0.334)"} 
0.9894¢ (double stage coupler, 8 element amplifier). (22) 


Corresponding to the single stage, eight element amplifier there is 
a 0.88 dB distortion at the lower band edge while, corresponding to 
the double stage coupler amplifier there is a 0.08 dB distortion, where 
these distortions occur relative to an amplifier output at eight times 
the saturation power of an individual elementary amplifier. The 
implications of these distortions must be further elaborated, how- 
ever, in relating the choice of couplers on the economics of the 
amplifier. 

The intent in design, from the beginning, was the use of a feed- 
forward compensation system able to “clean up” amplitude distortions 
of 1 dB at maximum power, namely, the use of a clean-up amplifier 
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having about one-fourth the output power of the main amplifier. 
While the 0.88 dB, corresponding to the single stage coupler, is within 
the 1 dB clean-up capability, it does form a reasonably large propor- 
tion of it. For many applications, however, amplifier power covers 
the entire spectrum of the bandwidth and the proportion of the power 
associated with the band edge is very much less than full power. Un- 
der these circumstances, the clean-up capabilities are quite adequate. 
If it is required to provide single frequency performance anywhere 
within the band, then about 10 per cent power output reduction is 
advisable. 

The use of double stage couplers forms more of a fail-safe situation 
and the feed-forward amplifier is barely exercised. This leaves the 
feed-forward amplifier with its capabilities intact to handle phase and 
gain distortions without the excess power demand introduced by pre- 
mature saturation within the main amplifier. The cost, of course, is 
in the use of couplers having roughly double the expense. 

A further basis of choice of the double coupler lies in the use of the 
amplifier without the feed-forward implementation. According to the 
tabulation for the eight element amplifiers the term with the highest 
degree of saturation is S?, which has an unsaturated value of 0.386 
for the two-stage coupler as opposed to the saturation level of 0.354, 
corresponding to a 0.72 dB saturation. Without feed-forward, this loss 
of capacity of the full amplifier, namely 0.72 dB, would be required 
to avoid any of the aspects of saturation. The corresponding penalty 
paid using the single stage coupler is 2.52 dB, which might be too high. 


2.4 Hybrid Power Divider Arrays 


We have discussed power divider arrays using only quadrature 
hybrids. It might be argued that it is required, at most, to be con- 
cerned only with paired amplifiers since this is all that is required 

to form adequate in-band match. After having performed such a 
- match then an arbitrary number of such pairs might be combined 
using in-phase, as opposed to quadrature-phase, power dividers. 

There are several objections. We list and discuss them in order. 

Stability: Quadrature couplers have excellent match characteristics 
from de to, typically, greater than three times their center frequency 
of operation. Since match and isolation are directly related, the am- 
plifier elements of a large quadrature coupler array stage generally 
do not interact over the active range of these elements, assuring sta- 
bility. The control of parasitics in an in-phase divider, particularly 
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in a high frequency ferrite core type of divider, is not of the same 
degree of excellence and interactions leading to oscillation are a more 
severe problem. 

Match: Assume that the array were composed of amplifier element 
pairs coupled via quadrature couplers, with the remainder of the 
array being in-phase dividers. The return loss from (13) and (14) 
would be 13.16 dB at the lower portion of a 1.5:1 band for a single 
stage coupler and 28.4 dB for a double stage, assuming the amplifier 
to be a perfect mismatch. The economic advantage of the single stage 
coupler is so major that, assuming an insistence on in-phase dividers, 
and all other considerations aside, a minimum single-stage quadrature 
coupling of four elements would be most desirable. This would lead 
to a 26.32 dB return loss, which is generally acceptable. 

Cost: By far the most major advantage of the lumped quadrature 
coupler is its extremely low cost, particularly as a single stage device. 
It is reasonable to anticipate a cost of a few cents, as opposed to 
several dollars for other types of hybrids. In the vhf range an in- 
phase divider is composed of a pair, or more, of multiply-wound fer- 
rite cores (Ruthroff transformers). Microwave distributed couplers 
are vastly more complex. A hybrid array of in-phase and lumped 
quadrature couplers is far more costly than an array of quadrature 
couplers only. 

Power Handling: Amplifier arrays under development are yielding 
typical average powers between 50 and 100 watts. It takes no great 
leap of imagination to conjecture such amplifiers in the kilowatt to 
multikilowatt range for any of assorted operations. Should any severe 
imbalances occur within the array owing to the failure of a transistor, 
an in-phase divider in the high power portion of the amplifier would 
experience large internally circulating power. If the divider contains 
high frequency ferrite transformers, burnout is virtually assured 
because of ferrite heating. Lumped quadrature couplers are air-core 
devices capable of tolerating kilowatt power levels without burnout; 
the only limitation is the dielectric strength of the insulation, which 
can be great. 


2.5 A Two-Rail Array 


While the authors believe the exclusively lumped quadrature cou- 
pler array to have prominent advantages over other forms of organ- 
ization, it is of interest, nevertheless, to show one specific situation 
in which a hybrid system of in-phase and quadrature couplers has 
certain advantages in construction. 
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Consider a situation in which ferrite core implementations are not 
realistic to provide the 180° phase shifts required for all-pass opera- 
tion. Such a situation occurs in the upper portion of the uhf region, 
and beyond, where the ferrite magnetization can no longer follow a 
rotation mode. Both in-phase dividers and differential 180° sections 
become much more expensive and cumbersome, but let us assume that 
the expenditure is warranted in providing two each of these afore- 
mentioned devices. 

Figure 5 uses the in-phase dividers in a fashion reminiscent of two 
rail logic. Two universes are formed 180° apart via the in-phase 
dividers and the differential 180° sections. Whenever a 180° phase 
shift is called for it is accomplished simply by a symmetric inter- 
change or transposition between the two universes. Thus, the cost 
of a large multiplicity of 180° sections is reduced to that of provid- 
ing only two, aside from the two in-phase dividers. 

Another advantage is gained in respect to saturation. Figure 5 is, 
in actuality, a doubled four-way quadrature divider to form an eight 
element amplifier and, as such, has a lower degree of compression 


IN 
PHASE 
POWER 

DIVIDER 





Tig. 5— Hybrid array of in-phase and quadrature couplers: doubled for way 
divider to form eight element amplifier. 
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than does a pure quadrature coupler array. Corresponding to (18) 
we have for n = 2 


fr, = (Sie a Sig) = Sto + 2828.3 + Se (23) 


Let us now consider only the single stage coupler for which we have 
tabulation shown in Table II. Under the model of saturation proposed 
earlier we have 


—T = t{0.5 + 2 X (0.487)? + (0.397)*} 
£{0.9362}. (24) 


The transmission value of 0.9362¢, or 0.48 dB, in the double-four 
amplifiers compares to the 0.88 dB compression relative to capacity 
of the full eight amplifier array. There is less demand on the feed- 
forward cleanup and power reduction is, therefore, unnecessary. Sin- 
gle stage couplers, therefore, produce acceptable performance and 
the resulting economy offsets the in-phase divider costs. 

Having discussed the limitations of hybrid arrays, is Fig. 5 a 
workable system? The two major questions one may pose are those 
of power and stability. The experimental sections of this paper dis- 
close that such a structure was successfully constructed for the 75 
Mc amplifier. The ferrite loaded in-phase dividers were adequately 
broad, extending beyond the range of significant transistor activity, 
and stability was not the problem. However, the means to provide 
power handling capability within the ferrite transformers was a 
problem. 

Assume that the amplifier were to be designed, not on a fail-safe 
basis, but on the assumption that all stages operate appropriately. 
Since the in-phase power divider (Ruthroff hybrid) has a conjugate 
antisymmetric port, can this port be used to save the need for the 
180° differential phase shifters? The answer is negative in that the 
antisymmetry port would require power transformation through 
the ferrite. Since the core must be small to avoid circuit parasitics, the 


I 


TABLE II—SINGLE STAGE COUPLER 





Term Multiplicity Value F. SAT 
Si? 1 0.608 0.5 
S12813 2 10.487 10.487 
2 


Si3 1 —0.397 —0.397 


ERROR-CONTROLLED AMPLIFIERS 667 


handling of large continuous power might cause burnout. The sym- 
metry port barely excites the ferrite and its use is indicated. 

The circuitry of a simple 1:1 transformer is much less complex 
than that of a Ruthroff hybrid and the 180° section admits construc- 
tion as a heavy core ferrite reversing transformer. In the event that 
the power might yet be too high to sustain in the reversing trans- 
former, simple nonferrite circuits are available which produce 180° 
phase differentials over a 3:1 bandwidth. 

As stated, this system does not contain the reliability feature dis- 
cussed earlier. Transistor failure produces severe imbalance in the 
array and ferrite excitation in the transformers would result, with 
possibly deleterious effects at high average power. The specific am- 
plifier design using the transformers with which the authors were con- 
cerned, related to pulsed operation, and this failure possibility was not 
significant. 


2.6 Out-of-Band Filters 


We have, up to now, considered in-band performance in which the 
return losses of the various stages were adequate to suppress multiple 
reflection effects. Out of band, the power split of the couplers becomes 
highly unequal and the input reflection into a quadrature-coupled pair 
of amplifiers approaches the amplifier reflection. 

A simple fix for this unstable multiple interactions situation is 
obtained through the use of band-pass multiplex filters. In band, these 
filters provide excellent transmissions between various stages of the 
filter; out of band, all energy is transmitted reflectionlessly into 
dummy loads, with high isolation between the previously well coupled 
ports. 

The multiplex, or out-of-band filters, are themselves quadrature 
couplers having a desired synthesis of |Srel. The design through 
quadrature coupler techniques® is made to assure the wideband match 
required for stability. In the construction of the amplifier system, 
the use of these filters was determined on a practical basis. If stages 
directly feeding one another tended to oscillate, an out-of-band filter 
was injected between them and a large margin of stability was 
restored in each case. 


2.7 Comparison with Distributed Ampliflers 


The distributed amplifier 1s a chain of active elements in contrast 
to the parallel-like arrangement of the binary complex tree. Since 
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the excitation varies along the chain one must reduce the gain of 
each successive stage to preserve equal power output per stage. The 
implicit assumption here is that each amplifier stage has been de- 
signed for power output to the limit of its intrinsic power capability. 
The gains of the various stages in sequence are clearly interdependent. 

There are four major advantages of the coupler array over the 
distributed chain: 


(1) The gain profile along the chain in a distributed amplifier is 
tailored to accommodate to a presumed performance of each stage, 
and the net energizing of any one stage in the chain is a function of 
all those preceding it. Fluctuations of performance in time, or simple 
tolerance effects, are cumulative. This is not true for the coupler 
array in which each element is orthogonally energized. 

(it) Imperfections along a chain produce multiple internal scat- 
tering with a necessary result of rapid gain variation with frequency. 
In a signal flow sense, each active element of the coupler array is 
“hit” only once with substantially no multiple reflection. 

(421) In the event of the failure of a single element of a distributed 
amplifier chain, the transmission properties are greatly modified, pro- 
ducing no fail-safe guarantees. Through the orthogonality properties 
of the coupler array construction, the failure of any one active ele- 
ment component simply causes its ordinarily contributed voltage to 
be subtracted from the common mode without any catastrophic inter- 
actions with the other properly functioning components. Thus, the 
coupler array provides a fail-safe redundancy. 

(zv) The return loss properties of a quadrature coupled array im- 
proves with complexity. The return loss of a chain, contrariwise, 
deteriorates with complexity. 


The intrinsic return loss property of a quadrature coupled array 
improves with increasing complexity through the presence of an ever 
greater number of parasitic-mode absorbers. A chain composed of 
imperfectly isolated stages becomes ever more difficult to match with 
an increasing number of stages. 


2.8 Dynamic Range 

It is evident that coupler distributed systems provide increased 
power in proportion to the number of active elements. Notice that 
this increased power output comes with little expense to the system 
noise temperature over that of any one of the active elements, aside 
from the copper loss of the couplers. 
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The noise emission of any one amplifier is spontaneous and, as 
such, is entirely uncorrelated with that of any other active element. 
The coupler system is a 2" = N mode system and, with this lack of 
internal correlation, each active element as a 6 function source equally 
stimulates each mode. Of the N modes, N — 1 are absorbed in the 
dummy loads, and only 1/N of this power finds egress. With N indi- 
vidual sources, the output noise emission is N X 1/N of the noise 
power of any one element, and so, is just the noise emission of one 
of the elements. 

The dynamic range of the system is given by the ratio of the power 
at saturation to the noise power emitted over the amplifier bandwidth. 
Since the upper limit of the power has increased by a factor of N, and 
the noise level is unaffected, the dynamic range also increase by a 
factor of N. It should be recognized that the dummy loads of the 
coupler system occupy only ports orthogonal to the power transmis- 
sion path, hence, none of these loads contributes to noise. It is con- 
cluded, therefore, that the coupler organized power multiplex in- 
creases dynamic range with insignificant noise figure deterioration. 


III. EMITTER FOLLOWER AMPLIFIERS 


3.1 General Considerations 


The desire to deliver large power requires, in turn, a choice of an 
active element compatible with that power, together with a capability 
of supplying the circuit demands for proper frequency behavior. High 
power transistors always exhibit large internal capacitances and are 
limited in bandwidth. 

The major problems arising in a general design that uses power 
transistors are the circuit complexities to realize the desired fre- 
quency characteristics, and the reproducibility from device to device 
for large variations of the internal parameters. The resolution of 
these problems may be carried out effectively by means of the emitter- 
follower amplifier. 

The emitter follower, the analogue to the cathode follower, is gen- 
erally used as a buffering device in transforming, at substantially a 
unit voltage gain, between a high impedance source and a low im- 
pedance load. While it is well recognized that power gain exists, a 
means for obtaining voltage gain between equal impedances has been 
viewed as demanding too complex a use of transformers in respect to 
parasitic elements. The eminent desirability of the emitter follower 
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comes, of course, from its large self-degeneration which makes its 
response fairly impervious to internal parameter variations. 


3.2 Marriage of Emitter Followers and Quadrature Couplers 


Let us review some pertinent features of the ideal emitter follower, 
where the connotation of ideal implies an infinite input impedance, 
a unit voltage gain, and very low output impedance. Figure 6 shows 
the insertion of an ideal emitter follower between a normalized 
Thevenin generator and an equal impedance load. Since the open 
circuit voltage is twice that delivered to a matched load, the generator 
is characterized by a voltage of amplitude 2. The insertion of the emit- 
ter follower, therefore, doubles the voltage ordinarily delivered by the 
generator to the load, and a 6 dB gain results. 





Fig. 6 — Thevenin equivalent of insertion of ideal emitter follower. 


This gain, however, does not lend itself to a simple transformerless 
cascade of emitter followers as seen in Fig. 7.* It is evident that hav- 
ing achieved a voltage level of 2, this amplitude is transferred from 
section to section without further enhancement. The fallacy in form- 
ing this simplistic cascade is that match has not been restored from 
stage to stage, remembering that it was precisely through the inter- 
face of a matched source and open circuit that the 6 dB gain was 
achieved. 

This situation is rectified by the use of quadrature couplers to pair 
emitter followers as shown in Fig. 8. The harnessing structure of the 
couplers forms an all-pass unity transmission network and the gain 
takes place via the mechanism of Fig. 6. Since the couplers do restore 
match, substantially, at the input as well as the output, this circuit is 
capable of iteration with return losses as given by (7). Specifically, 
incident waves of amplitude S,. and S,3; impinge, respectively, on each 
of the emitter followers. Since each emitter follower is an open circuit, 
a voltage reflection factor of unity doubles each of these incident 

* The lack of a transformer distinguishes this cascade from a Darlington net- 


work in which cumulative power gain occurs only because of impedance trans- 
formation. 
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Tig. 7— Transformerless cascade of emitter followers. 


voltages, forming emergent waves of 2S,,. and 28,3, respectively. The 
termination is excited by an amount 2(S{, — Sj.) which, from (5), 
demonstrates a 6 dB gain. 

Since the circuit of Fig. 8 shows an inheent 6 dB gain, no trans- 
formers are required to achieve this modest, but highly useful, gain. 
However, a mild input or output voltage transformation might be 
used to raise it somewhat beyond this value. From the point of view 
of power generation, the stages shown in Figs. 1 though 3 might in- 
crease at a four-fold rate. Such proliferation of circuit complexity 
seems not nearly so modest. 


3.3 Real vs. Ideal Emitter Followers 


The ideal emitter follower is a most imperfect representation of 
the real device. However, it is not too difficult to patch things up to 
obtain a creditable performance within the frequency range of our 
interest. Figure 9 is a simplified high frequency equivalent from 
which we may deduce the more significant features of operation. 

Define the following quantities: 


Y=1+imw0,=1+i7% 
1 
1 





Co = 


Wep 


{25:2 


ee. 


QUAD QUAD 
COUPLER COUPLER 











Fig. 8— Transformerless pairing of emitter followers through quadrature couplers. 
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Fig. 9 — Simplified high frequency equivalent of terminated emitter follower. 


E, = Base signal voltage 
E, = Output voltage at the emitter 
ae” 
Ca SS 
1+i> 
De 


7 = emitter to base drift time and is of order 0.6/w, . Analysis 
yields the results 


ik, = 1 ef e + rl — a) v} i, (25) 


E (1 oF 4s) ty dk =a) 
uA Sef eee res IN : Web 7 
an # iwC., + Y(1 — a) 
Since we are seeking nominally small gain we may, with good approxima- 
tion, take a = 1. Expanding, and dropping terms beyond w”, (25) 
becomes 


kh, = i F raf — cos wr — — (4 + 7) | 
Wr \We 
+ tw E ++ ro 3 -)} E,. 
Wes C2 


The quantity 7,, is small relative to unit loading and we may omit 
the term r,,(1/w, + 7) in contrast to 1/w,,. We have then 


B, = #,{[1-o(+ +2-4)| +24 (27) 
WoW, Wy 2 Wes 


We perceive in (27) a negative real part resulting from the Miller 
effect action with respect to C,. The tuning procedure demands flat 
power gain, and we may offset the decrease of the real part in (27) 


(26) 
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against the rising imaginary portion to maintain constant magnitude. 
Assuming small parasitics at the frequency of operation, the constant 
magnitude condition yields 


Vn fe oe ae 
—3- = Wp, alr wees 2 (28) 
Web WoW] Wy) a 
A value of w, always exists to satisfy (28). 

According to (27) and (28), the voltage appearing at the base emerges 
with unit gain from the emitter. The question, indeed, is what voltage 
appears at the base? To determine it, we examine the input impedance 
function given in (26) to within first order terms in w. We find 


1 . 
; 1 
ia Cs +—+ r) 
Wa 








Z, =T5 + (29) 


While (29) has a disconcerting form with respect to dimensions, as 
does (28), that a unit terminating impedance is dimensionally concealed 
in these results. Under these circumstances, 7;, is small compared with 
unity and we momentarily neglect it in comparison with the series 
capacitance C’,, + 1/w, + 7. 

Assuming this capacity not to be excessively large, it may be in- 
corporated as a terminating element of a bandpass filter that has a 
characteristic impedance equal to that of the generator. Over its 
bandwidth the generator incidence is maintained to within the time 
delay of the filter circuit and all the appurtenances of the ideal emit- 
ter follower are realized; an incident signal is delivered at matched 
impedance from the generator to the open circuit, and the voltage at 
the input is transmitted to the output in substantially frequency- 
independent fashion. 

In general, this bandwidth is far larger than the 1.5:1 frequency 
width desired and, corresponding to this gain-bandwidth availability, 
tradeoff may be made with the bandwidth to enhance the gain by 
using appropriate input circuit design. This circuit may also com- 
pensate for 7,, by providing high frequency peaking to offset it. 


3.4 Linearity 


While much has been said here of the linear behavior of the emit- 
ter follower, the feedback which causes the emitter to follow the base 
extends well beyond the region of linear characterization. In large 
signal excitation the interelectrode capacitances vary throughout the 
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cycle, causing nonlinear interactions with the circuit environment. 
A particularly significant term is C.,, but its variation is greatly 
damped by the generator impedance. Aside from the considerations 
of frequency characteristics, it is the desire to maintain high gen- 
erator damping that vitiates a greater use of impedance transforma- 
tion to yield higher gain. 

In class A operation the region of linearity is limited, in any event, 
by current cutoff. In class AB push-pull operation the distortion takes 
the form of a slightly different low and high signal gain. This last 
distortion may well be the more insidious of the two and it is only 
the use of the feed-forward system which makes it acceptable. This 
acceptability is extremely significant because the push-pull class AB 
operation, when operated near the class B state, has several advan- 
tages over class A: 


(1) Power supplies are substantially reduced. 
(it) Since all of the transistor current is used to produce RF power, 
_and since twice the number of transistors are used, about four times 
as much power is generated for the same amount of coupler hardware. 
(ww) Lower average currents yield wider band response because of 
lower average internal capacitance. In actual measurement the limit 
of the linear range seems closely identified with avalanche break- 
down, for which there is the danger of thermal damage to the tran- 
sistor. 


There is a notable disadvantage in class AB operation, but it is 
not too objectionable in a 50 per cent or less duty cycle amplifier. 
As the excitation varies, the average current varies, shifting bias 
points. Further, the operating temperature may shift the quiescent 
point in each transistor, significantly affecting low signal distortion. 
However, an active clamp may be used which samples the off-signal 
voltages and restores proper operating levels. This clamp is described 
in Section 6.8. 


3.5 Economics of Coupler Organized Emitter Follower Amplifiers 


While it is implicit throughout the technical descriptions, one 
general feature which must be stressed is the significant economy 
possible through emitter followers. Because of its large internal nega- 
tive feedback, reproduction of the transistor becomes a problem of 
relatively negligible proportions. 

In contrast to the situation in which refined transistor devices may 
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be improved in yield, we have a situation in which frankly cheap — 
devices are used to provide an adequate yield. The low potential 
costs of quadrature couplers and relatively inexpensive transistors, 
operating well beyond their traditional ranges, forms the prospect of 
much less expensive, high performance amplifiers. 

Acceptable within this framework may be a moderate aging of 
semiconductor parameters, usually a bugaboo in device manufacture. 
Its importance is diminished not only through the degeneration of- 
fered by the emitter follower circuit, but by the feed-forward provi- 
sions as well. Thus, cheaper encapsultations might be permitted 
than otherwise would be, admitting possibly consequential cost re- 
ductions in transistor fabrication. 


IV. FEED-FORWARD 


4.1 Feedback vs. Feed-Forward 


Much of the development of the high power amplifier stands or 
falls on the ability to devise a proper error correcting system. As 
stated in the introduction, the choice made in the design of this 
system, if indeed a choice was open, was a “feed-forward” construc- 
tion. 

In the late 1920’s, H. 8. Black conceived of feed-forward as well 
as feedback, and his patent disclosures show either grid or plate 
injection of the error signal.? The grid meant back injection whereas 
the plate meant forward injection and either was adequate to within 
180° of phase. Clearly, the plate injection was at a higher level and 
required subsidiary amplification, while grid injection ingeniously 
reused the same amplifier. 

With an occasional, but short lived, resurrection, feed-forward fell 
into a general disuse until independent revivals by McMillan and 
Van Zelst.1:? These were followed by further activities by Deighton, 
Cooke-Yarborough, Miller, and Golembeski.* 7 The essence of the 
McMillan and Van Zelst disclosures was to combine both techniques 
of Black by providing frequency shaping networks simultaneously in 
feedback as well as in feed-forward. The aim was to provide channel 
redundancy as well as parameter desensitizing. The Deighton paper 
uses this technique to stabilize a particular circuit of a pulse amplifier 
while Golembeski has represented a generalized parallel channel con- 
figuration and developed a synthesis procedure and realizable condi- 
tions, 
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Even though the basic inadequacy of feedback techniques were evi- 
dent, we developed the feed-forward compensation system because 
we needed to provide stable correction to an amplifier that had a rela- 
tively large time delay. Indeed, simple physical insights not only 
repudiate the use of feedback but force use of the alternative, for our 
case. 

Feedback, a technique which has been used with much success, 
attempts a causal contradiction: after an event has occurred, reshape 
its cause. This violation may be resolved only by time smearing the 
event to blur the distinction between “before” and “after” to an 
adequate degree. This smearing requires a prescribed control of spec- 
tral response; a craft too well known to presume further comment. 

Feed-forward evolves from entirely different premises, avoiding 
causal anomaly.’ Instead of seeking to reverse time, it recognizes the 
passage of time. Error is determined in relation to a time-shifted 
reference and corrected in a time sequence compatible with the main 
signal. No obscuration of time is permitted and frequency shaping 
to attain stability is irrelevant. 

The fact that frequency shaping (that is, control of the return 
ratio function) is irrelevant is of major importance. A feedback am- 
plifier consumes a good portion of its gain-bandwidth figure of merit 
in controlling the return ratio well beyond the band of interest to 
prevent instability. 

The power transistors demanded for economical construction of the 
power amplifier clearly are not geared to high frequency application. 
To obtain the needed bandwidths, gains must be kept modest. The 
amplifier, then, would be built of a sufficiently large number of stages 
to offset the gain consumed in a feedback correction. With this num- 
ber of stages, however, the excess phase would mount to such degree 
as to be beyond control in a feedback system. 

The feed-forward system has three strikingly important features: 


(1) Feed-forward correction does not reduce amplifier gain. 

(iw) Gain-bandwidth is entirely consumed within the band of in- 
terest. 

(wt) Feed-forward realization is independent of the magnitude or 
shape of the amplifier time delay and it exists in general. 


We demonstrate these features in the following section. The major 
point to be made here, however, is that it is not a question of feedback 
vs. feed-forward as a choice, but as a necessary engineering decision. 
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4.2 Detailed Description of Feed-Forward Operation 


Figure 10 is the general block diagram of the feed-forward cor- 
rection system. Input signal is divided into two portions; that going 
into the amplifier and that entering a reference path. The signal is 
amplified through the main amplifier A and a sampling is taken 
through a directional coupler with very flat coupling 1/G, where G 
approximates the in-band gain of A. This sample is compared with 
the reference path signal which has been delayed a time +, a quantity 
approximately equal to the time delay of the main amplifier. 

A time-shifted error is thus obtained and it is amplified by a sub- 
sidiary amplifier whose gain G’ is such that this error is restored to 
the appropriate level to cancel the error of the main amplifier. If 
the subsidiary amplifier consumes a time roughly equal to 7’, a time 
delay of that value is inserted in the main amplifier path. The injec- 
tion of the error back into the main line takes place through a reactive 
two-port, shown typically as a transformer. 

The choice of a reactive two-port, as opposed to the use of a coupler 
four-port, is made irrespective of the ordinary desire of keeping the 
outputs of both amplifiers decoupled. While a coupler would isolate 
the subsidiary amplifier from the main amplifier, the cost of energy 
losses of the subsidiary amplifier into the dummy loading of a coupler 
would be unacceptable. 

While we have emphasized the care to be taken with time delay, 
there is an implicit condition within the phase specification to produce 
error cancellation. Since time delay is the slope of phase with respect 
to frequency, assuming small phase curvature, it is adequate to adjust 
phase at but a single point in the band and have it track throughout. 
In the event that the main amplifier and the subsidiary amplifier had 
large phase curvatures, then 7 and 7’, respectively, would require more 
complex realization than simple delay lines; they would have to incor- 
porate the same curvatures as well. Evidently, this may be extended 
to yet higher phase derivatives with respect to frequency. 





Fig. 10— One loop feed-forward network. 


678 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


As indicated, the gain and phase of A and the gain and phase of 
G’, are known only to within some first order tolerance. Via the cou- 
pler 1/G and the time delay 7, strict gain and phase standards are 
set with which to compare the amplifier A. These form a Procrustean 
bed such that anything within A that does not fit these standards is to 
be treated as error, and chopped off or stretched. In processing just 
this error portion through G’ and 7’, first order errors of the latter 
two quantities only produce negligible second order quantities in 
error correction. 

Having “feed-forward compensated” the amplifier A in this man- 
ner, it is quite proper to consider this circuit complex as a new 
amplifier, B, which is to be “feed-forward compensated.” Since B 
has desired performance to within second order errors, this new com- 
pensation would reduce the errors to third order errors. This method 
of successive compensation may be iteratively compounded, control- 
ling errors to a very high degree. A two-loop system is shown in 
Fig. 11. 


4.3 Transformer Design for Error Injection Circuit 


As already discussed, the use of hybrid couplers for error injection 
is avoided, since any such matched device leads to excessive power 
loss to the system. If, for example, the coupler provided equal power 
division, then a 3 dB power loss would be sustained by the main 
amplifier. If a division ratio more favorable to the main amplifier 
were used, the consequences to the subsidiary amplifier might well 
prove disastrous. It could be called on to deliver a greater power 
than the main amplifier, in some instances, to compensate moderate 
errors in that main amplifier. 

- We will try to design a feed-forward compensation system that can 





Fig. 11— Two loop feed-forward network. 
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correct the main amplifier over a range of about 1 dB, using trans- 
former injection. 
Appendix 2 shows that the optimum transformer ratio is 


| 


N= wae 


(30) 
under the condition that the subsidiary amplifier sees a matched load 
when it is correcting the peak relative in-phase voltage error, 8,. At 
this point the power being delivered by the subsidiary amplifier is 


bp 


si aes ESE 





as (31) 


where P,, is the main amplifier peak power. It is shown that the 
mismatch introduced by the transformer results in no loss to the 
power output of the system, only a small loss of gain, so there appears 
none of the loss that would be associated with a hybrid coupler instru- 
mentation. In our case of 1 dB correction, §, = Y%, for which N? = 7. 

The analysis shows another interesting feature. If, with no feed- 
forward circuitry present, a relative distortion occurs at peak power 
approximately equal to 26,, then: 


(t) In the presence of feed-forward circuitry the power gain is 
depressed by a relative amount $,, so that the power necessary to 
produce an undistorted output is only 8). 

(w) This residual distortion power of 8, is contributed by the 
subsidiary amplifier. 

(viz) In this sense, a subsidiary amplifier power capability of 8, 
controls a distortion of 28). 


The above design takes into account the in-phase errors, mainly 
in-phase compression in the amplifier at peak power. Another considera- 
tion is the ability of the above design to compensate for phase distortion. 
It is shown that, if the maximum relative quadrature voltage at peak 
main amplifier power is 76 , the relation between this and the maximum 
permissible in-phase voltage error is 


of = 


Pp 


dle 


is (32) 


Then the same system which will correct a 1 dB amplitude error with 
6, = § will provide a phase correction of + arctan 6f = +3.6° at 
maximum main amplifier power. . 


The above might seem a small phase correction range but remem- 
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ber that the main amplifier consists of an array of emitter follower 
amplifiers. To within small variable transit time effects with high 
power, the emitter follower amplifier would tend to saturate as a func- 
tion of instantaneous excitation. Hence, satuation effects would tend to 
be in phase and would not seriously tax the range of quadrature cor- 
rection capability. Nevertheless, from the equalization point of view, 
it would be well to limit phase ripple peaks to below 2° before apply- 
ing feed-forward correction. 

We emphasize that the numerical results obtained above correspond 
to the worst case and should not be used inflexibly. For instance, there 
might be some best phase positioning of the subsidiary amplifier along 
its transmission path to optimize the phase correction range. Or suit- 
able judgment might suggest an optimum trade-off between in-phase 
and quadrature error correction in a specific construction. 

Appendix: 2 also considers two questions concerning match within 
the system. First it shows that, even though the subsidiary amplifier 
is mismatched at less than maximum excitation, reflections within 
the error correction loop can never be phased so as to demand either 
more current or voltage from the subsidiary amplifier than is de- 
manded at peak power. Second, it considers the output mismatch and 
the use of nonoptimal transformers. 

The term nonoptimal is meant to suggest that the transformer is 
not optimally designed to minimize the power demand on the sub- 
sidiary amplifier. As shown, a non-optimal transformer yields the 
advantage of simplifying the problems in providing very high output 
return loss from the amplifier system, but at major cost to the sub- 
sidiary amplifier. 


4.4 Circuit Modification for Optimal Noise Performance 


As discussed, the dynamic range of a coupler organized amplifier 
is raised at the saturation power level in proportion to the number 
of active elements without deteriorating the noise level of the system. 
This promises a dynamic range increase of only a factor of N. It is 
possible to obtain much more spectacular ratios by a noise exchange 
with a low noise subsidiary amplifier. This improvement is of major 
significance since the active elements are individually high power 
devices emitting large noise powers and advantage is to be gained, 
in providing dynamic range, by working at the low end of that range. 

We emphasize that the subsidiary amplifier in the feed-forward 
system has an entirely different function in noise reduction than does 
an amplifier used as a low noise preamplifier, In the latter use, the 
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amplifier must sample the entire signal range and, intrinsically, limits 
dynamic range. The subsidiary amplifier, on the other hand, handles 
only the spontaneous error and, consequently, a negligible portion of 
the coherent signal range. Therefore, the dynamic range of the sub- 
sidiary amplifier may have virtually no relationship to the dynamic 
range of the system it controls. 


4.4.1. Major Considerations 


The feed-forward or subsidiary amplifier has a considerably lower 
dynamic range than does the main amplifier. With multiple feed- 
forward looping, as indicated in Section 4.2, the dynamic range of 
the ultimately controlling amplifier may be extremely small, yet 
adequate. The low power requirement of the subsidiary amplifier in 
such use may be turned to advantage by designing for relatively low 
noise figure performance. 

The noise of the main amplifier is an erroneous output which the 
feed-forward system corrects. However, noise generated by this sub- 
sidiary amplifier is introduced in the process. The advantage of such 
an exchange would be lost if the signal-to-noise ratio of the subsidiary 
amplifier were to be lessened by the use of significant power division 
at the input to the reference path. 

Figure 12 shows a modification of Fig. 10 to achieve optimal noise 
performance. Couplers C,; and C3 are made to direct most of their 
power to the subsidiary amplifier G’. While the gain of amplifier A 
must be made larger to offset the coupling loss of C,, this increase is, 
in principle, irrelevant. The worsened noise figure of A is of no con- 
sequence since the error cancellation of the feed-forward system sim- 
ply treats the additional noise as additional error. 


4.4.2. Nonideal Behavior 


Of course, what we have just stated is an ideal view; the cancellaticn 
is imperfect. The primary source of this imperfection lies in the varia- 
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Fig. 12— Modification of feed-forward network to improve noise performance. 
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tions of the subsidiary amplifier, which is uncontrolled. A realistic 
quantitative assumption might be a power transfer variation of +5 
per cent over the frequency range or a peak voltage transfer variation 
of 24 percent. This latter number limits the cancellation to a power 
ratio of (0.025) or 32 dB. Let us assume further that C, is a 12 dB 
coupler, increasing the excess noise of the main amplifier by that amount. 
Nevertheless, even with a reduction of the 32 dB cancellation possibility 
by 12 dB, the net noise of the feed-forward compensated system is the 
sum of the excess noise of the subsidiary amplifier plus a 20 dB reduced 
excess noise of the main amplifier. Since the subsidiary amplifier noise 
is presumed minor, relative to that of the main amplifier, this represents 
a major gain in performance. We have neglected the small transmission 
loss of C, which deteriorates the noise performance slightly. 

In multilooping along the lines of Fig. 11, successive reductions of 
the above order may be accomplished until the noise of the lowest 
excess noise subsidiary amplifier stage dominates the process. Because 
of the rapidity of convergence it is questionable whether more than 
three subsidiary loops would ever be meaningful with respect to noise 
minimization. 


4.4.3. Subsidiary Amplifier Stabilization 

Since the limit of noise cancellation of the high power amplifier 
relates to gain variations of the subsidiary amplifier, we suggest an- 
other form of multilooping, differing from Fig. 9, in which the subsid- 
lary amplifier itself is feed-forward stabilized. We believe that this 
variant, shown in Fig. 18 is inferior to that of Fig. 11 in that a much 
more restrictive tolerance is placed on couplers Cz and C3. This is so 
because the error determination is now only “one shot” with reference 
to the main amplifier as opposed to the progressive error sampling of 
the earlier system. 







MATCHING 
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Fig. 18— Modification of low noise system in which subsidiary amplifier is 
stabilized. 
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Nevertheless, the looping progression shown in Fig. 13 requires 
that some weak subsidiary amplifier only be confronted with a next 
stage merely one degree higher in power. In contrast, Fig. 11 requires 
that all subsidiary stages work into the highest power level of the 
system. That this might be a consequential difference is not yet clear. 


4.5 Subsidiary Amplifier Gain 


In the development to this point we have given the relationship 
between subsidiary amplifier gain and main amplifier gain in various 
implicit forms. Here we explicitly state this relationship in terms of 
the insertion power gain G,,, of the main amplifier and G, of the sub- 
sidiary amplifier. 

Let us form a skeletal schematic of the levels of power flow in the 
full amplifier, which is shown in Fig. 14. We define the following 
quantities: 


L,,—Loss through input power splitter to main amplifier 
L,—Loss through input power splitter to subsidiary amplifier 
L‘—Loss through the error coupler to subsidiary amplifier. 


If we particularize the circuit between the subsidiary amplifier 
output and the load, 1 + 1/N?, we find the equivalent shown in Fig. 
15, whereas the output of G,, to the load is shown in Fig. 16. The 
amplifiers G,, and G, are each normalized to a unit impedance and, 
through Figs. 12 and 18, we define the transfer functions 7’, and T’ to 
the load, 1 + 1/N*, which we determine to be 








N 
‘i 1+? (33) 
and 
N? 
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Tig. 14— Feed forward amplifier schematic. 
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Fig. 15 — Equivalent circuit as seen by the subsidiary amplifier. 


The nature of feed-forward correction is such that the overall 
transfer power from the input through the main amplifier path 
equals that through the subsidiary amplifier path. This provides 


y A has 
eae ni 
which, through (33) and (34), yields 





G, — N’L,L 
Gr Lm 
Equation (35) determines the differential gain of the subsidiary am- 
plifier over the main amplifier. 
We may quantitatively evaluate the differential gain for a par- 
ticular, but typical, low noise figure amplifier. For this amplifier we 
have 








(35) 


L, = Li = 1.07 — 0.3 dB, 
L,, = 20 — 13 dB, 

and 
N* = 9 — 9.55 dB. 


The differential gain, in dB, of the subsidiary amplifier is, in this case, 


9.55 + 0.3 + 0.38 — 18 = —2.85 dB. 


In this calculation we find, as a practical matter, that the subsidiary 
amplifier actually has a smaller gain than does the main amplifier 
in a low noise configuration. 


ais 
‘ N2 
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Fig. 16 — Equivalent circuit as seen by the main amplifier. 
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4.6 Sampling Via Directional Coupling 

The causal arguments, which motivate the development of this 
form of amplifier, imply that all sampling be done in a strictly or- 
dered time sense. This requires sampling devices which are well 
isolated from reflections within the system which appear as events in 
improper time sequence. 

Figures 12 and 13 shows the use of directional couplers for sam- 
pling. Let us assume for the moment that the output of the main 
amplifier, A, is perfectly matched. Under the circumstances error 
injection by the subsidiary amplifier G into the final transformer has 
only backward component of wave travel through the coupler Ce. 
and no feedback occurs to that amplifier. 

If the output port of A is mismatched, any reflections of the in- 
jected error signal create a new forward wave representing undesired 
signal. However, the undesired signal is sampled by Ce, treated as 
an error, and corrected via the correction loop. 

Hine reflected back into the system from a mismatched ee is 
similarly handled. Thus, in the steady state, the directional character 
of Co, has acted to preserve the prescribed net emergent signal. To 
the extent that a coupler has finite directivity this ideal is blurred, but 
couplers with high directivities are available in practice and the 
idealization may be met to high degree. 

Aside from its use to eliminate parasitic feedback paths, the sam- 
pling coupler stabilizes the main amplifier output impedance and 
provides for a high output return loss. As we have already discussed, 
the fixed amplifier system has an unusually high degree of linearity 
and, as such, is entirely describable by a Thevenin equivalent. Let us 
then, in gedanken fashion, measure its output impedance, which is a 
physically meaningful quantity. 

To measure the system output impedance, a wave emanating from 
an external generator is made incident upon the output terminals, 
and its reflection measured. From Figs. 12 and 18 we can see there 
is a first reflection from the series impedance imaged through the 
1:N transformer, and a second reflection from the main amplifier 
itself. However, as we have seen, this second reflection is negated 
via the coupling through Cz back to the subsidiary amplifier. 

By design, the effect of the transformer is small and any mismatch 
of the subsidiary amplifier is a thoroughly negligible second order 
effect. The net input to the output terminals therefore appears as a 
resistor of normalized value 1/N? in series with a perfect absorber; 


686 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


a total resistive impedance of 1 + 1/N?. Since N? = 7, the mismatch 
to unity is relatively small and a simple matching section may be 
used to produce a large return loss across the band of interest. Most 
importantly, this return loss relates only to the quality of passive 
circuits and is essentially independent of the active system. 

A simple catechism helps us form an effective perspective in rec- 
ognizing the role of directional coupler sampling in error correction 
in relation to the two other, better known, types. 


(t) Voltage sensing produces an unvariable emergent voltage, mak- 
ing the source impedanceless. 

(1) Current sensing produces an unvariable emergent current, mak- 
ing the source admittanceless. 

(wit) Incident wave sensing produces an unvariable emergent wave, 
making the source reflectionless. 


Within this context, the use of directional couplers to provide in- 
cident wave sensing is not limited to only the feed-forward system 
described here to produce a highly matched generator. 


4.7 Parallel Literature on Feed-Forward Technique 


It seems clear that the presently described structure is closely 
related to the previously mentioned class of redundancy and desensiti- 
zation networks.' 7 Nevertheless, we feel that we emphasize certain 
important points in this paper that are contained, but not obvious in 
the general treatments: 


(1) The significant role of matched delay is cited, from causal con- 
siderations, as the basic means of providing a degree of correction 
realization unavailable in feedback. The degree of correction is 
arbitrary. 

(1) The absence of stability considerations in feed-forward. 

(1) Realization of full gain-bandwidth, fully consumed within the 
band of operation, without restriction on the degree of correction. 

(iv) The advantages of the use of unequal amplifiers with respect 
to noise figure. 

(v) Use of coupler sampling both for stability and high output re- 
turn loss. 

(vt) Transformer injection of the error correction to minimize loss 
of power capacity. 
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V. SPECIAL COMPONENTS CONSIDERATIONS 


It is intended to review the coupler art elsewhere; we are not con- 
cerned with it here. Indeed, the couplers that we used were com- 
mercially available items* and were not, themselves, the source of any 
general concerns. 

There were, however, certain specific concerns in the experimental 
developments which did not materialize. While they did not manifest 
themselves for our particular equipment, they still might show up in 
other circumstances. 


5.1 Error Injection Transformer 


While not in exact correspondence with the specific values chosen 
for the mathematical model, a 3:1 transformer is close enough to the 
value needed and generally can be constructed over the necessary 
bandwidth range at the frequencies of interest. The major concerns 
were the practical one: can a ferrite loaded transformer sustain the 
power levels flowing through it?; and the hypothetical one: what if 
the frequency range were beyond the capabilities of ferrites? 

In seeking other realizations of the transformer, one must ap- 
preciate that the transformer shown in Fig. 14 is little more than a 
schematic format of the circuit. Specifically, the transformer is a re- 
active reflection isolating the subsidiary and main amplifiers. Any 
other reactive two-port having the same insertion loss would, accord- 
ing to Weissfloch,® be equivalent to the transformer to within the in- 
sertion of phase corrections at each of the ports. Further, the error in- 
jection might take the form of a dual shunt current feed as opposed to 
the series voltage form and, as such, there are quite flexible options in 
design. 

While there is hardly a dearth of networks to produce the desired 
circuit function, we shall develop one which seemed simply suited for 
synthesis. Figure 17(a) shows a tee of inductances. The power trans- 
mission |T?| vanishes at zero and at infinite frequency without any in- 
tervening singularity so that, as shown in Fig. 17(b), a stationary 
transmission value exists. If the circuit is modified as shown in Fig. 
18(a), by introducing the shunt capacity C”, a high frequency propa- 
gating region is introduced as observed in Fig. 18(b). With appro- 
priate proportioning of the capacitance an inflection region is intro- 
duced at the point of stationarity and a flat insertion loss region is 
obtained. This flat insertion loss may be corresponded to the trans- 


* Merrimac R&D Co., West Caldwell, N.J. 
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Fig. 17 — Network with stationary transmission. 


former N within the Weissfloch equivalent and the various component 
values are shown in Fig. 19. 

For completeness we provide the relevant formulations. If L is the 
numeric power insertion loss, we have 








n=1+3(w-4); (36) 
= 4 N/’ 
a : V/3L” ( w vaL') ae |: 
L=1+ E = 1+ 6L’? — 3L' /3L’ ‘s Lol) |; 
(37) 
pi oe ie 3L”" : 
L = 6L ’ (38) 
and 
CN" = an (39) 


t+ on — an 

The net results of these machinations are shown in Fig. 20. Since 
the network has a relatively high inductive impedance at midband, a 
short transmission line 6 is used to present a resistive input. Thus input 
impedance is a minor perturbation and frequency effects are corre- 
spondingly small but, if desired, a low Q resonator matching network 
in tandem with the load may be employed for higher degree matching. 





Fig. 18 — Modification to produce an inflection point. 
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Fig. 19 Component values for equivalent transformer. 


The load conductance G; is equal to 1 + 1/N? in dual fashion to series 
injection form shown in Fig. 21. 


5.2 180° Phase Differential Networks 


The 180° phase differencing networks are staples of the interconnec- 
tion formats we have described here. They can be realized to a good 
degree in the tens of magacycles range as a pair of mutually reversed 
ferrite transformers together with capacitive compensations to account 
for the differences in the way in which the winding capacitance is en- 
ergized.*° There was fear that the power handling capacity might prove 
excessive for the ferrite or that the parasitics might dominate. These 
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Tig. 20—Shunt error injection using transformer network equivalent. 
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Fig. 21— Thevenin equivalents of subsidiary amplifier feed into main amplifier 
output circuit. 


were exactly the anticipatory fears with relation to the injection trans- 
former and, in precisely like fashion, the fears proved groundless at the 
power level and the frequency range of operation. 

Nevertheless, there is the question of what to do in the event that 
ferrite transformers are unrealistic. We use a single section quadrature 
coupler together with a technique credited to E. M. T. Jones to 
produce the desired realization. Jones’s technique recognizes that 
the transmission between the input and its conjugate port, in a sym- 
metrically mismatch terminated quadrature coupler, is directly propor- 
tional to the reflections terminating the coupled ports. The transmis- 
sion in such a device whose coupled ports are short-circuited differs by 
180° of phase from that same device coupled into open circuits. 

Equations (10) and (11) show a single section coupler used in such 
application to have only a relatively narrow band for high return loss 
and that the transmission is more akin to that of a resonant section. 
Further analysis shows, indeed, that the identification with a resonator 
is a very good one and that the conventional technique of adding prior 
and subsequent resonators to form a filter chain is entirely applicable. 

With the coupled ports open-circuited, the equivalent circuit is ex- 
actly a series resonator with both the inductance and capacitance cor- 
responding perfectly to unit impedances at band center. A maximally 
flat three section filter chain is constructed by bounding this resonator 
on both sides with shunt resonators composed of elements possessing 
exactly double unit impedance. 

If the coupled ports are short circuited, a dual situation results. The 
coupler equivalent is a shunt resonator composed of unit admittances 
and the bounding networks are series resonators composed of elements 
of twice unit admittance. 
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Since the open and short circuited couplers are perfectly dual struc- 
tures, the 180° phase difference is perfectly maintained irrespective of 
frequency. The bandwidth of the structure pair is about 3:1 before the 
onset of any significant transmission loss; this results in an eminently 
practical component. 


VI. REALIZATIONS 


6.1 Motivation 


The above philosophies and design considerations have been realized 
in the construction of two amplifiers intended for testing dispersive 
ultrasonic delay lines in the 60-90 Me and 25-35 Mc ranges. The large 
loss inherent in these devices made necessary large power output while 
the dispersive nature of the delay imposed the requirements of linear 
flat gain with little mtermodulation distortion. Since the delay line 
represents a nearly lossless capacitive load on the amplifier, the latter 
must be a well matched device in order to prevent large interaction ef- 
fects between amplifier and delay line. 


6.2 60-90 Mc System Performance 


The first development undertaken was the 60-90 Mc amplifier. 
During the course of this project we learned much about the operation 
of a feed-forward system. This increased sophistication was applied to 
the subsequent 25-35 Mc system. The first amplifier is shown in Fig. 22, 
and it is seen to be a single stage feed-forward system with its input 
sampling loop and its output error correction loop which is capable 
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Fig. 22 —75Mc feed forward system. 


692 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


of supplying one quarter of the main amplifier power. This amplifier 
delivers 10 watts average power in CW operation (20 watts peak) to 
the 0.1 dB compression point. 

At this power level the phase distortion is 1.25 degrees (see Fig. 23) 
at 75 Mc. The swept response across the band is shown in Tig. 24 both 
for small signals and at maximum power output. (The dashed curve 
shows the system response with the correction circuit disabled.) A 
sensitive measurement was made of the small signal gain and phase 
across the band (Figs. 25 and 26) which showed +-0.1 dB gain ripples 
superimposed on a 0.5 dB slope* and +2° phase deviation from con- 
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Fig. 23 —System gain and phase at 75Mc as a function of input power. 


stant time delay (linear phase vs. frequency) out of a 1300° phase 
variation across the band. 


6.3 Couplers in 60-90 Mc System 


The couplers, used to pair amplifiers in the preamplifiers and in the 
more general harnessing in the main power amplifiers, were double 
stage quadrature hybrid couplers purchased from Merrimac Research 
and Development Co. (see Sections 1.2 and 1.3, and Fig. 27. The trans- 
fer characteristics of this coupler are shown in Fig. 28 and the input 


* The +0.1 dB ripples may have been partially due to a procedural error in 
the operation of the test set used for the measurement. 
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match with two identical total reflections in Fig. 29. This last. figure 
points up one of the main advantages of the quadrature coupler type 
of organization; the fact that input match is maintained for arbitrary 
port termination as long as the two terminations are identical. 
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Fig. 24 — 75Mc system output as a function of frequency. 
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Fig. 25 — Gain variation with frequency at low power output. 


6.4 Amplifiers in 60-90 Mc System 


The preamplifiers shown in Fig. 22 are composed of an input pair 
of low-level high-gain amplifiers and enough stages of intermediate 
power emitter follower pairs emitter follower pairs (Fig. 30) to pro- 
vide sufficient power to drive the main power array of emitter fol- 
lowers (Fig. 31).* All emitter followers are biased to class A operation. 


* The amplifiers used were based on designs and models provided by R. V. 
Goordman and R. C. Petersen. 
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Fig. 26— Phase deviation from linear phase vs. frequency behavior (Constant 
time delay). Average time delay = 2.5usec. 


The high gain amplifiers are of a broadband design with a flat 52 dB 
gain from 0.2 Mc to 100 Me and capable of 100 milliwatts output be- 
fore compression at the higher frequency (greater power at lower 
frequency). 

The drive stage emitter followers in the preamplifiers used 2N3375 
transistors they have sufficiently high base impedance at these fre- 
quencies to permit some modest input transformation in order to 
enhance the 6 dB theoretical gain made available by driving the high 





Fig. 27 — Quadrature hybrid coupler, 
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Fig. 28 — Coupling characteristics of 2 section coupler. 


impedance emitter follower input from a matched source. The driver 
stage amplifiers had gains of 8.5 to 9 dB and were capable of 0.6 to 
0.7 watts output. The 2N3632 transistors used in the power stage 
amplifiers, however, did not allow any transformation and, after 
some input tuning, had 6-7 dB gain but were capable of 1.25 watts 
output. 

In order to realize the matching properties of the coupler organi- 
zation it was found necessary to select matched pairs of transistors. 
The transistors were inserted in a test jig and tuned for best reponse 
with calibrated reactances. The values and final response were noted. 
Pairs of transistors were selected by comparison of these data. The 
amplifier pairs constructed using them showed typically > 25 dB 
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Fig. 29 — Return loss of 2 section coupler with open circuit terminations on 
coupled ports. 
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Fig. 30-— 60-90 Me class A, 8.5-9.0 dB gain, 600-700 mW linear driver stage. 


input and output return loss across the band. Interaction between the 
amplifiers outside of the band is prevented by using “out-of-band” 
filters as described in Section 2.6 (see Fig. 22). 


6.5 Szgnal Processing Components in 60-90 Mc System 


In the feed-forward system, as developed in the preceding text, 
there are six essential signal processing operations (see Fig. 22): 


(1) Splitting off a reference signal from the input 

(wi) Delaying this reference signal for a time equal to the time of 
passage through the main amplifier 

(wz) Splitting off a sample of the amplified signal 

(tv) Comparison of the sample with the reference 
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Fig. 31 — 60-90 Mc class A, 6.0-7.0 dB gain, 1.25 watt linear power stage. 
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(v) Delaying the amplified signal for a time equal to the time of 
passage of the error through the subsidiary amplifier 
(vi) Injecting the correction into the output path. 


In the 60-90 Me system operations 7 and w are carried out within 
in-phase power dividers (Merrimac Research and Development Co. 
PD-75). Operations 7% and v are carried out in approximately 60- 
and 30-foot lengths of coaxial cable. (These lengths were first calcu- 
lated from measurements of phase vs frequency within the system.) 
Operation 722 is carried out in a two-section quadrature hybrid coupler 
similar to the ones used within the amplifiers (see Figs. 27, 28 and 29). 
However, instead of a 3 dB power split, this coupler was synthesized 
to couple out a signal 15 dB down from the input and to present only 
0.27 dB loss to the main power path (see Fig. 32). The final opera- 
tion, (vz), uses the reactive network described in Section 4.1 and 
Fig. 18. 

The other components appearing in Fig. 22 are of secondary im- 
portance to the feed-forward concept but are quite essential to sys- 
tem operation. The fixed pad in the main amplifier input and the 
variable attenuator in the sample path insure that the power levels 
are appropriate to a proper differencing in the power divider. The 
variable attenuator in the error correction loop performs a similar 
function. These variable attenuators were constructed based on the 
transmission and match properties of the quadrature hybrid coupler. 
Figure 33 shows the transmission between the input and hybrid ports 
of a two-stage quadrature coupler both under open circuit termina- 
tion of the coupled ports and with 10 © terminations. The use of 
ganged carbon rheostats as terminations (Fig. 34) produces a well- 


POWER 





Fig. 32 — Coupling characteristics of hybrid directional coupler used as sam- 
pling coupler. 
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Fig. 33 — Transmission of two section coupler under open circuit and 10 Q 
termination. 


matched variable attenuator which is quite flat over the modest at- 
tenuation ranges needed. 

The “sample coupler equalizer” is a low Q shunt resonant circuit 
which has a transmission characteristic that compensates for the bow 
in the sampled signal from the coupler, as shown in Fig. 35. The 
“phase shift networks” shown in Fig. 22 are also resonant equalizers. 
When the cable length corresponding to 7; was adjusted such that 
the time delay through the main amplifier equaled that through the 
reference path (equal phase vs frequency slopes), it was found that 





Fig. 34— Variable attenuator. 
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Tig. 35—Sampling coupler equalization. (a) Sample coupler characteristics. 
(b) Transmission of equalizer. (c) Equalized characteristics. 


there was a constant 70° excess phase shift through the amplifier out 
of a total variation of 1500° across the band (see Fig. 36). In order 
to equalize this difference, use was made of the fact that a resonant 
circiut has zero phase shift at the resonant frequency, at which point 
there is an inflection in the phase vs frequency characteristic (see 
Fig. 37). The latter produces constant time delay at the resonant fre- 
quency. An extra length of line was then added to the reference path 
to produce equal phase shift at center frequency and the two-reso- 
nator network shown in Fig. 38 added to the sample path to equalize 
the phase slope (see Fig. 36). It was found that a two-resonator net- 
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Fig. 836 — Phase equalization of reference loop. 
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Fig. 37 — Amplitude and phase of transmission through a resonator. 


work was needed, minimally, to produce sufficient flatness across the 
band. The insertion of the phase shift network, ye, in the correction 
loop relative to rz is similar to the above. 


6.6 System Checks 


With the system assembled, each of the two loops was put through 
an operational check, and final adjustments were made. The ref- 
erence loop is checked by applying a small signal to the input and 
observing the output of the comparison power divider at B. The 
variable attenuator in the power sampling line was adjusted to pro- 
duce the best null across the band. Small adjustments of 7; were 
made and the best null was obtained by using a sweep generator 
and an oscilloscope. 

The error correction loop was then adjusted by feeding a swept 
frequency signal into the main amplifier at A and observing the sys- 
tem output. Both the variable attenuator at the input to the sub- 
sidiary amplifier, and r2, were adjusted to produce a null output. 


0.065pH 7OpF 
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Fig. 38 — Phase shift section (values shown for 75 Me, 70° correction) in 50 
Q system. 
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The final experimental proof of system performance lies in breaking 
(and terminating) the reference line at C in order to “fool” the error 
correction loop into operating as if there were no input signal, and 
observing the entire system between input and output ports. The ex- 
cellence of the correction system may be judged by the degree of 
degeneration of the output. The results of this observation with the 
subsidiary amplifier turned off, and with the subsidiary amplifier re- 
stored are shown in Fig. 39 along with a simplified schematic indicat- 
ing the measurement condition. 


6.7 25-85 Mc Amplifier—Differences from 75 Mc System 


After experience with the 75 Me system, the 30 Mc amplifier was 
designed to increase the power output several fold over the 10 watts 
attained above, without increasing the complexity of the system 
(that is, the number of coordinated branches in the power amplifier). 
This was done by going to class B push-pull operation rather than to 
a single-ended class A operation as in the previous system. In addi- 
tion to doubling the number of transistors, the de-rf power conversion 
capability of each was doubled. The push-pull pair (see Fig. 40) 
delivered 7.5 watts, linearly. 




















Fig. 39— System output with reference line open, feed forward switched on 
and off. 
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Fig. 40 — Class B push-pull amplifier for 30 Me system (20-40 Me, class A-B, 
push-pull, 7.5 watts linear). 
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This transition to class B operation made another design change 
possible. It was decided to design for low duty cycle pulsed operation 
and to use an “active clamp” biasing circuit which samples and re- 
stores bias voltages between pulses. 

In addition to extending the dynamic range of the amplifier by 
increasing its power output capability, an attack was made on the 
system noise figure by the use of a directional coupler at the input 
favoring the reference path. As a further improvement the series in- 
put matching resistors (~43 ©) of the low-level preamplifier in the 
subsidiary amplifier were removed. The input match was restored by 
the amplifier pairing operation using quadrature couplers, which 
produced a matched amplifier having approximately a 4 dB noise 
figure. 

Other changes were in the use of a series 3:1 injection transformer 
instead of the previously used shunt injection network, and in the 
use of limiting in the low-level preamplifier sections of both amplifier 
chains to prevent destruction of the power transistors under accidental 
overdrive. A schematic diagram of this system appears in Fig. 41. 


6.8 ‘Active Clamp” Bias Circuit 


The active clamp biasing used in the 25-35 Mc amplifier (see Figs. 
42 and 48 provides a means of stabilizing the amplifier against changes 
in bias condition with variations in the signal power level and duty 
cycle. At high signal levels an unrealistically large bypass capacitor 
(see Fig. 44) would otherwise be required to perform this function. 
The clamp circuit also serves to make the bias levels insensitive to 
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Fig. 41— 30 Mc system. 
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Vig. 42 — Active clamp schematic. 
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Fig. 43 — Active clamp biasing of power transistor. A, B, and C are connections 
to the other transistor of the push-pull pair. 
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Fig. 44 — Conventional transistor biasing arrangement with DC stabilization. 


thermal drift in the controlled transistor and to parameter variations 
between substituted transistors. 

The high gain differential amplifier, shown in Fig. 42 provides the 
self-correcting bias condition since its output is fed into the base of 
the power transistor in such a way as to keep the potential at A of 
Fig. 3 equal to the reference potential at B. The gain of the differential 
amplifier also allows the use of a much smaller value of R,; (15 2 as 
in Fig. 48) than would be used in a conventional current feedback bias- 
ing scheme (shown in Fig. 44) to sense variations in quiescent current 
through the power transistor. 

Ordinarily, as stated, a large bias bypass capacitor C, shown in Fig. 
44 and dotted in Fig. 42, must be used to provide base current to 
the transistor during the conducting part of the cycle. It must be 
made adequately large to prevent substantial discharge, independently 
of duty cycle or drive amplitude. Transistor Q4 of Fig. 44, allows the 
use of a much smaller capacitor, C,, since this capacitor need only 
supply a control current to its base. The transistor, Q3, prevents the 
discharge of C, into the differential amplifier during the pulse and 
allows a rapid recharge of C, after the pulse because of its low emit- 
ter follower output impedance. 


6.9 Performance of the 25-35 Mc Amplifier 


The characteristics of the completed amplifier at low power are 
shown in Figs. 45 through 48. These figures show that the gain of the 
amplifier slowly oscillates +0.1 dB about a 0.4 dB slope across the 
band, the total amplification being 39 dB. The phase varies +1° from 
a constant time delay phase characteristic. Further, the spot noise 
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Fig. 45 — Deviation from flat gain at low power. 


figure is better than 12.5 dB across the band and better than 9.5 dB 
over 85 percent of the band. The integrated noise figure is 8.6 dB. 

Figures 49 through 51 show the gain and phase as functions of 
power level with the system fully functional and with the feed- 
forward correction loop disabled. These figures show the marked 
improvement obtained by the use of feed-forward correction. Notice 
that, as implied by the analysis in Section 3.8, when the subsidiary 
amplifier is called on to deliver large power to compensate for gain 
variation, it can no longer provide phase compensation. (It was 
shown that the compensation of only 3.6° at maximum power output 
requires all the power capability of the subsidiary amplifier.) The 
40 watts indicated on the figures and the 8.6 dB noise figure shown 
above yield a dynamic range for the amplifier of 103.5 dB. Its output 
return loss was > 25 dB (see Fig. 52). 


1200 





PHASE SHIFT IN DEGREES 
o 
° 
°o 





QO 
25 26 27 28 29 30 31 32 33 34 35 


Fig. 46 — Phase shift through 30 Mc amplifier system as function of frequency 
(low power). 
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Fig. 47 — Deviation from constant time delay (low power). 


VII. CONCLUSIONS 


The amplifier systems that were constructed at 30 and 75 Me had 
an excellent correlation with design and showed the engineering prac- 
ticability of the methods of design construction. The primary features 
were the use of: 


(7) Quadrature coupled corporate amplifier arrays to phase co- 
ordinate power emission of the individual active elements 

(it) Emitter follower amplifiers as the basic elements to provide 
broadband, low distortion gain with transistors, otherwise incapable 
of producing this type of performance with a high yield 

(iii) A time-compatible or “feed-forward” error injection tech- 
nique to equalize the transfer function of the amplifier, reduce distor- 
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Fig. 48 — Amplifier noise figure, 
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OUTPUT POWER IN DECIBELS 


Fig. 49 — Gain (a) and Phase (b) Variation at 28 Mc. 


tion and diminish noise in a system having relatively large time delay 
(iv) Directionally coupled error sensing to produce a high return 
loss at the amplifier output 


The 75 Mc amplifier was operated over a band of 60-90 Me pro- 
ducing a 10 watt output average power with an amplitude compres- 
sion there of about 0.2 dB. Its transfer function showed slow amplitude 
variations within 0.1 dB and phase deviations from linearity of about 
1.5° out of a total phase spread of greater than 1500° across the band. 

The 30 Me amplifier had a band of operation between 25 and 35 
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Me. It produced a 40 watt average power with an amplitude compres- 
sion of 0.2 dB. Transfer characteristics showed slow amplitude excur- 
sion of about 0.2 dB with phase distortion about 1°. Output return 
loss was about 25 dB across the band and lent itself to simple 
equalization. 

The 30 Me amplifier had an 8.6 dB noise figure across the 10 Me 
band of design interest which, together with its 40 watt average un- 
distorted output capability, provided a dynamic range of 103.5 dB. 

While correction did not extend beyond that of a single loop in 
these amplifiers, the textbook-like adherence of the results to the 
design elements makes it clear that further feed-forward looping 
would provide an arbitrary degree of correction over the prescribed 
bandwidth, probably with some further noise improvement. This is 
to be contrasted to the inabiilty to use a feedback control for these 


f 40.8 WATTS 
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Fig. 50 — Gain (a) and Phase (b) Variation at 32 Me. 
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Tig. 51 — Gain (a) and Phase (b) Variation at 85 Me. 


amplifiers specifically because of excessive delay over the bandwidth, 
and to the limited ability of feedback, in general, to degenerate 
error cumulatively because of stability concerns. 

The feed-forward correction scheme does not impose any concern 
for active element characteristics beyond the band of use, as does 


RETURN LOSS IN DECIBELS 





20 22 24 26 628 30 32 34 36 38 
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Fig. 52 — Output return loss of 30 Mc amplifier. 
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feedback, providing the simultaneous ability to use higher power and 
cheaper transistors than would ordinarily be acceptable at the fre- 
quency used. Using the self-degenerating emitter-follower mode of 
amplification, experience with the two amplifiers constructed showed 
a high degree of insensitivity to transistor selection aside from those 
failing catastrophically. Nothing more than a de transistor checker 
was needed to construct the 30 Me unit, while a rough rf sweep selec- 
tion procedure was employed at 75 Mc to form slightly preferred 
transistor pairs. 

The use of time-compatible error cancellation is not limited to the 
specific circuits described here either in frequency or in the specific 
active element circuitry used. It is generally applicable and, in con- 
junction with the corporate power accumulation available through 
coupler arrays, and smooth (no guard-band) frequency multiplex ar- 
rangements, there are available at least in principle, controlled, ar- 
bitrary (bandwidth, arbitrary power, capabilities virtually through- 
out the entire rf spectrum. 
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APPENDIX A 
Conditions for Directional Coupling 


In Section 2.2 we indicated that a sufficient condition that a pas- 
sive, reciprocal, and reactive, fourport possess directional coupling 
properties was that it be matched and that it possess two axes of sym- 
metry. We shall show the sufficient conditions to result from the ap- 
plication of simple time reversal arguments. 

Figure 53 shows the fourport excited by a unit incidence into port 
1 without reflection but with scatterings a, 8, and y, into ports 2, 3, 
and 4, respectively. If this figure is time-reversed, being reactive, the 
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Fig. 53 — Prototype network. 


fourport remains unchanged, but the scattered quantities are reversed 
in flow direction and their phases inverted. This result is shown in 
Fig. 54 which then must be an excitation format compatible and 
simultaneous with that of the first figure. 

Use in turn incident amplitudes, —8*, —y*, and «*, into ports 3, 4; 
and 2, respectively, to cancel the various incident quantities of Fig. 
54. The respective scattered results are shown in Fig. 55, where the 
two axes of symmetry and the scattering format of Fig. 53 were used. 
If the results of Fig. 55 are superposed, then the excitation scheme 
of Fig. 56 is yielded. This last figure shows only emergent waves. Since 
there are no incident waves, each of the emergent waves must vanish 
and we have 


1 — aa* — pp* — yy* = 0, (40) 
ay* + ya* = 2 Re ay* = 0, (41) 
aba + Ba* = 2 Re Ba* = 0, (42) 


and 


I 
© 


yB* + By* = 2 Re y6* (43) 


Equation (40) simply confirms the conservation of energy in a 
reactive fourport. Equations (41) through (48) state that a, 8, and 
y are all mutually perpendicular phasors. Since these three phasors 





Fig. 54— Time reversed description. 
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Fig. 55 — Excitation of ‘8’ Port (a), “y’ Port (b) and ‘e’ Port (c). 
are in a two-dimensional space this mutual perpendicularity is absurd 
unless one of the quantities vanishes. Therefore, 


(1) There exists one port “hybrid” to the incident port 
(iz) All the power is transmitted to the remaining two ports. 


These two statements define a directional coupler. 


APPENDIX B 


Transformer Injection 


B.1 First Considerations 


B.1.1. In Phase Errors 


Let us first see what is required of the circuitry of Fig. 14 to ac- 
complish in-phase error correction then consider the correction of 





Fig. 56—Superposed excitations. 
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quadrature phase errors. Our specific concern is the specification of 
the 1:N transformer shown as the means for final injection of the 
error cancellation signal. 

Figure 21 shows a first Thevenin equivalent of the insertion of the 
subsidiary, error correction, amplifier into the output circuit of the 
main amplifier. Both amplifiers are shown with normalized internal 
impedances and the external load matched to the combined imped- 
ance. We assume proportional error 8, of the main amplifier open 
circuit voltage 2V, and we assume further that the open circuit voltage 
of the subsidiary amplifier is 2v. 

As seen from Fig. 21(b) the reduced equivalent, the condition that 
the load be excited independently of 8 is given by 


v = NVs. (44) 


If we isolate the subsidiary amplifier, as in Fig. 57, the amplifier has 
a voltage v’ at its load terminals corresponding to a load current 
V/N(1 + 1/N?). At the maximum voltage of V, defined as Vp, with 
corresponding subscript notation on all other quantities, we would like 
to see the subsidiary amplifier match terminated to deliver maximum 
power. Accordingly, v’ is half the open circuit voltage and it is equal 
to the current at match in the impedance normalized description. Then 





Up = scale 5 (45) 
va + 3) 
and, using (44), we have 
N = fice (46) 
and 
vp = V5 V 5,(1 _ 6,). (47) 
1 Caterers 
N(1+ xe) 
2vVv v' 


Fig. 57 — Thevenin equivalent of loaded subsidiary amplifier. 
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By the above construction, at peak power the subsidiary amplifier 
works into a match and furnishes a power P, which compares to the 
main amplifier P,, , via (46), by the relationship 


bp 
Go: 





Pia SCO a P, (48) 
where V2(1—6,)”=P,,, is the peak available power of the main amplifier. 

We may consolidate our results to this point under the approxima- 
tion that 8, is small. 


(1) The main amplifier power deviation is approximately [ (1—8,)? 
— 1]P pn =~ — 28,P,, at its peak power capability. 

(it) The subsidiary amplifier furnishes a correction power of ap- 
proximately 8,Pn. 


There seems a curious divergence of a factor of 2 in providing the 
necessary compensation. If, for example, 8, = 4, implying a peak- 
power deviation of 4 P,,, a subsidiary amplifier having only a power 
capacity of 4 P,, would appear to be adequate for compensation. 

A second feature of transformer injection may also be discerned 
which at first would appear unrelated to the missing factor of 2. The 
error injection is not without some cost to the main amplifier. As Fig. 
21 shows, the subsidiary amplifier images a resistance of 1/N? into the 
main amplifier circuit, causing an effective generator characteristic im- 
pedance of 1 + 1/N? to be presented to the load. 

The subsidiary amplifier provides the power both to equalize the 
irregular main-amplifier transfer characteristic and to cancel distor- 
tions produced by nonlinearities. Assume the main-amplifier power to 
be small such that no nonlinear distortions are observed and assume 
further that equalization is unnecessary. There is, therefore, no er- 
ror and, consequently, no contribution from the subsidiary amplifier, 
and the open circuit voltage is that of the main amplifier only. The 
available power is decreased by an amount equal to the changed 
impedance and an insertion loss of 


4 Atl 

“1 6, 

is sustained where the results of (46) have been used. 
The question exists whether this is indeed just a loss in gain or a 

loss in power capability. We shall now show that it is only a loss in 

gain and that we add to the power capability of the main amplifier 

by the amount of the subsidiary amplifier peak power, 


1 
D=1+ N? (49) 
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From Fig. 21 we see the delivered voltage of the main amplifier to 
be, at peak, 





Vir 
1 
1+ ap 


2V,(1 — 6,) — 


Since 1 + 1/N? = 1/1—8, by (46), the voltage delivered reduces to 
V, (1—8,). Since the current drawn is 2V,/2(1 + 1/N?) = V, (1-8), 
the total power delivered by the main amplifier is 


Pa= Vid 6: (50) 


However, (50) is, identically, the available power from the main am- 
plifier without the subsidiary circuit, so that there is no difference in 
the power actually delivered by the main amplifier into the load. 

The power delivered by the subsidiary amplifier is the product of 
its output voltage referred to the main loop, 2v,/N — [(1/N? x V,/ 
(1 + 1/N*)] and the main loop current V,/(1 + 1/N?). The use of 
equations (46) and (47) shows this power to conform to (50) produc- 
ing a power 


P, = Vi(1 — 4,)6,. 
The total power reaching the load is P,, + P,; = Pz which is 


Ps 
Lae 


Equation (51) shows most succinctly that the power reaching the 
load exceeds the power capability of the main amplifier and that the 
transformer interconnection, as characterized earlier, is vastly super- 
ior to a power divider interconnection which would have produced a 
3 dB loss in power capacity. Further, the insertion loss shown in (48) 
is not a power loss, but a small loss in gain only. 


P, = V3(1 — 6,) = 





(51) 


B.1.2. Quadrature Errors 


Let us assume that the transformer N has been specified by (47) 
and the subsidiary amplifier by (48) in terms of some desired real 
peak voltage error 6,. We now seek to determine the peak capability 
of the system to correct a quadrature error 76/ . 

Given this quadrature error, the subsidiary amplifier output voltage 
is v = 1NV,6, corresponding to its current 7 = V,/N. The subsidiary 
amplifier perceives, therefore, a terminating reflection factor 
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sock ie iV 
pg ay 
pear iNV, 8, +4 


which, obviously, has a unit magnitude. 

Given a unit reflection, the subsidiary amplifier, depending on its 
electrical position relative to the load, might be called on to develop 
twice its rated voltage or twice its rated current. To safeguard against 
this possibility, the open circuit voltage must be down rated by a fac- 
tor of 2, and corresponding to (47), we have 


of = Te V5 — 8) = NV, 8. (52) 


So that, by (46) 
5, = 36, - (53) 


Quantitatively, if 6, = 4, corresponding roughly to a +1 dB cor- 
rection, (53) yields 64 = 0.0625. Since the phase correction has a mag- 
nitude of the arctangent of 6/ , the phase cleanup range is, minimally, 
3.6° at high power. 


B.2 Further Considerations 


B.2.1. Subsidiary Amplifier Mismatch 


At peak power, by design, the subsidiary amplifier is matched. Let 
us now go to the extreme of low power. Using Fig. 21(b) the main 
amplifier supplies a current V/(1 + 1/N?) which develops a voltage 








- 1 (Z)o -— 1 (54) 
1 + 72 N+ 


across the step-up side of the injection transforms via the subsidiary 
amplifier image impedance 1/N?. The main current reduces by 1/N, 
and a current flow also equal to V/(N + 1/N) flows through the volt- 
age step-up side of the transformer. . 

As at high power the subsidiary amplifier output voltage is equal 
to its current, but with a major difference. The voltage sense is re- 
versed and now, instead of finding a positive resistance match, the 
termination is a negative match. In transition between high and low 
power, the termination of the subsidiary amplifier has passed through 
violent shift and a question exists whether, as power decreases, some 
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requirement existed on the subsidiary amplifier worse than that at 
peak main amplifier power. 

Having aroused this concern, we shall now dismiss it. Equation (54) 
shows the negative voltage bound to be less than the positive voltage 
bound, since V, by hypothesis, is less than its peak voltage V,. Simi- 
larly, the subsidiary amplifier current is always less than its peak 
value V,/(N + 1/N). The subsidiary amplifier, in actuality, is sepa- 
rated from the terminating transformer by a transmission line. Along 
this line, peak current and voltage amplitudes interchange. Since 
neither of these peaks ever exceeds the value supplied by the subsidiary 
amplifier at peak main-amplifier distortion, the initial design of the 
subsidiary amplifier covers the worst case. 


B.2.2. Nonoptimal Transformers 


While the transformer design chosen optimizes the powers delivered 
by the main and subsidiary amplifiers, nevertheless one may wish to 
forego this optimization somewhat to reduce the contribution of 1/N? 
to the output matching problem. While this problem is not particu- 
larly serious, it is generally, a nuisance in the following sense. Very 
high return loss is guaranteed by the directional coupler sampling 
arrangement exclusive of the 1/N? contribution, and the very inser- 
tion of any matching circuit at all is a severe perturbation on the 
goodness of the system. The less the mismatch to be corrected, the 
higher the degree of return loss realization across the band. 

Let us now consider a nonoptimal transformer n replacing the op- 
timum value N deduced in earlier expressions. If 2v’ is the open cir- 
cuit error voltage applied to the high side of the injection transformer 
then 


dv, = 2nV,6, 
and the current referred to the high turns side is 
ecole 
1+” 
n 


The terminating transformer, together with the correlated current 
flow in the low turns side, provides an output impedance Z for the 
subsidiary amplifier which must conform to the relationship 


2u), 


2 2S ngs 1 
l+Z2=7 = mi{i +3 on 
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and 

Z = 2(n? + 1)6, — 1. (55) 
I.quation (55) shows that for (n?+1) > 1/8, the subsidiary amplifier, 
as expected, operates into a high impedance mismatch. 

The subsidiary amplifier must be made to cope with the standing 
wave, to be capable of supplying either an excess open circuit voltage 
or an excess short circuit current. We assume a power capacity to be 
required from the subsidiary amplifier such that it can tolerate that 
peak through its internal impedance. Since this peak, relative to 
matched loading, is given by the ratio of 2Z2/(Z + 1), namely 1 + [k|, 
where k is the reflection factor, we have 








ve Ww+tis—in 
in QE Dk (56) 
with the corresponding peak subsidiary amplifier power ratios 
Pe, 2 + De =A) we 
Poe. ( (n® + 1) 6, N? ep) 
Combining (57) with (46) we have, finally, 
, 2 a 2,2 
Pie CG ed) on (58) 


Py @e ely le) - 

The formulation of (58) shows, of course, that when n? = N? = 
1 — 8,/8,, this ratio reduces to unity. As a quantitative determination, 
let us increase n such that it introduces roughly one-half the mis- 
match it inserts under optimal conditions for 8, = 4, our 1 dB design 
center. Since N? = 7, let us choose n? = 16, where the latter number 
is chosen to form an integral turns ratio. Equation (58) demands 7.3 
times as much power as that for the matched case, a level of power 
for the subsidiary amplifier virtually equal to that of the main am- 
plifier. 

The quantitative example chosen, of course, represents a fairly 
undesirable situation; nevertheless, it shows a tradeoff possibility 
which is not to be ignored. In multiple looping feed-forward systems, 
each loop is called on to handle progressively smaller and smaller 
deviations. In using these progressively smaller values of 8,, we form 
a sequence of progressively higher values of N, corresponding in turn, 
according to (48) to lower values of P,. Below some value, the value 
of P, is no longer limiting and we may wish to use correction power 
as a tradeoff against mismatch. 
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APPENDIX C 


Schmidt-Briicken Patent 


During the final stages of preparing this article we found a German 
Federal Republic patent, issued July 14, 1960 (document 1085194), 
to Heinrich Schmidt-Briicken, entitled “A Method For Compensating 
Distortion In Amplifiers.” 

Clearly this patent, in spirit, precedes our own activities in attempt- 
ing to provide time compatible error cancellation. While it suggests a 
time flow correction system seemingly much similar to our own, its 
actual embodiment severely compromises the intent of the patent. We 
have reproduced Figures 1 and 3 of the Schmidt-Briicken patent as 
our own Figures 58 and 59. We claim the disclosure shown in Fig. 58 
to be operative only in principle, and that of Fig. 59 to be generally 
inoperative. 

We feel it relevant to pursue the dissection of these devices since 
they illustrate, negatively, important considerations in the format of 
implementation of feed-forward amplifiers. It is not, however, intended 
by so doing to deprive Mr. Schmidt-Briicken of a well earned recogni- 
tion of his priority in recognizing time as a controllable factor in 
error cancellation. 

A major consideration lacking in both the configurations of Figures 
58 and 59 is that of directional wave sampling. Let us observe the 
consequences of a simple voltage sampling as used in Fig. 58 where 
such sampling follows the usual practice in feedback amplifiers. 

If all components in Fig. 59 are matched, that system will operate. 
In general, however, there is no disclosure of the degree of match 
required in the patent, and the effect of multiple reflections are 
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Fig. 58 — Schmidt-Briicken Fig. 1. 
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Fig. 59 — Schmidt-Briicken Fig. 3. 


ignored. Let us assume, however, that all devices are matched only 
to within reasonable specification of return loss, and that the ampli- 
fiers, particularly, might have return losses of the order of 20 dB. 

The error reference is established at point A in Fig. 58 through 
the use of a voltage divider network, and, consequently, the error volt- 
age is a function of the input impedance seen looking into the delay- 
line labeled 6. The phase of the 10 per cent reflection factor of the 
distortion amplifier 5 may rotate through more than one cycle as 
a function of the delay line length, producing a peak to peak imped- 
ance variation of 50 per cent. The reference sampling is modulated 
by this order of magnitude, making a shambles of the compensation 
system. 

Inherent in this system, therefore, is an inclusion of an active de- 
vice, namely the amplifier, as part of the reference standard, greatly 
limiting the general practicability of this class of device. In our own 
embodiments, we have carefully employed directional coupler sam- 
pling, stressing the incident wave, and not the voltage, as the reference.. 
With properly designed match, third tame around incidence is not 
greatly significant since it is down by twice the return loss, the order 
of 40 dB or greater, making amplifier mismatch a secondary con- 
sideration. 

The same form of argument applies to point B with respect to the 
adder 2. Notice that both boxes 4 and 2 are bridges or directional 
couplers, and that they are used to prevent spurious coupling between 
chains. Awareness is manifested, therefore, within the patent itself 
of the existence of directional coupling devices, but nowhere are these 
couplers applied to the crucial problem of error sampling. 

Figure 59 shows a device containing not only the above fault, but 
adds fallacy as well. The intention is to show alternative means for 
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providing chain isolations by the use of constant current sources, such 
as pentodes, operating in parallel or, as actually shown, by the use of 
series connected secondaries of transformers 13 and 14 working into 
a relatively infinite output impedance. 

The implications to multiple reflections are obvious and the conse- 
quence to the current sensing resistor 19 is disastrous. If, for example, 
delay line 7 passes through a half wavelength, 19 is open circuited 
and develops no sense voltage. The paralleling of pentode chains is 
equally disastrous since one of the pentodes, if viewed looking back 
into the delay line, may prove to be anything but a current source. 

There are no considerations by Schmidt-Briicken to use his device 
to transfer power. All considerations of the relevant requirement of 
the output combiner are missing, as are considerations of output im- 
pedance. Pointedly, we differ with Schmidt-Briicken in calling for an 
output bridge or directional coupler, since these must be power wasting 
four-ports. We have expressly demanded amplifier interactions in an 
output three-port to conserve power. In like fashion we disagree with 
Schmidt-Briicken in demanding an open circuit output to avoid inter- 
action when he does, indeed, chance to use an output three-port. 
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Effects of Precipitation on Propagation at 


0.63, 3.5, and 10.6 Microns 
By T. S. CHU and D. C. HOGG 


The attenuation and scattering of laser beams by rain, fog, and snow 
have been calculated and measured at 0.68, 3.5, and 10.6y. Attenuation 
of the infrared wavelengths by light fog is up to one order of magnitude 
less than at 0.63y. But for dense fog, calculation shows that the attenuation 
at 10.6u can exceed 40 dB per km. It ts found that attenuation by rain 
can be calculated to good accuracy from average path rain rates provided 
that forward scattering is taken into account; this scattering reduces the 
attenuation. 

Measurements of propagation through precipitation over a 2.6 km path 
are discussed in detail and are found to be consistent with predictions. 
The wavelength dependence of attenuation is found to vary from one fog 
to another because of different drop-size distributions. Attenuation of 0.63y 
an rain showers is about 20 per cent less than at 3.5u and ts, of course, 
much less than the attenuation in fog. But even for extremely heavy rain 
showers, the 0.63 attenuation never exceeded 20 dB per km, which ts 
less than the attenuation of millimeter waves under such conditions. Both 
the attenuation and foreward scattering properties of snow appear to be 
between those of fog and rain. 


I. INTRODUCTION 


It is well known that optical waves suffer much attenuation in 
propagation through precipitation in the atmosphere. Penetration of 
light through a dense fog is much more difficult than through a heavy 
shower. One would also expect similar phenomena for near infrared 
waves. However, quantitative data are necessary in order to evaluate 
the potential use of lasers for short range communication through 
the atmosphere. Attenuation by precipitation is the dominant factor 
in determining the feasibility of open-air laser communication sys- 
tems. The interesting fluctuations of laser beams caused by turbulence 
in the clear atmosphere are not discussed here. 
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The gas lasers which provide significant amounts of power within 
the infrared windows of the clear atmosphere are 3.5» and 10.6z; 
3.39 is not suitable because it is attenuated by the small amount of 
methane in the air.t The wavelengths of visible gas lasers with suitable 
output power are 0.49» and 0.63. Since the attenuation of optical 
waves in a fog is largely the result of scattering by numerous small 
drops, it is not likely that transmission of 0.492 and 0.63» through a 
fog would differ appreciably. The literature reveals that the attenua- 
tion of blue light (~0.49.) through a light fog is somewhat greater 
than that of red light (~0.63), while their attenuations in a dense 
fog differ little from one another.’ In this study, 0.63» is the visible 
wavelength used. 

Some computed results of attenuation by precipitation at the wave- 
lengths of 0.63, 3.54, and 10.6 are given first. The effect of forward 
scattering is examined using a Gaussian model for the scattering pat- 
tern of a rain drop and a single scattering model is used to describe 
the broadening of a laser beam by precipitation. Then the results of 
measurements of attenuation by precipitation over a 2.6 km path at 
these three wavelengths are described. Some of the measured data 
have been discussed briefly in previous notes.* * 

The measured attenuation by rain is compared with that calculated 
from the rainfall measured on the path. Optical attenuation by fog 
is compared with total liquid water content derived from path-loss 
measurements taken on a radar operating at a wavelength of 3 mm. 
Similar quantitative correlations have not been made for attenuation 
by snow because of the lack of an independent physical measure- 
ment for this form of precipitation. 


Il. ATTENUATION CALCULATIONS 


2.1 Method of Computation 
The total power removed from a laser beam by a precipitation par- 
ticle, say a spherical drop of radius a, is given by 
E = SC(a, d) (1) 


where S is the magnitude of the incident Poynting vector, and C(a, X) 
is the total or extinction cross section of the drop. The change in S 
upon passing through a layer of precipitation Az is 


AS = —S hz [ n(a)C(a, d) da (2) 
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where n(a)da is the number of drops per unit volume within the range 
(a,a+da). Integrating equation (2) yields 


S a See S48" (3) 


where 


ae / CR ere | (4) 


The total cross section is often normalized with respect to the geo- 
metrical cross section 


Q(a,») = “2. (5) 


Since the energy removed is partly scattered and partly absorbed by 
the particle, the total cross section is divided into scattering and 
absorption cross sections: 


Qext = Qaecat oe Qavs . (6) 


The Q’s are sometimes called efficiency factors. Notice that part of 
the energy represented by the scattering cross section is scattered 
forward; this leads to a decrease in attenuation as discussed in Section 
IIT. . 

The attention coefficient in equation (4) may be calculated if the 
particle-size distribution and the cross sections are given. Equation 
(4) gives the extinction, scattering, or absorption coefficients when the 
appropriate cross section is substituted.* We take the shape of the 
water drops in fog and rain to be spherical. With this assumption, it 
is not a prohibitive task to use Mie’s rigorous expressions to calculate 
the extinction and absorption cross sections of the drops. However, 
very simple approximate formulas are quite applicable here in view 
of the uncertain drop-size distributions and the doubtful sphericity. 
Using a combination of geometrical optics and Huygen’s principle, 
Van de Hulst® obtained the following formulas for the extinction and 
absorption cross section: 


cos 8 


Qext = 2 — 4e° woa( 208.8) sin (p — 8) 





ade" wea(oos8 is ) cos (p — 28) + ( e) toe 2e> ta) 


_ * Extinction is often used in optical terminology for attenuation; however, it 
is misleading. It is more exact to designate the extinction coefficient as the at- 
tenuation coefficient without forward scattering corrections. 
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—4Azn’ —4Azrn’ 
e — | 
Que = 1+ San’ + 3 (7b) 
where 
4ra n’ 2ra 
p= = 1), tan 8 =~: ak ae 


and n and 7’ are the real and imaginary parts of the refractive index, 
respectively. The above two expressions are derived under the assump- 
tions (n — 1) <1, n’ <1, 2ra/d > 1. Numerical comparison indicates” 
that the essential features of the cross section curves using equation (7) 
remain correct for | — 1 — in’ | S 1 and 2ra/y 2 1. When p <1, 
equation (7) reduces to 


Qexe = 3p tan B + 3p°(1 — tan’ 8) (8a) 
Qavs = $p tan B ne p. tan” B. (8b) 
2.2 Refractive Indices 


Figure 1 shows the complex refractive index (n — in’) of liquid 
water in the visible and near infrared region. For the wavelengths 
below 18» the data have been taken from the tabulation by Centeno® 
except for those of n’ in the 2.8u band where Centeno’s values are too 
low and have been replaced by the Bayly results." No serious changes 
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Fig. 1— Complex refractive index (n — tn’) of liquid water, 
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at other wavelengths are reported in more recent literature. In par- 
ticular, for the three wavelengths we investigated, the complex re- 
fractive indices of liquid water are well agreed upon: 


\ = 0.63u, p=1338; n’~ 10“ 
A= B.0n, n = 1.42, n’ = 0.018 
A = 10.6n, n= 1.18, n’ = 0.08. 


The curves beyond 18 in Figure 1 have been obtained from Kislov- 
skii.2 Most measurements of the refractive index of water have been 
conducted at 20°C, but it is known that as the temperature of the 
water is increased, the position of maximum absorption is shifted to 
slightly shorter wavelengths.® 7° However, detailed information on 
temperature dependence seems to be lacking. Much less data have 
been published for the complex refractive index of ice than for water. 
Figure 2, ice at —10°C, is reproduced from Kislovskii for the sake of 
completeness. 


2.3 Cross Sections and Attenuation Coefficients for Monodisperse Drop- 
Size Distribution 


Using equations (7) and the above values of the complex refractive 
indices, extinction and absorption cross sections have been computed 
as a function of radius of water spheres for the three wavelengths 
we mentioned. The data are shown in Figure 3. When the radius of 
the water sphere becomes large, say 1 mm (that is, about the size of 
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Fig. 2— Complex refractive index (n—in’) of ice. 
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Fig. 3— Approximate extinction and absorption cross sections. 


a rain drop), the extinction cross sections approach two, while the 
absorptive cross sections at 3.54 and 10.64 approach unity. The ab- 
sorption cross section of 0.63. does not appear in Figure 3 because 
pure water is practically lossless at visible wavelengths. The extinc- 
tion cross section for 3.54 remains substantially less than that of 0.63 
only when the radius of water sphere is smaller than 1p; also, in that 
range, the extinction cross section of 10.6% almost entirely results from 
absorption. The longer wavelength and the lower refractive index push 
the first maximum of the 10.6, extinction cross section to an “a” of 
about 19. 

Although fogs and rains never have drops of a single radius, it is 
instructive to calculate the attenuation coefficients for this ideal case. 
If we substitute 


w 

= ths C 

n(a) aa? 5(a) (9) 
into equation (4), where w is the liquid water density and 8(a) is the 
Dirac delta function, then the extinction coefficient of liquid water 


becomes 











_ 3 Qexe nepers / Mi gta /™s 
~~ 4 a km my: Pe a km/ m?° (10) 
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where the radius of water sphere, “a” is in microns. The extinction 
and absorption coefficients obtained using equation (10) are plotted 
in Figure 4. 

The vast difference in attenuation by rain and fog is very clear 
from equation (10). Although the liquid water content of a typical 
heavy shower, say 1000 mg/m’, is ten times that of a typical dense 
fog, say 100 mg/m, the fog drop radius is about one thousandth that 
of a rain drop. This is the primary reason why attenuation by rain 
is two orders of magnitude less than that of fog at optical wave- 
lengths. For rains with drop radii about 1mm, the extinction coef- 
ficient is only 0.0065 (dB/km)/(mg/m*) (see Figure 4), thus a 
fairly heavy shower with a liquid water content of 1 g/m* (rain rate 
about 25 mm per hour) has an extinction of about 6.5 dB/km for 
both optical and near-infrared wavelengths. 

At least half of the attenuation resulting from rain at 3.54 and 








60 { 
7773S OL EXTINCTION 
0.63 ph 10.6 ph 
Gs EXTINCTION ie. 
myn Lop) To) 
= | & o = 
a — 


de. 
km 

dB 
km 











= 
° 
ke 
a 
(aa 
Oo 
m9) 
o 
< 
fay 
Pol 
< 
Zz 
Oo 
ke 
O 
Zz 
| ed 
x 
in 


COEFFICIENTS ( 


YX 
| SX = 
. SS 


i | =, 
0.2 Pee 


———— : 
W135 4 ABSORPTION 
er a ea 
“0 2 4 6 8 10 12 14 16 18 20 
RADIUS OF WATER’ SPHERES IN MICRONS 














EXTINCTION AND ABSORPTION COEFFICIENTS ( 














_Fig. 4— Extinction and absorption coefficients of precipitation with mono- 
disperse drop-size distribution. 
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10.6 is absorptive loss while at 0.63,» it is caused almost entirely by 
scattering, as seen from Figure 4. 

Attenuation by fog is contributed largely by water drops less than 
a few microns in radius, therefore much of the 10.6n attenuation is 
absorptive loss, as implied by equation (8) and Figure 4. The extinc- 
tion coefficient of fog is expected to be substantially less (up to one 
order of magnitude) at 10.6 than at 0.63 where scattering predomi- 
nates. However, heavy fogs with liquid water contents exceeding 
0.1gm/m* will attenuate 10.6% more than 40 dB/km. 

Attenuation by fog at 3.5n, largely caused by scattering, remains 
significantly less than that of 0.63» only if the fog drop radii are 
smaller than Ip, as is also true for haze. Fogs which consist of drops 
with an average radius of 5u are expected to attenuate about the 
same amount at optical and shorter near-infrared wavelengths. 


2.4 Calculations for Realistic Drop Size Distributions 


Experimental difficulty and natural variations have caused gross 
inconsistency among the published data on measured fog drop-size 
distributions. Most of the measuring instruments fail to respond to 
droplets less than one micron in diameter, yet there are indications 
that the number density continues to increase as the diameter of 
fog drops decreases toward one micron. The relatively consistent 
information seems simply to be that the number density decreases 
for the larger drop radii. 

For example, Eldridge’s measurement shows that the number den- 
sity decreases approximately with the inverse square of the drop 
diameter.7* In searching for a realistic drop-size distribution for fog, 
the following model?” for the number density, appeared to be plausible: 


n(“) = a(#)" exp | -a()" |. (11) 


Here 
(atl)/y¥ 
4 = (2) ieee: cae Nr me 
Y (2 + 1) Qin Y 
T ee A 
Y 


N is the number of drops per unit volume (say cm®), a, the radius 
of the drops with the maximum number density, and a and y, param- 
eters which can be adjusted to fit an observed fog particle size dis- 
tribution. The liquid water density appropriate to the above model 
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is easily obtained by integrating equation (11). 


ret + 4) 
, ay oe da - Y 3 
lv = nana da = > NV (12) 
0 3 \a r(2 +1) 
Y 


In Figure 5 four plots of equation (11) are given for several indi- 
cated combinations of parameters. Three of these exhibit a slope on 
the larger drop side of the maximum (a > a,,) similar to that of 
Eldridge’s measured distribution. The other combination of param- 
eters (a = 2, y = 3) shown in Figure 5 corresponds to the asymptotic 
distribution of coagulating particles*® and has a sharp decreasing rate 
on both sides of the maximum number density. 

Using these model distributions, extinction coefficients have been 
plotted versus the radius at maximum number density, a,,, in Figure 
6. These curves merge at somewhat smaller values of a, than do the 
corresponding curves for constant drop radius in Figure 4 because 
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Fig. 5— Drop size distribution of fog models. n(r) = Artexp[(— a/y)r"] A = 
(N/am) (a/y)P” y/Tla+1)/y], and an = radius of drops with maximum 
number density. 
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the drops larger than a,, contribute more to the extinction than do 
smaller ones.* The similarity between Figures 6(a) and 6(c) 1s 
attributed to the similar average slopes on the larger drop side of 
the corresponding distributions in Figure 5. 

It is interesting to compare the above model with the Junge dis- 
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Fig. 6 — Extinction of coefficients for various fog models. 


tribution!* which states that n(a) = ca? where 2 < p < 4.1 1° This 
simple relation points out the importance of decreasing number den- 
sity for the larger drops; however, more parameters, such as those 
offered in equation (11), are needed to describe both ends of the 

* Although more extinction is obtained for the same amount of liquid water 


in the form of smaller drops than in the form of larger drops, the extinction cross 
section of a larger drop is, of course, greater than that of a smaller drop. 
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distribution. Since the extinction cross section is of the order of the 
geometrical cross section 7a?, the decreasing rate of larger drops can 
be approximated by the inverse square power (p ~ 2) only in a 
limited range beyond which the number density must decrease much 
more rapidly to insure convergence of the integral in equation (4). 

Raindrop size distributions measured at ground level seem to be 
fairly consistent among the published data. The Laws and Parsons** 
distributions which have been virtually confirmed by later workers 
are il‘ustrated in Figure 7 and are readily available in tabulations.’® 
Substituting these well-known distributions into equation (4) yields 
the extinction coefficient versus rain rate shown as the upper solid 
curve in Figure 8. Notice that the total extinction is the same for 
the visible and near-infrared wavelengths. Half of the extinction 
coefficient at 3.54 or 10.6% is contributed by the absorptive loss as 
illustrated in Figure 8. 


III. FORWARD SCATTERING 


The energy intercepted by precipitation particles is partially scat- 
tered forward, toward the receiver, and this reduces the attenuation. 
Here is an approximate calculation for predicting the reduction of 
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Fig. 7 — Raindrop size distributions by Laws and Parsons. 
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attenuation at 0.63 and 3.54 by forward scattering during rain, 
and a single-scattering model for beam broadening. 


3.1 Approximate Correction Factor for Attenuation by Rain 


The forward scattering pattern of a rain drop is assumed to be a 
Gaussian beam with an effective aperture area equal to the geo- 
metrical cross section of the water sphere.*® Using the geometrical 
parameters illustrated in Figure 9, the scattering pattern of a rain- 
drop of radius a situated at (p, 9, Z) is taken to be 


a eek t/we? 
3 


wr & wif + (Mi 2) | 


and the effective aperture area is 2rw?, = ma’. The beam from the 
transmitter, 7’, is taken to be 


where 


—2p2/w? 
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Fig. 8 — Rain attenuation coefficient for a 2.6 km path. 
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Fig. 9— Geometry for forward scattering. 


w = wie| + (24) | 
TW o7 


and W r is the minimum beam radius, located at the uniphase trans- 
mitting aperture. Then the attenuation in a path of length ZL is: 


where 


L 
al(t — 8) = | a dz} — Gt 
TWr 


fr =e Ce ae eRe" orp dp |: (13) 
0 rw TW 


a is the extinction coefficient (the attenuation coefficient without 
forward scattering correction), 6 the reduction factor resulting from 
forward scattering, Q.../Q-x. the fraction of the intercepted energy 
available for forward scattering, 2/rw7 the maximum flux of the 
unattenuated transmitted beam at the receiver 


(ut = whl + Zz) )). 


and p,/(L — Z) = 6,4 the acceptance angle of the receiver. It is easily 
shown from the geometry of Figure 9 that R/p = L cos 6/Z; thus 
when @ is small, as it is in practice, R = pL/Z. If the acceptance angle 
of the receiver is large compared with the scattering beamwidth of 
the rain drop, ps can be replaced by ~. 

The above approximations simplify evaluation of the inner integral 
in (13) and lead to the following expression for the forward scattering 
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correction factor: 





: L 2 
B es : Oa / Wr dZ. (14) 
os 0 


Using equation (14), & versus raindrop radius has been computed for 
L = 2.6 km, \ = 3.5u, wor = 0.55 em, and \ = 0.63n, wo7 = 0.25 cm. 
Table I gives the data. The two values of w,7 are the clear weather 
beamwidths in our experiments. Upon numerical integration over the 
raindrop size distribution one obtains the dashed curves in Figure 8. 
The above approximate calculation is essentially a crude multiple 
scattering model for very narrow beams because equation (13) implies 
that the forward scattering energy effectively rejoins the original 
transmitting beam. 


3.2 A Single Scattering Model for Beam Broadening 


Measurements discussed later show that various degrees of beam 
broadening occur when a 0.63» laser beam propagates through precip- 
itation. Beam broadening resulting from rain is more pronounced 
than that of snow and fog. Let us consider a single scattering model 
describing the effect of forward scattering on the attenuation, and 
beam broadening of a narrow laser beam. We are interested in the 
power distribution at the receiving plane. 

Using cylindrical coordinates with origin at JT as illustrated in 


Taste I—Vorwarp ScaTrerRiInGc Correcrion Facror B FoR A 
2.6 KM PROPAGATION PATH 





Raindrop »X = 0.63z A = 3.5u 
Radius (cm) wor = 0.25cem wey = 0.55em 
0.025 0.106 0.025 
0.050 0.202 0.049 
0.075 0.289 0.071 
0.100 0.366 0.093 
0.125 0.434 0.1138 
0.150 0.495 0.13838 
0.175 0.548 0.152 
0.200 0.595 0.169 
0.225 0.636 0.186 
0.250 0.672 0.202 
0.275 0.7038 0.217 
0.300 0.730 0.232 


Wr is the beam radius at which the field amplitude has fallen to 1/e 
of its maximum value in the collimated transmitting aperture. 
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Figure 9, the intersection of 7P with the receiving plane is (pL/Z, 
¢, L) whereas an arbitrary point in the receiving plane is designated 
by (r, ¢’, ZL). The Cartesian coordinates of these two points in the 
receiving plane are (pL/Z cos ¢, pL/Z sin ¢, L) and (r cos ¢’, r sin 
¢’, L) respectively. Since 6 is always very small for narrow laser 
beams, the scattering pattern of a raindrop will show negligible obliq- 
uity at the receiving plane and may be written as 


BE gi) @ _ \ 
5-2, exp{ 2 | i; cos ¢ — 7 cos@ 
pL ] es, ee ieee \ 
+ (2! sin @ rsing 


| (22 B BPE eels 4 
Ree es dda 


TW, 


(15) 








2 2 
Ws We 


Then the scattering by a uniform rain of constant drop radius a over 
a path of length Z is: 


= [a8 Qren gy) 1 ah ie <5 oo" Sp inde 06 


The definitions of the parameters in equations 15 and 16 are the 
same as in the preceding section. The integration over ¢ is recognized 
as a Bessel function of imaginary argument 


A BES cont 2) 


[ex |—4° do = ae 2). 
Using the well known identity 
1 2 
- J,(at) exp (—p’P)t dt = Op exp (- “,) (18) 


where a is an arbitrary complex number, one can perform the p-inte- 
gration in equation 16: 


. -4oLr 1 yas 
2m [ J s(-i sel) exp pe Hs) |e dp 


Bee Vag (ee 
~ 2 \Bw, + wl? 8 \Zw) \Putw) | 
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Finally the Z-integration can be reduced to the following form for 
numerical computation: 





Qaca r Wr 27° 
BS | a | ae dZ. (20) 
ut + (3)u w+ (Z) 


Both a and w, in the above equation are functions of raindrop radius. 
Upon numerical integration of P; over the raindrop size distribution 
one obtains the scattering pattern in the receiving plane. 

Using L = 2.6 km, »X = 0.63n, Wor = 0.25 cm (equivalent to our 
experimental clear weather half-power beamwidth of twenty seconds) 
and the Laws-Parsons distribution, scattering patterns for a light 
rain 2.5 mm per hour and a heavy rain 50 mm per hour have been 
computed; these are shown in Figures 10(a) and 10(b), respec- 
tively. The direct beam pattern, the total power pattern, and the 
normalized total power pattern are also shown there. It has been 
tacitly assumed in the computations that the scattering and the direct 
beam are attenuated by the same factor, e-“”, and that the total 
power pattern is the arithmetic sum of the direct beam and the scat- 
tering. For light rain, (Figure 10a), the scattering modifies only 













































































10 
_-TOTAL POWER PATTERN 
5 |—,.4—— 
5 
2 
1.0 ) 
a 
% © 0.5 
3 |, NORMALIZED TOTAL _ = ig 
& POWER PATTERN 2 sf 
0.2 Ww 
Ww w a 
= > oO 
= F OM -10 Ww 
< : ~ < fa) 
o PATTERN oy: NORMALIZED 
TOTAL POWER -4-15 
fe be PATTERN 
0.01 -20 
es ween & = 
0.005 | SiRECT BEAM 0.005 |S iRECT BEAM 
PATTERN ~~ PATTERN —-~ -25 
0.002 + Wor = 0.25 em 0.002 |- Wor =0.25 cm 
0.001 0.001 | -30 
) 10 20 30 40 50 ) 10 20. 30 40 50 
SECONDS SECONDS 


Fig. 10— Calculated scattering pattern and beam broadening. (a) Light rain, 
2.5 mm per hour. (b) Heavy rain, 50 mm per hour. 
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slightly the clear weather or “direct” pattern. However, in heavy rain, 
(Figure 10b), the scattering becomes predominant. 

One might expect more beam broadening the heavier the rain; 
however, since the average drop size increases with rain rate, the 
scattering pattern of a heavy rain is narrower than that of a light 
rain. This effect tends to offset somewhat the beam broadening by 
heavy rain. On the other hand, side lobes of rain drop scattering 
which have been neglected in the Gaussian approximation and mul- 
tiple scattering effects which are inherently absent in the single 
scattering model will tend to increase the beam broadening. Also 
notice that beam broadening will be quite sensitive to inhomogeneity 
of the rain. 


IV. EXPERIMENTAL EQUIPMENT 


Figure 11 is a block diagram of the transmitting and receiving 
equipment for 0.63 and 3.5u. Care was taken to stabilize the trans- 
mitters both thermally and mechanically. The equipment for 10.6. 
has been described by Wilson and Penzias.* Figure 12 gives the 
profile of the 2.6 km propagation path at Holmdel, N. J. Table II 
lists the parameters of the transmitting and receiving systems. 


4.1 0.6384 Equipment 


The 0.63 transmitter is a 50 mW multimode He-Ne laser; the 
cavity consists of two 10m-radius of curvature multilayer dielectric 
mirrors separated by 1.25m. The output is chopped mechanically at 
330 cycles per second and passed through a refracting telescope which 
has an input lens of 5 cm focal length and an output lens of 50 cm 
focal length; the transmitting beam has a half power beamwidth of 
about 20 seconds in clear weather. 

The received signal passes through attenuators, a 3A bandwidth 
filter, a double-convex lens with an iris in the focal plane, and a 
polarizer to a photomultiplier. The receiver acceptance angle, about 
10 minutes, is determined by the iris opening. Using a reference 
signal transmitted by telephone line, the synchronous detector gives 
a signal to noise ratio of 65 dB on a clear day; the limit of this 
sensitivity is set by the scattering of sunlight. 


4.2 3.54 Lquipment 


The 3.5y transmitter is a 5 mW multimode He-Xe laser; the 1.1m- 
long cavity consists of a plane reflecting mirror and a calcite polarizer 
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Fig. 11— Sketeh of the equipment. 
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Fig. 12 — Profile of transmission path. 


which acts as an output mirror. The output passes through the 330 
Hz chopper to a bisected nearfield Cassegrainian telescope, a com- 
bination of two bisected confocal paraboloids. The main reflector is 
one half of a 38.1 cm diameter spin-cast paraboloid of 30.48 cm focal 
length while the subreflector has a diameter of 1.27 em and a focal 
length of 1.02 em. 

Another bisected near-field telescope collects the signal for the 
receiver. The concentrated beam from the subreflector is amplified by 
two He-Xe tubes, each of about 18 dB gain. The good impedance 
match of the telescope plays an important role in successful applica- 
tion of these high-gain laser amplifiers. After the two amplifier tubes, 


TasLE IJ — PARAMETERS OF THE TRANSMITTING 
AND RECEIVING SYSTEMS 








Wavelength (n) 0.63 3.5 10.6 
Laser oscillator He-Ne He-Xe CO, 
Transmitter power (mW) 50 5 500 
Transmitting beamwidth (rad) 10-4 2 xX 10-4 10-3 
Receiving aperture diameter (cm) 5 12.5 15 
Laser amplifier None He-Xe None 
Detector Photo- In-Sb (77°K) — Thermopile 
multiplier 


Measuring range (dB) 65 45 50 
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the signal passes through a calcite polarizer, a quartz lens, and an iris, 
to a liquid-nitrogen-cooled indium-antimonide photovoltaic infrared 
detector whose output is fed into a synchronous detection system. The 
receiver acceptance angle, about 1 minute of arc, is determined by the 
dimensions of the amplifier tubes and the telescope. A signal-to-noise 
ratio of 45 dB is obtained on a clear day. The receiver sensitivity is 
limited by the laser amplifier noise. . 


4.3 10.6n Equipment 


A COz laser with an output of 0.5 W and a thermopile at the focal 
plane of a 15 cm collecting mirror comprise the 10.6 transmitting 
and receiving systems. The laser is operated with a continuous flow 
of gas. Therefore, virtually uninterrupted measurement has been pos- 
sible for a period of more than one year. 


4.43 mm Wavelength Radar 


The 3mm wavelength radar, based on a method used earlier for 
study of absorption by atmospheric gases at millimeter wavelengths,”° 
provides an independent measure of the liquid water content of fogs. 
The system uses a klystron, modulated linearly in frequency. Two 
precisely constructed trihedral corner reflectors are placed on the 
propagation path, one near the center of the path and the other at the 
optical transmitter terminal (see Figure 12). Measurement of the 
relative power reflected from the two reflectors results in attenuation 
data, accurate to about 0.2 dB per km, for that portion of the path 
between the reflectors. 


4.5 Rain Gauge System 


Three gauges of the tipping-bucket type which register a pulse for 
each 0.25 mm of rain, are located at the transmitting site, path mid- 
point and the receiving site (see Figure 12). Data from the three 
gauges are recorded at the site of the optical receivers. No quantita- 
tive measurements are made during snow. 


V. MEASURED RESULTS 


5.1 Fog 


In spite of available signal-to-noise ratios of 65 dB for 0.63, 45 dB 
for 3.5n, and 50 dB for 10.6n, the received signal at all three wave- 
lengths has often fallen below their noise levels during dense fog. Thus 
our ability to measure the wavelength dependence of attenuation by 
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fog is limited to relatively light fogs with attenuations less than 25 
dB per km. Figure 13(a) shows that in a typical morning fog (Sep- 
tember 21, 1965) the attenuation (in decibels) of the 3.54% wave is 
about half that of the 0.634 wave. This result may be explained by 
fog droplets of about 1p diameter. 

The measured data occasionally reveal a much smaller ratio as in 
Figure 13(b) where the 3.54 attenuation (September 15, 1965) 
amounts to only one sixth of the 0.63» attenuation; this result implies 
a fog with many droplets smaller than one micron. 

The measured attenuation of 10.6% in a light fog is lower than at 
3.54, as shown in Figure 13(c). At the beginning of this record (1.30 
p.m. January 4, 1967), the signal at 0.63 and 3.5 are equal at a 
level of about —15 dB while the 10.6 signal is at about —2 dB. 
This result can be explained, for example, by entering Figure 6(d) 
at a maximum drop radius of 1.34; there, one sees that the attenua- 
tions at 0.63» and 3.54 are equal and exceed the 10.6n attenuation by 
a factor of eight, roughly the observed value. But at 3 p.m. in Figure 
13(c), the situation has changed significantly, and the ratios of the 
0.63, 3.54, and 10.6» attenuations are now about 15:6:1. However, if 
we enter Figure 6(d) at a maximum drop radius of 0.7», the attenua- 
tion ratios are found to be very close to the observed values; thus it 
is found that small changes in the drop sizes can drastically affect 
wavelength dependence. 

A percentage time distribution for the ratio of attenuation by light 
fog at 3.54 and 0.63, is plotted in Figure 14. The data, taken during 
seven separate fogs, were accumulated over sixteen hours during 
which the attenuation at either wavelength over the 2.6 km path was 
between 5 and 50 dB. It is seen that the median value of the ratio is 
0.62 and the attenuation at 3.5y is less than at 0.632 85 per cent of 
the time. The 15 per cent of the time that the attenuation at 3.5p ex- 
ceeds that at 0.634 may be caused in part by equipment instability; 
but it can be real for certain fog parameters, such as, in Figure 6(d), 
a,, = 2 microns. 

The data of Figures 4 and 6 predict that attenuation of 10.6 is 
insensitive to the drop-size distribution in light fog where drop radii 
are about 1p. Therefore, one expects close correlation between the 
liquid-water density of a light fog and the measured 10.6n attenua- 
tion. 

During measurement of attenuation of the laser beams by a morn- 
ing fog May 20, 1966 (Figure 15), the liquid water density of the 
fog was determined by measuring the attenuation of 3 mm waves. At 
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Fig. 14— A distribution for the ratio of attenuation by light fog at 3.5u and 
0.63. 


3 mm, where the wavelength is certainly much greater than the fog- 
drop radius, the liquid water content is given by M = 0.25y gm/m?, * 
y being the 3 mm attenuation coefficient in dB/km. 

As shown in Figure 15, at 9:30 a.m. the laser beams could not be 
detected and the 3 mm data indicate a fog density of 0.8 gm/m?. As 
the fog gradually lifted, the optical signals emerge in sequence above 
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Fig. 15— Correlation between attenuation of laser beams and liquid water 
content measured by a 3 mm wavelength radar in fog. 
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the noise level but because dense haze lingered throughout the day, 
the 0.63 wave remained below the clear weather level. There is good 
correlation between the attenuation of the laser beams and the liquid 
water content of the fog obtained by the 3 mm wavelength measure- 
ment. The zero reference of the liquid water density is taken at the 
vanishing of the attenuation of the 10.6 laser beam. However, precise 
agreement between the measured attenuation of 10.6» and that calcu- 
lated from the liquid water content obtained at 3 mm is difficult to 
establish because of the very small attenuation at 3 mm for light fogs. 
It should also be kept in mind that the fog density measured by the 
3 mm radar is an average over only about half of the propagation 
path (see Figure 12). 


5.2 Rain 


Transmission through some light rains is difficult to interpret because 
the rain is mixed with fog; however, summer showers are very suitable 
for evaluating attenuation by rain. Figure 16 shows typical measured 
data during a summer shower; attenuation of 3.5u is consistently 
about 20 per cent greater than that of 0.63». The attenuation is well 
correlated with the rain rate averaged over the path. The average 
path rain rate is taken to be the arithmetic average of the readings 
of the three rain gauges described in Section 4.5, with that of the 
midpoint gauge carrying a double weight. 

Owing to the variation in raindrop terminal velocity and drop size 
distribution, the relationship between the attenuation coefficient and 
the rain rate can be hardly expected to be linear. However, the pre- 
dicted attenuation coefficient for the rain rate between 12.5 and 100 
mm per hour is almost linear as shown in Figure 8. Therefore, it is 
justifiable to determine a least square line for measured data of at- 
tenuation by rain above a rate of 12.5 mm per hour. Figure 17 shows 


a least square line 
B 
o( 88 | - 0.1559(™) + 2.66 


for the measured 0.63» data of Figure 16 and two other showers. 
The calculated curves for 0.63» attenuation by rain with and with- 
out forward scattering correction (from Figure 8) have been included 
in Figure 17 for comparison. The good agreement between the meas- 
ured line and the curve calculated with forward scattering may appear 
somewhat fortuitous in view of the rms deviation 2.8 dB/km. How- 
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Fig. 16— Measurement of 2.6 km transmission loss in rain showers on Septem- 
ber 24, 1965. 0 dB-signal level in clear weather. 


ever, only peak attenuations of summer showers which are relatively 
free of foggy conditions have been taken as measured data in Fig- 
ure 17. 

The relatively large rms deviation is explained as follows: the sam- 
pling of a spatially nonuniform shower by three rain gauges does not 
provide a satisfactory average rain rate over the 2.6 km path; the 
attenuation depends upon drop size distribution which may vary from 
one shower to another; since optical attenuation by fog is orders of 
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magnitude higher than that by rain, any slight mist or haze within 
the rain can appreciably change the measured values. 

Simultancous measurements at all three wavelengths were carried 
out in an unusually heavy rainstorm of September 21, 1966, during 
which more than 125 mm of rain accumulated. A sample of the 
measured data from this storm is given in Figure 18. 
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Fig. 17— Measured rain attenuation coefficients of 0.634 for a 2.6 km path; 
standard deviation is +2.8 dB per km at a rain rate of 50 mm/hr. 


Points corresponding to the peak attenuations of the separate 
showers in the storm are plotted in Figure 19 and least-square lines 
for each wavelength have been drawn for rain rates between 10 and 
60 mm per hour. The attenuation by rain at 10.6 and 3.5 is con- 
sistently greater than that at 0.63. At 0.63 the agreement between 
calculated and .measured values is again rather good. Notice that the 
measured points at 95 mm per hour fall below the predicted curve.* 
When we compare Figure 19 with Figure 13(c), say, it is obvious 
that the wavelength dependence in propagation through rain is in- 
verse to that through fog. 
~-* The 0.63n data in Figures 18 and 19 were not taken with the same receiving 
telescope as those of Figure 16 (5 cm diameter, see Table 1). Rather, a near-field 
cassegrain telescope with a collecting aperture of 600 cm was used. This change 
was made to determine whether or not a large receiving aperture would intercept 
more of the energy forward scattered by the rain (see Figures 10 and 21a) and 
thereby further decrease the attenuation. However, the attenuations (relative 


to the clear-day value) measured with the large aperture did not turn out to 
be significantly less than those measured with the small one. 
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Another sample of measured data is shown in Figure 20;7 it is of 
particular interest because the attenuation at 0.63 is much lower 
(about 18 dB) than in previous measurements and calculation for the 
same rain rates. The attenuation at 3.5u and 10.6p is only slightly 
below the prediction. Examination of the raw data shows that the 
rain distribution along the path was characterized by higher rates 


SIGNAL LEVEL IN DECIBELS 





125 


100 
75 
50 


25H 


AVERAGE PATH 
RAIN RATE (mMm/HR) 


0 
1:00 11:30 NOON 12:30 
E.S.T, 


Fig. 18 — Measurement of 2.6 km transmission loss in a rain storm. 0 dB-signal 
level in clear weather. 


near the receiving site. This observation suggests that this very low 
attenuation at 0.63» is caused by a strong forward scattering con- 
tribution which is most effective when the scattering particles are 
near the receiver. This result is also confirmed by the measured data 
of the rain showers in Figure 16. There the relatively low 0.63, at- 
tenuation coefficient measured in the second shower is also associated 
with a rain distribution with high concentration near the receiving 
site. 


+The data in Vigure 20 were taken with a 0.634 transmitter beamwidth of 
5 & 10-* radians. 
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Fig. 19 — Measured attenuation cocfficients of a rain storm for a 2.6 km path. 


5.3 Snow 


Figure 21 shows measurements of transmission loss through a light 
snow shower. The accumulation never exceeded about 12.5 mm per 
hour. Snows heavier than this have caused the signals at all three 
wavelengths to fall below the noise levels. The 10.62 wave often ap- 
pears to suffer greater attenuation in snow as illustrated in Figure 21; 
however, no definite wavelength dependence is consistently observed. 


5.4 A Cumulative Distribution of Attenuation at 10.6p 


In Section 4.3 we mentioned that the 10.6% equipment operated 
continuously. Signal levels for the period June 1966 to May 1967 
have therefore been analyzed and reduced to a cumulative distribution 
of attenuation as shown by the central curve in Figure 22. For ex- 
ample, one sees that attenuation on the 2.6 km path, resulting from 
all causes, exceeded 30 dB (11.5 dB/km) three per cent of the time. 
Also shown in the figure are distributions for the months with the 
highest (March) and lowest (July) outage times. Clearly, during 
July, a system operating at 10.62 would be quite reliable; however, 
during March, heavy attenuations, mainly caused by fog, occur all 
too often on a 2.6 km path. 
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5.5 Beam Broadening at 0.63p 


Typical measured results on the beam broadening effects of rain, 
snow, and fog are presented in Figures 23, 24, and 25, respectively. 
Curve 1 in each figure is a typical clear-weather beam (about 20- 
second half-power beamwidth); this broadening is essentially deter- 
mined by the inhomogeneities of the refractive index in the atmos- 
phere. Notice that the clear weather beam itself is subject to modest 
variation but this variation is generally much less than the beam 
broadening observed under conditions of precipitation. The beam 
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Fig. 20 — Measured attenuation over 2.6 km during a rain shower on July 11, 
1967.0 dB-signal level in clear weather. 
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Fig. 21— Measurement of 2.6 km transmission loss during a light snow shower 
on January 26, 1966. 0 dB-signal level in clear weather. 
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Fig. 22—A cumulative distribution of attenuation at 10.64 over a 26 km 
path. O = Average attenuation from June 1966 to May 1967. A = March 1967, 
highest month. [1] = July 1966, lowest month. 
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Fig. 28 — Beam broadening caused by rain. 
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Fig. 24 — Beam broadening caused by snow. 
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Fig. 25 — Beam broadening caused by fog. 


broadening is measured by moving the receiver on a horizontal track 
transverse to the axis of the incident beam. 

The other measured curves show beam broadening resulting from 
precipitation and the corresponding attenuations along with descrip- 
tions of precipitation; the curves are normalized at 6 = 0, the direc- 
tion of the beam axis. In general, the highest attenuations correspond 
to the greatest beam broadening. In Figure 24, we see that the onset 
of beam broadening caused by snow occurs when the attenuation is 
fairly high. Broadening effects in the case of rain are observable at 
much lower attenuations (Figure 23). This is because raindrops pro- 
duce more forward scattering than snow flakes which suffer from poor 
geometry and random orientations. Thus, the X points show the same 
amount of beam broadening in Figures 23 and 24; but the correspond- 
ing attenuations are 5.4 dB/km for the rain and 15.6 db/km for the 
snow. 

Although the beam broadening, caused by fog also is readily meas- 
urable, the fog drop, being comparable in size with wavelength, is even 
less efficient in forward scattering than the snow flake. The Y points 
in Figures 24 and 25 show that the attenuation by fog is about twice 
that of snow for the same degree of beam broadening. 

The measured beam broadening due to rain exceeds the prediction 
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of a single-scattering gaussian model as discussed in Section 3.2. 
Quantitative correlations are difficult because of the time variability 
and inhomogeneity of rainfall along the path. 


VI. DISCUSSIONS 


A simple picture of the vast difference in attenuation of micron 
waves by fog and rain has been substantiated by theory and experi- 
ment. The relatively low attenuation by rain is, in part, the result of 
forward scattering. Owing to the difficulty of obtaining independent 
physical data on snow, quantitative correlation between physical 
measurements of precipitation and transmission has been made only 
for rain and fog. The predictability of attenuation of optical waves 
by rain is comparable to that for microwaves.”! Indeed, visible wave- 
lengths are attenuated less than millimeter wavelengths. | 

Attenuation by deénse fogs often exceeds 20 dB/km. The wave- 
length dependence of attenuation by light fog, that is, progressively 
less attenuation with increasing wavelength, has been observed as 
expected, but the dependence is a strong function of the drop-size 
distribution. Similar results for incoherent light have been reported 
in the literature.*® 2? Although the predicted wavelength dependence 
for heavy fogs which consist of drops with average radius of 5,, is 
weaker than for light fogs, the attenuation at 10.6% would still be 
significantly less than at visible and shorter infrared wavelengths 
(see Figure 4). However, a typical dense fog of 0.1 g/m® liquid 
water density is calculated to attenuate 10.6% by at least 40 dB/km. 
Attenuation by fog at 10.64 is characterized by strong absorption 
and relative insensitivity to drop-size distribution. 

The measured attenuation by rain, about 0.2 (dB/km)/(mm/hr) 
at 0.63, is consistently less (about 20 to 40 per cent) than at 3.5h 
and 10.6n. Rain drops, which are absorbing spheres at 3.54 and 
10.64, offer less scattering than at visible wavelengths and also 
have a less directional scattering pattern owing to the longer wave- 
lengths. 

The scattering patterns of fog drops have little directivity com- 
pared with the laser beams themselves therefore the forward scat- 
tering correction for the fog attenuation of narrow laser beams is 
small. If the beamwidth of the incident laser beam is comparable 
with the width of the scattering patterns of precipitation particles, the 
forward scattering will reduce the scattering loss significantly which, 
indeed, is mentioned in Section 5.2 for the propagation of 0.63» through 
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rain. Naturally the absorptive loss is not subject to any reduction. 
But the scattering loss of a laser beam propagating through fog can- 
not increase indefinitely as the fog density increases.* In particular, 
for a visible source radiating into a homogeneously distributed dense 
fog, the energy must ultimately diffuse outward in all directions, thus 
the free space loss of an isotropic point source appears to be the 
limit of the scattering loss. 

When the scattering by precipitation particles is added to the at- 
tenuated original beam, a broadened beam is observed at the receiv- 
ing plane; this has been illustrated by a single scattering gaussian 
model in Section 3.2 and by measured results in Section 5.5. 

The apparent attenuation on the beam axis is reduced by the scat- 
tered component. For heavy rain, the received signal is dominated 
by the scattering as illustrated in Figure 10(b) and the reduction in 
attenuation becomes significant. When the scattering correction is 
small, estimates by the single scattering model (Section 3.2) and the 
multiple scattering model (Section 3.1) are practically the same. For 
relatively large corrections the multiple scattering model is expected 
to give a better estimate. 

The precise prediction of attenuation of laser beams by precipita- 
tion in the atmosphere is hampered by the uncertainty of particle 
size distributions. A better knowledge of the drop size distribution 
is certainly desirable to provide a firm ground for the prediction of 
attenuation, but even a measured drop-size distribution might lead to 
ambiguous interpretation when the variability and inhomogeneity of 
the natural fog or rain is unaccounted for. Summations over unreliable 
drop size distributions introduce complicated-looking problems but 
provide little more information than a simple relation such as equa- 
tion (10). 

Although it is convenient to distinguish between rain and fog by a 
classification of drop size, one should be aware of the frequent occur- 
rence of the intermediate states in the form of a drizzle or a mixture of 
rain and fog. The attenuation under these conditions has been observed 
to exhibit transition characteristics; however, little can be said beyond 
this trivial observation. 

The prediction of attenuation by snow has not been touched upon 
here, but the experimental data indicate that it lies intermediate be- 


* A common experience is to observe a heavy fog during daylight hours; al- 
though the sun is completely obliterated, the observer still has considerable 
scattered light in his environs. 
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tween that of rain and fog for a given amount of equivalent liquid 
water. 

In order to evaluate relative merits of systems using different wave- 
lengths, it is interesting to compare the attenuation by precipitation 
at optical wavelengths with that in other parts of the electromagnetic 
spectrum. Figure 26 presents the approximate attenuation coefficients 
for a typical dense fog and a representative shower in the frequency 
range from 10% to 10% cps.1%?1 The absorption for a water layer 
containing the same amount of liquid water as the rain is also shown 
in Figure 26 for reference. The attenuation by rain decreases slowly 
from millimeter to visible wavelengths. However, the attenuation by 
fog is about two orders of magnitude less than that of rain at upper 
microwave frequencies, while the optical attenuation by fog is about 
two orders of magnitude greater than that of rain. Notice that, even 
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Fig. 26— Attenuation due to liquid water in various forms as function of 
frequency. 
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at 10, where attenuation by fog passes through a weak minimum 
(see Figure 26), the attenuation is still rather large as far as a prac- 
tical communication system is concerned. 
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The Stream Guide, a Simple, Low-Loss 
Optical Guiding Medium 
By P. KAISER 


Heat generated by viscous friction in a high speed gas flow in a pipe 
causes a transverse refractive index variation that can continuously focus a 
light beam passing through this pipe. The comparatively simple construction, 
together with the low losses inherent in gaseous guiding media, makes the 
stream guide attractive as an optical beam waveguide. This article pre- 
sents theoretical and experimental results. 


In thermal gas lenses, temperature gradients and ensuing refractive 
index variations are caused by heat transfer from an external heat 
source. The tubular gas lens power is supplied by resistive heating; 
the figure of merit, that is, the focusing strength per unit of applied 
power is a criterion for its efficiency.2 The efficiency may be improved 
by using Peltier heat as it has recently been suggested in connection 
with the proposal for an alternating gradient, thermo-electric beam 
waveguide.® In still another way, heat is generated by viscous dissipa- 
tion which exists in a high-speed gas stream in a pipe. The temperature 
rise creates a refractive index variation capable of continuously focus- 
ing a light beam passing through this pipe. 

The adiabatic flow of a compressible fluid is treated in the litera- 
ture.* > The axial density variation resulting from frictional pressure 
loss and the transverse density variation caused by viscous dissipation, 
together with the assumption of a parabolic transverse velocity and 
temperature profile, lead, after simplification, to the following expres- 
sion for the refractive index n(r, 2): 

2 
nr,2) =1—3z(w — 1) eg Pa M,A1(z). (1) 
a “ Po 
In this equation, a is the radius of the pipe, no the refractive index at 
atmospheric pressure po, x the ratio of the specific heats (c)/cy), R 
(0.90) the so-called recovery factor,* p,, M, and M(z) the pressure 


* See p. 211 of Ref. 4. 
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and the Mach number at the beginning and along the pipe, respec- 
tively. 

Provided that the absolute velocity is small compared with the 
velocity of sound, and that simultaneously the relative change of the 
velocity along the guide also is kept small, it can be shown that the 
eikonal equation yields as ray path r(z): 


re) = 1 + 895] co(22 0 + 694) + (2S. + 69°) (2) 
with 
M, 


$ 

= | — lI) —DR 2] and $8 = ie 
whereby J1/3(v) and Ni;3(x) are the Bessel and Neumann functions 
of fractional order 1/3, C; and Cz are integration constants and f 
represents the average friction coefficient. Equation 2 describes an 
undulating trace with steadily decreasing amplitude and period as the 
result of the ever increasing velocity and, consequently, focusing 
strength downstream. 

If one considers the flow parameters to be locally independent of 
z one obtains for the period of the beam size variation Ly(z) (with 
R = 0.90; see Ref. 6): 


a= 


142) = 8.8a{ (m — Doe — Daa |” (3) 


The 1/, half width of the beam which does not change its size in this 
pseudouniform guide is given by wy = (ALy(2))1/?/m (A = light wave- 
length). One notices that small spot sizes can be achieved by choosing 
a monatomic gas with high refractive index, by selecting a small pipe 
diameter and by increasing either the initial velocity or the pressure, 
or both. However, a characteristic feature of frictional flow is a com- 
paratively large axial velocity gradient for high initial velocities, a 
fact which directly affects the total pipe length employable for a 
specified exit speed. 

The experimental setup is shown in Fig. 1. The gas was injected into 
the pipe through a large number of radial jets. The gas flow and 
pressure were adjusted in such a way that a light beam which en- 
tered the guide collimated (w ,;, ~ 0.9 mm) left it either collimated 
or with a focal point coinciding with the end of the pipe. Thus the 
total pipe length represented an even or odd multiple of L,/2, respec- 
tively. For a pipe length of 3.07m and at air flow rates between 40 
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Fig. 1— Experimental setup for measurements with the stream guide. 


to 55 liters per minute (23.2 < p, < 64.7 psi, po = 14.7 psi) the 
averaged period of the beam size variation ranged between approxi- 
mately 6.14 and 1.53m. Typical values for the Reynolds number were 
in the order of 40,000. 

In Fig. 2 the pressure and flow rate is varied in such a fashion as to 
keep the averaged focusing strength constant (ZL, = length of pipe L). 
For a given focusing strength, by increasing the pressure the required 
velocity decreases and the final velocity approaches more and more 
the initial velocity. The use of argon instead of air increased the focusing 
strength as expected. The corresponding theoretical periods of the beam 
size variation at the beginning and the end of the pipe as computed via 
the experimental pressure and flow rate from (8) are also shown (Jy, 
and L,,, respectively). Their average value Z,, is roughly constant 
and the results for both air and argon are in close agreement. However, 
for L = 6.42m, D,, was approximately 36 percent lower than the ex- 
perimental value. In general, different pipe lengths and varying opera- 
ting conditions resulted in values of Z,, which varied around the experi- 
mental results. 

The axial velocity distribution depends largely on the friction coef- 
ficient f. Its value may be determined theoretically, but owing to its 
involved dependence on various factors, an experimental evaluation 
is more reliable. From Fig. 2 follows that f is smaller for smaller 
absolute velocities. Thus, for comparable focusing action, an increase 
in pressure not only results in smaller velocity gradients because of 
a reduced initial velocity but also because of a reduced friction 
coefficient. 


764 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


0.25 25 
A AIR 


O ARGON 
~ | 20 

















Lp, AND Lb, IN METERS 














MACH NUMBER M; AND Moz 

















6 7 8 
P,/Po (Po =1 ATM) 


Fig. 2— Experimental Mach number and theoretical period of the beam size 
variation Ly, and Ly, (according to equation 3) at the beginning and the end of the 
pipe as a function of the initial pressure p,/po for air and argon (temperature 22°C). 
The experimental period of the beam size variation was constant and corresponded 
to the total length of the pipe (6.42m). 


The eventual application of the stream guide for long distance 
communication requires that the pressure be reestablished periodically, 
whereby the period is expected to be in the range of a few meters. 
This can be accomplished by alternately injecting gas from a high 
pressure supply line and exhausting it into a lower pressure pipe 
running parallel to the guide. 

Focusing properties which are close to those of a uniform guide 
can be achieved if the stream guide is operated at high absolute 
pressures together with small pressure drops along the guide. Rayleigh 
scattering losses impose an upper pressure limit of approximately 500 
psi for a maximum transmission loss of 1 dB per km (for air) if losses 
resulting from turbulent scattering are neglected. Also, eventual mode 
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conversion losses caused by a non-parabolic refractive index varia- 
tion are not taken into account. 

Whether the stream guide has actual applications as a simple, 
continuous guiding medium depends on the results of further studies 
of its optical properties, especially the effects of velocity gradients 
and turbulence on beam fluctuations and power loss. 

I gratefully acknowledge valuable comments by C. A. Fritsch, T. J. 
Kessler, and R. J. Port, about the characteristics of compressible fluid 
flow, and I appreciate the assistance of H. W. Astle in carrying out 
the experiments. 
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Pulse Shuttling in a Half-Mile 
Optical Lens Guide 


By D. GLOGE and W. H. STEIER 


To investigate the distortion of a gaussian light beam in a long lens 
guide a light pulse was shuttled back and forth up to 150 times in a 
straight underground conduit with lenses (quality \/10) positioned con- 
focally every 140 m. The pulse was formed by suddenly ‘‘dumping’’ the 
light stored in a He-Ne laser cavity into the guide. Part of 
the round-trip loss was replenished by a He-Ne amplifier. In scanning the 
intensity in the guide cross-section on various round trips of the pulse it 
was found that the gaussian profile was virtually undisturbed. 


I. INTRODUCTION 


Lens guides with widely separated glass lenses for long distance 
optical communication have been proposed in several varieties* * 
and investigated in various laboratory and field experiments.*> The 
use of dielectric coatings reduces the loss to one per cent per lens 
or less,? an estimate which has been confirmed by the experiments 
reported here. This loss figure permits 1000 lenses to be placed in 
the beam path before the beam intensity is down by 40 to 50 dB 
and a repeater or amplifier is necessary. 

An average distance of 100 m between lenses would lead to a 
repeater spacing of 100 km. There is legitimate doubt that a light 
beam can be kept bundled and centered on the guide axis over this 
distance. Measurements in a 140 m underground air-filled pipe showed 
a slow beam drift of several mm per week correlated with the ground 
temperature field. Thermal insulation of the pipe could reduce this 
effect, but eventually along the line there will have to be some beam 
position control to cope with this and with an expected slow ground 
drift. 

Careful investigation in the same pipe revealed that density fluc- 
tuations or any short time influence of the gas in a well-shielded 
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conduit 5 feet below the surface can be neglected even over distances 
of 100 km.’ 

However, these measurements give no information about any time- 
independent distortion of the beam caused, for example, by imperfect 
lens surfaces. Abberations are most likely negligible for the low power 
lenses used in the lens guide, but random irregularities may lead to 
severe damage of the beam profile in a long sequence of lenses.* ° 
This will impede detection, increase the diffraction loss, and hamper 
any active position control of the beam. 

The difficulty of measuring these distortions lies in the sensitivity 
required. If 10 per cent deviation from the correct intensity profile 
of the beam is tolerable after 1000 lenses a check in one section has 
to be sensitive to 0.01 per cent deviation. In addition, this has to be 
done for many different lenses to gain a reasonable average. The shut- 
tle pulse technique is one method for measuring these small deviations. 


II. THE SHUTTLE PULSE LINE 


A straight 3-14-inch-diameter iron pipe 5 feet below ground with 
lenses roughly every 140 m was used as a test line. The lenses had a 
focal length of 70 m, a diameter of 60 mm, a quality of 4/10, and 
antireflection coatings on both sides. There were six 140 m sections 
with a total length of about 14 mile. For further information about 
this line see Reference 5. 

After a passage of 6 lenses no distortion of the beam could be 
measured reliably. To increase the minute effects involved the beam 
was folded back on itself, passing the same lenses many times. To 
distinguish between the various transits a light pulse shorter than 
the round trip travel time had to be launched. This “pulse shuttling 
technique” was used before? and was adopted here with some varia- 
tions which are shown in Figure 1. 

A polarization switch consisting of a KD*P crystal and a calcite 
Wollaston prism was positioned inside the cavity of a He-Ne laser 
oscillating at 6328A. Activated by a voltage pulse this switch 
“dumped” the energy stored in the oscillator cavity into the test line. 
This light pulse passed through a He-Ne amplifier and a mode match- 
ing two-mirror telescope before entering the lens guide. When the pulse 
returned from its first round trip through the line it found no voltage 
on the switch, passed it once, and another time upon reflection from 
the mirror at the back of the switch, and entered the line once again. 

On every reflection from one of the telescope mirrors about one 
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per cent of the light pulse was transmitted and fed into a receiver. 
This loss and the loss in the switch were compensated by the He-Ne 
amplifier.* By counting the number of pulses in the receiver the 
number of round trips in the line was determined. The total loss in 
the line was calculated from the decay of the pulse train. 

The use of a polarization switch to generate high power pulses from 
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Fig. 1—Shuttle pulse line showing pulse injection, receiver, and processing 
apparatus. 


a He-Ne laser has been described before.1° To reduce the losses in 
the switch both the KD*P crystal and the calcite prism were im- 
mersed in a matching oil chamber with antireflection coated quartz 
windows. Figure 2 shows the refractive indices for the polarization of 
the oscillator light. The z-cut KD*P crystal was 3X3xX12 mm, the 
longitudinal field was applied by silver-painted rims around both 


* Its use was suggested by Mr. E. A. J. Marcatili. 
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ends. A loss of one or two per cent per transit may be expected for 
this switch. 

The 2 m long oscillator cavity was almost confocal. An iris inside 
the cavity limited oscillations to the fundamental mode which had 
a 1/e diameter of 1.25 mm at the mirrors. The gas tube was 1 m 
long with a bore of 4 mm. 

The 2 kV pulse from a commercial unit* had a rise time of 50 ns. 
During this time about half the oscillator light was rotated into the 
perpendicular plane of polarization and deflected into the shuttle 
pulse line by the Wollaston prism. The cavity was emptied com- 
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Fig. 2— Optical switch with numbers giving the refractive indices for light 
polarized in the plane of the figure. 


pletely after about 150 ns. The deflected pulse is shown in Figure 38. 
The peak power of the pulses entering the line was about 100 mw.?° 

The amplifier tube was 1 m long with a 4 mm bore. The telescope 
magnified the 1.25 mm diameter beam to the 10 mm diameter of 
the lens guide mode. The spacing of the telescope mirrors was 1.25 m, 
large enough to neglect spherical aberration and astigmatism. All 
mirrors used in the experiment had dielectric coatings. 

As a necessary shield against turbulence the beam had to be en- 
closed everywhere. In the laboratory acetate tubes were used. In 
addition a wooden cover provided a shield against dust and acoustic 
vibrations. All parts of the launching system were mounted on a 
solid concrete table. The laboratory was an underground concrete 
room of 5X5 m. Air currents and convection in the underground 
conduit were prevented by an air-tight seal closing the tube where 
it entered the laboratory. One of the antireflection coated lenses al- 
ready described was used for this purpose. 

The mirror at the far end of the line had a 140 m radius of curva- 


* Pulse Engineering Inc., Velonex Div. Model 350. 


PULSE SHUTTLING 771 


TIME IN MICROSECONDS 
ie) 0.2 0.4 0.6 





Fig. 3 — Pulse injected into the lens guide (amplitude not calibrated). 


ture and was adjusted by remote control. It could be mounted in 
place of any of the lenses thus changing the guide length from 140 m 
to 6X140 m. The change in the decay rate of the pulse train was 
observed as the guide length was changed. From this it was concluded 
that each section added about 2 per cent to the round-trip loss of the 
pulses. Consequently there was about one per cent section loss per 
transit of which a large fraction can be attributed to the lenses. 
Notice that these lenses had been in use for about a year in previous 
field experiments where moisture and dust might have affected their 
coatings. 

With the amplifier turned off the round-trip loss in the switch, 
the mirrors and the amplifier tube was about 18 per cent. The switch 
and the mirrors probably account for the largest part of the loss. The 
loss in the switch will be somewhat higher for the shuttle pulse polari- 
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Vig. 4 — Pulse decay in the 3-section line with the amplifier turned off. 


zation than for the oscillator polarization. For the perpendicular 
polarization of the shuttle pulse the indices in the calcite are inter- 
changed from those shown in Figure 2, producing a larger mismatch. 
Figure 4 shows a pulse train received from a line with three sections 
corresponding to a round-trip loss of 24 per cent. 

With the amplifier turned on, the 8-section line went into oscilla- 
tion with modes up to the third order being observed. Though not 
as bright as in a cavity with one or two sections these modes were 
stable in time and showed no noticeable distortion of their intensity 
profile. Oscillating modes were a good indication for the optimum 
alignment of all components. In a line with more than three sections 
the gain was not sufficient to cause oscillations. Figure 5 shows a 
pulse train in a 38-section line after the amplifier gain was reduced 
enough to prevent oscillation. 

The detection of the pulses was not limited by detector noise but 
by cw-light leaking from the switch into the cavity during the in- 
terval between pulses. If this light is at a resonant frequency of the 
shuttle pulse cavity a relatively large amount of power can be 
coupled from the oscillator and large field strengths will build up in 
the line. Fortunately this long cavity changes its length randomly 
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during the round-trip time thus avoiding resonance. However, in this 
random case some power will still be coupled into the line and the 
oscillations will build up to an extent determined by the line losses. 
Therefore the amplifier gain had to be set to balance the background 
light and pulse decay so that a maximum number of pulses could be 
seen. This condition is studied quantitatively in the appendix. 

Figure 6 shows the optimum round-trip loss, resulting in the maxi- 
mum number of round trips N, for a given extinction ratio of the 
switch and a desired signal to noise ratio at the detector. The lines 
show the idealized situations of perfect resonance and completely 
random phase buildup. The experimental condition was somewhere 
in between, but close to complete randomness. To get the 150th pulse 
with an S/N of 20 dB, the amplifier gain was set for a round-trip loss 
of about one per cent. This suggests an extinction ratio of about 
50 dB (see Figure 6). This extinction ratio has to be understood as 
the ratio of the peak pulse power to the power that leaks in between 
pulses and propagates in the line. Only the small part of the light 
leakage that is in the lower order modes can propagate. 

The initial high amplitude pulses in Figure 5 suffered a larger 
attenuation than one per cent as they drove the amplifier into satura- 
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Fig. 5 — Pulse train from the 3-section line with the amplifier working. 
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Fig. 6— Round-trip loss versus extinction ratio of the pulse injector for the 
maximum numbcr of pulses detectable. The lines indicate the limits of perfect 
resonance and of complete phase randomness. 


tion. Some of the background light in Figure 5 is caused by acoustic 
resonances in the KD*P crystal following the high voltage pulse. 
Since this did not obstruct the detection of the initial pulses and 
vanished rapidly no measure was taken to prevent it. 


IH. THE SCANNING OF THE BEAM PROFILE 


The light transmitted by the larger telescope mirror in Figure 1 
was used to measure the lateral intensity profile of the light pulses 
impinging on this mirror. The technique was to measure the intensity 
at a certain point of the mirror with a pinhole receiver when, say, 
the 50th pulse comes by, then move to another point, launch a new 
pulse, and wait again for the 50th round trip. This can be repeated 
point by point for any detectable round trip. Since in the experiment 
the round-trip time was several microseconds a new pulse could be 
launched every millisecond. 

The pinhole receiver consisted of a photomultiplier with a small 
iris and a diffusing ground glass plate behind it. The hole diameter 
was a fraction of a millimeter. The light pulses had a spot diameter 
of 10 mm at that mirror. 

The multiplier received all pulses arriving from the far end of the 
line. The transmitted intensity from pulses traveling in the opposite 
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direction was deflected to the side and did not reach the pinhole. A 
diode gate, triggered after the proper delay time, picked the wanted 
pulse out of the amplified multiplier signal. The signal from the gate 
was amplified by a narrow band resonance amplifier which was 
tuned to the 1-kHz repetition rate of the injected pulse. The rectified 
amplitude of the 1-kHz signal was displayed on the y-axis of an x-y 
recorder. The x-deflection corresponded to the location of the pinhole 
receiver which was moved slowly across the back of the mirror. One 
sean took several minutes, the scanning speed being determined by 
the bandwidth of the rectified signal. This bandwidth was reduced 
to a fraction of 1 Hz to suppress time dependent effects as far-as 
possible. The accuracy of the scan was set by variations affecting 
the tremendously sensitive alignment of the system. This limited the 
reproducibility to about 5 per cent of the total signal. 

Figures 7, 8, and 9 show scans taken in a line with three sections. 
The intensity units in these figures are arbitrary and different for 
each scan as the amplifier settings were changed between scans. All 
these scans were taken along a horizontal line across the center of 
the light beam. The dotted lines show a gaussian intensity distribu- 
tion calculated for the theoretical 1/e-spot diameter of 10 mm for 
comparison. 

Figure 7 represents a check of the profile launched into the line. 
The bumps and spikes on this scan are attributed to variations with 
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Fig. 7— Intensity scan across the guide cross-section for light after one round 
trip in the 3-section line, 
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time and the injected beam may be considered gaussian within the 
accuracy of the scan. Figures 8 and 9 are scans taken after 75 and 
150 round trips, respectively. The 3-section line was about 400 m 
long. Consequently, after 150 round trips the light pulse had traveled 
120 km and passed 900 lenses. 

The slight tilt of the profile may be interpreted as light which 
accompanied the beam in higher order modes.2 Many more scans 
were taken in this and longer lines for various numbers of round 
trips. Figures 7, 8, and 9 are representative of the results. 


IV. DISCUSSION OF THE EXPERIMENT 


A question that has to be investigated in more detail is whether 
a 120 km lens guide would actually cause the same negligible dis- 
tortion on a light beam as found for the shuttle pulse experiment. 
There are three important points to be considered: 


(1) The distortion added to the beam by the optical components 
other than the lenses. 

(it) The finite apertures of these components. 

(wu) The periodic nature of a shuttle pulse experiment. 


In the shuttle pulse line there were, in addition to the lenses under 
investigation, several other optical components, such as the switch, 
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Fig. 8— Intensity scan across the guide cross-section for light after 75 round 
trips in the 3-section line. 
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Fig. 9 — Intensity scan across the guide cross-section for light after 150 round 
trips in the 3-section guide. 


the amplifier, and the telescope, which the light passed twice on 
every round trip. To find their contribution to the beam distortion, 
the nature of the expected irregularities on an optical surface has 
to be considered. 

Any irregularity on an optical surface will distort the beam, but 
not all of this distortion will accumulate in a lens guide, because 
the finite aperture of the guide smooths out the intensity and phase- 
front of the beam. This is equivalent to saying that only low order 
modes will be transmitted in the guide. Consequently these modes 
are the only growing source of distortion. They are generated by ir- 
regularities correlated over areas comparable to the beam cross 
section: smooth bumps whose heights determine the amount of light 
coupled into these modes.® 

The polishing process will normally cause those bumps to be the 
highest which extend furthest over the surface. Hence it is expected 
that a large beam will be distorted more than a small one. This 
implies that those optical components on which the beam had only 
a 1 mm diameter contributed very little to the distortion. So only 
the two large mirrors facing the line are left as major additional 
sources of beam distortion. Since they had the same optical quality 
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and beam size as the lenses, their contribution may be considered 
approximately the same as an additional lens in the line. 

The second point deals with the effective beam aperture of the 
shuttle pulse line. The amplifier had an inner diameter of 4 mm, and 
the KD*P crystal had a cross section of 3X3 mm. Both were in the 
2 m section between the small telescope mirror and the end mirror. 
For a wavelength of 6328A, an aperture of 3 mm, and a spacing of 
2 m, one calculates an effective Fresnel number of 1.8. In this sec- 
tion the beam had the same shape and phase front as in the confocal 
2 m oscillator cavity from which it was originally ejected. 

Applying Slepian’s propagation coefficients for modes in a diffrac- 
tion-limited confocal cavity one finds, for a Fresnel number of 1.8, 
a diffraction loss of about 10-°, 10°, 10-*, 10-° per transit for the 
first, second, third, and fourth parasitic mode, respectively.7?_ More 
significantly, all these modes except the fourth one went into oscil- 
lation in the 3-section line. Hence the aperturing in the system should 
not suppress the mode buildup for modes up to at least the third 
order. 

Further investigation of the individual modes shows that the first 
and second parasitic modes are not only generated by irregularities 
on the lens surfaces, but also by additional lens effects, for example, in 
the gain medium of the amplifier tube or by lens displacements, later- 
ally and along the axis. In turn these modes can be suppressed, even 
in a line with imperfect lenses, by proper alignment of these lenses.® 
This means that the first and second parasitic modes cannot be used 
to determine the lens imperfections in the experimental guide. They 
could be and were avoided by proper alignment. 

Let us look therefore at higher modes, say, the third one. Figure 
9 in Reference 9 shows the profile of a gaussian beam with 5 per 
cent of its power in the third mode. Comparing this with the recorded 
profile in Figure 8 one may assume that its distortion was of the 
same order. Without pushing the quantitative analysis of the experi- 
ment too much, it can be said that the power collected in the third 
mode after 900 lenses was apparently not larger than a few per cent 
of the beam power. 

Gloge® predicts this amount (and only a fraction of this in the next 
higher modes) for the lenses used, if the irregularities on the surfaces 
are correlated over distances about twice the beam diameter. On the 
basis of these calculations it may be assumed that the detection of 
the higher order modes would be beyond the experimental accuracy. 

The third point has to do with the iterative structure of the lens 
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TIME IN MILLISECONDS 





Fig. 10— Pulse train from the 3-section line after one lens was misaligned. 


guide simulated by the shuttle pulse line. Obviously on every round 
trip the light pulse suffers the same kind of distortion, whereas in a 
real lens guide this distortion would be completely at random. 

Consider for example the most simple irregularity: a displaced 
lens. Figure 10 shows a pulse train similar to Figure 5 but after one 
lens was misaligned. Apart from the expected decay there is a pe- 
riodic fluctuation of the signal. This can be explained by an oscilla- 
tion of the beam center about its correct position or, which is equiva- 
lent, by periodically generated and annihilated parasitic modes. One 
can show that the modes, because they travel with a propagation 
constant slightly different from that of the gaussian beam, should 
do exactly that in a periodic structure. . 

Only if the round trip is an integral number of beat periods for a 
specific mode will the power in that mode grow continuously. To 
multiply the minute effects expected per section in the experiment 
this condition was arranged for the third order mode. This was done 
by changing the distance of the telescope mirrors and watching the 
periodic fluctuations of the pulse train. One can calculate that the 
third mode builds up if the fluctuation period is three pulses on the 
pulse train. In this case we should observe the maximum third order 
mode distortion of the beam. Scans taken under this condition indi- 
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cated no major deviation from the scans shown. The explanation 
might be that the power generated per round trip in the third mode 
was smaller than, or at most, equal to the attenuation of that mode. 

In a lens sequence with independent random irregularities, the 
relative phases are completely arbitrary and the average power in 
the parasitic modes would grow proportional to the number of lenses 
passed, at least as long as reconversion can be neglected. If, as an 
upper limit, the generation rate per round trip is assumed equal to 
the loss of 2X10-‘, the third mode would increase to 3 per cent of 
the beam intensity after 900 equivalent lenses. 


V. CONCLUSIONS 


A 1%-mile underground lens guide with six confocally spaced lenses 
was investigated by sending a light pulse back and forth in this 
guide. The lateral intensity profile of the light beam could be checked 
after any round trip by scanning the intensity in the cross section. 
It was found that the gaussian profile after passing 900 lenses was 
still virtually undistorted. 

The lenses had a quality of 4/10, a diameter of 60 mm, and were 
antireflection coated to give a reflection and absorption loss of one 
per cent per pass. The effective aperture of the guide was determined 
by an amplifier tube and an optical switch within the path of the 
light. The diffraction losses for modes up to at least third order were 
negligibly low. Even for the condition of critical buildup of the third 
mode the distortion of the light profile after 900 lenses was negligible. 
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APPENDIX 


Assume that the peak pulse power entering through the switch is 
I, and that the pulse loses a fraction p of its intensity on every 
round trip. Consequently, after N round trips, the intensity 


Le Ay Ee (1) 


is left in the pulse. At the same time a power J, leaks in through the 
switch continuously and increases the background light J, to a level, 
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where J, just replenishes the loss per round trip. If J; adds in phase 
to the light already circulating, as in a perfect resonator, the com- 
parison of amplitudes yields 

Vi —p)+ VL = Vi. (2) 
If the phases are completely at random because either the resonant 


frequency or the entering light frequency changes within the round 
trip time one has instead 


Ld-p+li,=f,. (3) 


By introducing a stability parameter a = 1 or 2, both (2) and (38) 
can be expressed by the relation 


Dyed ep) (+) 
which for p < 1 reduces to 


1,(2) 
a 


It is convenient to define a signal to noise ratio 


k= 10 log = 
b 


required at the receiver to analyze the Nth pulse and an extinction 
ratio 


KE = 10 log = 


inherent in the switch. By combining (1) and (5) and solving for 
N, one now finds the number of round trips which can be properly 
detected: 


R — E — 10a(log p — log a) 








oe 10 log (1 — p) (6) 
The derivative dN /dp vanishes for 
N = qi— Bun. (7) 


Popt 
This is the maximum number of round trips that can be achieved 
with the optimum round-trip loss, for example by adjusting the 
amplifier gain. Equations (6) and (7) were used to plot Figure 6. 
It shows pope as a function of the difference H — R for a = 1 and 2. 
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These two lines present the limits of perfect resonance and completely 
random phase buildup. For example, to receive the 150th pulse 20 dB 
above the noise level in a line that does not resonate, one needs an 
extinction ratio of 46 dB in the switch and a loss of 0.67 per cent in 
the shuttle pulse line. 
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Linear Beam Position Control 


in Optical Waveguides 


By JAMES C. DALY 
(Manuscript received September 21, 1967) 


This paper studies the effect of an active control system that introduces 
corrective transverse lens displacements on the propagation of beams in 
optical waveguides. The induced displacement of the nth lens is a linear 
function of the beam displacement at the (n + 1)th lens. Beam displace- 
ments from the ideal design position are caused by transverse displacements 
of the lenses and sensors from their design positions. Expressions have 
been obtained for: 


(t) The conditions for spatial and temporal stability. 
(22) The rms beam displacement resulting from uncorrelated displace- 
ments of the lenses and sensors. 
(iz) The beam displacement caused by sinusoidal displacements of the 
lenses and sensors of arbitrary spatial frequency. 
(iv) The spatial rate of return to the guide axis of a beam injected into 
the guide off axis. 


We also show that the positions of beams, other than the sensed beam, 
are controlled by the system and that the system can simultaneously guide 
many different sensed beams. We calculate an upper limit on the prob- 
ability that the beam will not be contained within a given aperture. 


I. INTRODUCTION 


The ability of a beam waveguide to guide optical frequency elec- 
tromagnetic energy is severely limited by the tolerances held on the 
transverse positions of the lenses. Hirano and Fukatsu* have shown 
that uncorrelated transverse lens displacements will cause the ex- 
pected deviation of the beam from the guide axis to grow as the 
square root of the number of lenses through which it has passed. 
Steier? has shown that random longitudinal lens displacements and 
variations in the focal lengths of the lenses couple with random 
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transverse lens displacements to cause the expected value of the rms 
beam displacement to grow exponentially with the number of lenses 
through which it has passed. 

This paper describes the performance of linear self-aligning beam 
waveguides that are subjected to transverse disturbances of the lens 
and sensor positions. It is quite likely that nonlinear control systems 
may turn out to be more effective or practical in practice but an 
understanding of linear systems will facilitate the more involved 
analysis of nonlinear control systems. The idea of using active guid- 
ing media for optical communications was proposed long ago by 
R. Kompfner? and L. U. Kibler,t and has been discussed by E. A. J. 
Marcatili.» The improvement that can be obtained with self-aligning 
beam waveguides has been demonstrated experimentally and with 
computer simulations by Christian, Goubau, and Mink.® 

In the guiding systems considered in this paper, the lenses are 
physically moved an amount that is determined by sensing the posi- 
tion of the beam in the guide. 

The axis of the guide is defined as the line joining the centers of 
the beam sensors. In general in an actual installation this axis will 
not be straight. 

The displacement of paraxial beams from the guide axis is a linear 
function of the transverse lens displacements. An important conse- 
quence of this linearity is that the superposition principle can be 
applied to the beam displacement. Incremental displacements of the 
lenses produce incremental displacements of the beam. These incre- 
ments are independent of any bias positions resulting from-an inten- 
tional bend. Therefore, when we study the effects of imperfections in 
a straight guide, the results apply equally well to the effect of im- 
perfections in a curved guide. For example, in studying the effect of 
the control mechanism on the propagation of a beam through a bend 
we would proceed as follows: 


(1) Specify the guide axis. This is the path of propagation for the 
beam. The sensors are located along this axis. 

(iz) Find the proper set of lens positions that will cause the beam 
to propagate along this axis. In an analysis these are calculated. 
In an actual installation these positions would be determined by 
moving the lenses until the beam is aligned along the sensors. 

(121) Determine the effect of imperfections by assuming that the 
lenses and sensors are displaced from their proper positions. The 
consequence of the linearity of the system is that the beam displace- 
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ment resulting from these imperfections is the same as it would be if 
the guide axis were a straight line. 


This paper emphasizes the case where the design axis for the lens 
centers and sensors are straight lines and coincident with a reference 
line. The beam is injected into the guide along the reference line and 
the center of the first lens of the guide lies on this line. 

The optical waveguide consists of a sequence of identical positive 
lenses, of focal length, f, and spacing, d, as shown in Figure 1. In 
Figure 1, gn, s,, and r, are the sensor, lens, and beam displacements 
at the n™ lens, respectively. We shall consider paraxial ray propaga- 
tion in two dimensions for simplicity. It has been shown that the three- 
dimensional problem can be split into two independent two-dimen- 
sional problems? and that the center of the beam follows a path 
described by a paraxial ray.” 8 

When the guide is subjected to disturbances in time, the transverse 
displacement of the n lens relative to the reference line is given 
by s(t),. Similarly, the displacement of the beam center at the n™ 
lens from the reference line is given by r(t),. Using geometric optics,* 
and considering the guide to be time-dependent, results in the follow- 
ing difference equation for the paraxial beam in the guide 


PO)nv2 — (2 COS wd )r)nv1 + rOn = /f)sOn+i (1) 


where cos w,d = 1 — d/2f. 

In any lens waveguide there will be transverse displacements of 
the lenses, owing to the impossibility of holding exact tolerances or 
the movement of the earth after the guide has been constructed. 
These lens displacements are unavoidable. They are denoted by N,. 
Here we study the effect of a superimposed corrective lens displace- 
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Pig. 1— The beam waveguide and schematic diagram of the feedback system. 
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ment, 8,, on the performance of the beam waveguide. The total dis- 
placement of the nt lens is the sum of N, and 8, 


Sn = N, + bn rd (2) 


For most of this paper, the response of the system in time is deter- 
mined only in the steady-state when the disturbances are either step 
functions in time or at low enough temporal frequencies to be con- 
sidered constant functions independent of time. These conditions are 
relaxed and general time functions allowed in Section 2.2 where tem- 
poral stability is considered and in Section 2.4 where the high gain 
case is considered. 

Although the beam control is achieved here by inducing transverse 
lens displacements, the results can be applied to other beam control 
mechanisms such as mirrors and prisms, when such mechanisms can 
be interpreted as an equivalent transverse displacement of the lens. 


Ii. THE SYSTEM 


Figure 2 is a schematic representation of a mechanism for lens 
position control. 

The beam position is sensed at each lens. This beam position signal 
is used to control the transverse position of the preceding lens. The 
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Fig. 2— Schematic diagram of a mechanical device for lens positioning. 
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difference between the beam and sensor positions at the nth lens we 
call the error in the beam position. The corrective displacement of 
the nth lens, 5,(¢), is a function of time and the error in the beam 
position at the (n+1)th lens. 


6, = FCn+1 ? t) (3) 


where €,+1 is the position of the beam relative to the sensor at the 
(n+1)th lens. The form that the function F(e,41,t) takes will depend 
on the properties of the control mechanism. 

In order to study the stability of the beam in the guide it is con- 
venient to transform equation (8) into the temporal frequency domain. 
Assuming that the control mechanism can be represented by a transfer 
function H,(s), the Laplace transform of the corrective displacement 
of the nth lens for a system at rest at time zero can be written: 


da(s) = A.(s)E n+, (s) (4) 
where s, 5,(s), and #,41(s) are the complex temporal frequency, the 
Laplace transform of the corrective displacement of the n™ lens, and 
the Laplace transform of the error in the beam position at the (n+1)® 
lens. 

In general the transfer function will have the following form 


FF eee nL a 
s(sr, + 1)(s72 + 17,3 + 1) -:: 


It is customary to classify systems in terms of the number, JN, of 
factors of the form 1/s which appear in the transfer function. This 
number indicates the number of integrations performed by the com- 
ponents. For example: if VN = 0 the system is referred to as a type 0 
system, if N = 1 the system is referred to as a type 1 system, and so on. 

The number of integrations determines some important properties 
of the steady state response. Since s ~ 0 as t > oo, in the steady state 
equation (5) yields 


(5) 


_* 


H,(s) |, = A= V (6) 





iv) 


when s — 0. There are two possibilities for A. A is either equal to a 
real constant K (type 0 system) or A goes to infinity (type one or 
more system). We consider both of these possibilities in the following 
paragraphs., 

For the control mechanism shown in Figure 2, H,(s), the control 
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mechanism transfer function, will be of the form 


K, 
Ms’ + Bs+ K, (7) 
where K, is the gain, K, is the spring constant, M is the mass of the 
lens, and B is a friction constant. The steady-state value of H,(s) = 
A = K,/K,. This system is reduced to a type 1 system by removing 
the spring, that is, letting K, = 0. 


Hs) = 


2.1 The General Description 


The properties of the linear active guiding medium in time and 
space can be studied using a combination of Laplace and Z transform 
techniques.°"* Laplace transforms are used to transform from the 
time domain to the temporal frequency domain. Quantities that de- 
pend on position along the guide, given by the lens number 2, are 
transformed into quantities that depend on axial spatial frequencies 
using Z transform techniques. 

Taking the Laplace transform of equation (1) results in the follow- 
ing recursion equation for the Laplace transform of the beam dis- 
placement. 


[1 — (d/f)H.(s)|RG)n+2 — (2 cos w,d)R(S)n+1 


Gase 
§ 


+ RG), = (a/f) *** — (4/NH) (8) 


where R(s)n, Nn, and G, are the Laplace transform of the beam 
displacement at the nth lens, and the constant displacement of the nth 
sensor, respectively. Taking the Z transform of equation (8) results 
in the following equation when the boundary conditions are zero, 
that is, the beam is injected into the guide on axis and the first lens 
of the guide is on axis 


Ris, 2){{1 — (d/f)H.(s)]e” — (2 cosw,d)z + 1} 
((z) 


= d/) *2 — 2a/puy"2  @ 


where R(s,z) is the Z transform of the Laplace transform of the beam 
position, and N(z) and G(z) are the Z transforms of the constant 
lens and sensor displacements, respectively. Rewriting equation (9) 
and introducing the transfer function H(s,z) results in the following 
expression. 


R(s, 2) = H(s, 2)[N(z) — 2H.(s)G(z)]/s (10) 
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where 





2(d/f) 

Me) Gay © esas te 1 (11) 
Equation (10) is a simple relationship between functions that contain 
the information about the performance of the guide. The transfer 
functions H(s,z) and H,(s) contain the properties of the guide such 
as d/f and the characteristics of the control mechanism. The function 
R(s,z) contains all the information that is necessary in principle to 
determine the position of the beam at any lens as a function of time. 
N(z) and G(z) are determined from the constant lens and sensor 
displacements. 

When r,(t) can be considered a constant function of time denoted 
by 7,, then R(s,z) is the Z transform of r,/s. Multiplying equation 
(10) by s and letting the temporal frequency go to zero yields R(z), 
the Z transform of the steady-state beam position 


Riz) = H(@)[N(z) — AzG@)] (12) 


where H(z) is equal to the value of H(s,z) for s equal to zero. 


2.2 Temporal Stability 


A function is defined as stable in time if it is bounded for all posi- 
tive values of time. The position of the beam in the guide will be 
stable in time if the poles of the Laplace transform of the beam posi- 
tion are in the left half of the complex temporal frequency plane (s 
plane). The Laplace transform of the beam position is obtained by 
taking the inverse Z transform of equation (10). The following con- 
volution summation’ is a general expression for the inverse Z trans- 
form and gives the Laplace transform of the beam position. 


RAS) = 5 2 Ha)Noom — He)G rma] (13) 


where H,,(s) is the inverse Z transform of H(s,z). Assuming that the 
transfer function H,(s) (that describes the control mechanism) has 
no poles in the right half of the s plane, the only source of right half 
plane poles would be in the function H,,(s). In order to obtain an 
expression for H,,(s) it is convenient to factor the denominator of 
equation (11) and rewrite the expression for H(s,z). 2; and zg are the 
two roots of the denominator in equation (11) and are given by the 
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following expression 


cos w.d + [cos’ wd — 1+ (d/fH.(s)]'. 











®1and2 = 1 — (d/pH.As) (14) 
Therefore 

_ a/f 2 : 

Sa F = fra ee Ere ae) 


The inverse Z transform of equation (15) can be found in tables and 
is given by 


___ @/2fet =) 
Hal’) = Te? wd = 1-2 @/PH.ol 


A sufficient condition for temporal stability is that the poles H,,(s) 
are in the left half of the temporal frequency plane. The s plane poles 
of H,,(s) are the roots of the following two equations 


1 — (/f)H.(s) = 0 (17) 





(16) 


and 


sin’ wd — (d/f)H,(s) = 0. (18) 


Equation (17) is the characteristic equation of a single sensor lens 
loop. Therefore if a single section of guide is stable the roots of equa- 
tion (17) will be in the left half of the s plane. Equation (18) is more 
restrictive. For example, if a single section becomes unstable when 
the gain, K in equation (5), is increased beyond a value K,, it follows 
from equation (18) that when many sections are tied together to form 
the optical waveguide, the value of K at which the system becomes 
unstable is reduced by sin? ,d. 

When the H,(s) is given by equation (7), the system is stable if 


A<1- 4d/4f (19) 
where A = K,/Ksj. 


2.3 Spatial Stability 


A guide will be defined as spatially stable if the beam remains 
within a bounded region around the guide axis for all values of n. 
A sufficient condition for the steady-state time solution to be spatially 
stable is that the poles of H(z) are inside the unit circle in the Z 
plane. H(z) can be obtained by setting s = 0 in equation (11). 
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2(d/f) 

Nis [1 — (d/f Ale — (2 cosw,d)z + 1 om) 
Beginning with equations (19) and (20), using standard techniques 
it can be shown that for the type 0 system described by equation (7) 
only a negative value of A will satisfy the conditions for both tem- 
poral and spatial stability. For type 1 or higher systems A > o and 
the poles of H(z) approach the origin of the Z plane assuring spatial 
stability. 


2.4 The High Gain Case 


By high gain we mean large magnitude of H,(s). This is of par- 
ticular importance, since the performance of the system improves 
when the magnitude of H,(s) becomes large. When the magnitude of 
H,(s) is large, equation (11) becomes 


peas 
2H_(s) 
Using the above approximation for H(s,z) and relaxing the restriction 
that the disturbances are constant functions of time and substituting 
into equation (10) the general expressions N(s,z) and G(s,z) for the 
Z transforms of the Laplace transforms of the lens and sensor dis- 
placements yields 


H(s, 2) & (21) 


R(s, 2) — G(s, 2) & Sea (22) 


Division of a Z transform by 2 corresponds to subtracting one from 
the index of the inverse transform. The inverse Z transform of equa- 
tion (22) is 


Nn-1(8) 
H.(s) 
when H,(s) is large the Laplace transforms of the beam and sensor 
positions will be approximately equal. It follows that the functions in 
the time domain will also be approximately equal. In the steady-state, 
when the disturbances are independent of time, the displacement of 
the beam from the guide axis is 


R,(s) — G,(s) & — (23) 


ee ie A, (24) 


For a type 1 or more system A — o and in the steady-state the beam 
follows the sensors. 


792 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


2.5 The Propagation of Unsensed Beams in the Guide 


An unsensed beam will differ from the sensed beam primarily be- 
cause the active guiding mechanism will not respond to deviations of 
the unsensed beam from the guide axis. The unsensed beam may also 
be injected into the guide off axis and follow an entirely different 
trajectory in the guide than does the sensed beam. | 

The effect of the corrective displacements on beams other than the 
beam that is sensed by the control system, can be scen by noticing 
the form of the total lens displacement. The Z transform of the total 
lens displacement, S(z), is the Z transform of the constant lens dis- 
placement, N(z) plus the Z transform of the corrective lens displace- 
ment 8(z2) = Az[R(z) — G(z)], 


S(z) = N(z) + Az[R(z) — Ge)]. (25) 
Substituting equation (12) into equation (25), the Z transform of the 


total lens displacement becomes 


2” — (2 cosa,djze + 1 


Be) (1 — Ad/f)z” — (2 cosw,d)z + 1 





[N(z) — AzG(z)}. (26) 


Let P,, be the displacement of an unsensed beam propagating in the 
guide with the boundary conditions P;, and P%, just to left of the 
first lens which is assumed to be on the reference line. The difference 
equation that describes the propagation of any beam sensed or unsensed 
traveling in either direction through the guide is equation (1) rewritten 
here to describe the unsensed beam. 


Pase — (2 C08 wo) Pass +P, = (A/f) Sasi - (27) 


Operating on equation (27) using Z transform techniques?" yields 
the following expression for the Z transform of the position of the 
unsensed beam 


2(d/f)S(2) + (2? —z2]Pin + 2dP in 
2’ — (2 cosw,d)z + 1 
The substitution of equation (26) and (20) into equation (28) yields 


the following expression for the Z transform of the displacement of the 
unsensed beam. 





P() = (28) 


* = 2)Pin + 2dP iy 

— (2 cosw.djz + 1 
The displacement of the unsensed beam, P,,, the inverse Z transform 

of (29), is the superposition of two components. The first is represented 





P) = H@IN® — AzG(2)] +# (29) 
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by the first two terms in equation (29). A comparison of these two 
terms with equation (12) shows that this component of the displace- 
ment is identical to the displacement of the sensed beam. Comparison 
of the remaining terms of equation (29) with equation (28) shows 
that they are identical to the Z transform of the displacement of a 
beam traveling in a guide with zero lens displacement, S(z) = 0. 

The second component is identical to the displacement of an* un- 
sensed beam injected off axis into a perfectly straight guide. The 
total displacement of the unsensed beam is the sum of these two 
components. The unsensed beam will oscillate about the position of 
the sensed beam just as a beam injected off axis into a perfectly 
straight guide would oscillate about the guide axis. This result differs 
from equation (28) of Reference 5, which predicts a growth in the 
displacement of the unsensed beam as it propagates. In Reference 5 
substitution shows equation (28) to be inconsistent with equation 
(22). 

For stability, the sensed beam must travel in the direction of increas- 
ing n, but an unsensed beam may travel in either direction. 


2.6 The Propagation of Many Sensed Beams in the Guide 


Where there are many beams propagating in the positive direction 
of the guide, the signal obtained from the sensors is the average dis- 
placement of all the beams. The Z transform of the beam position 
signal is 


D(z) = (P@)) — G@). (30) 


Where D(z), (P(z)), and G(z) are the Z transforms of the position 
signal, the average beam displacement, and the sensor displacement, 
respectively. Different beams are injected into the guide differently 
and follow different trajectories through the guide. When the Z trans- 
form of the total lens displacement is S(z), the Z transform of any 
beam propagating in the guide is given by equation (28). The Z 
transform of the average displacement is given by the average of 
equation (28) over all beams. 


pig) = alS@) + be — e\Pa) + 2dlPn) 31) 


2’ — (2 cosw.d)ze + 1 


Where (P;,) and (P%,) are the average displacement and slope of the 
beams just before propagating through the zeroth lens. It follows 
from equation (31) that the average displacement of the beams in 
the guide is the same as the displacement of a beam that is injected 
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into the guide with slope and displacement equal to (P%,) and (Pin), 
respectively. Therefore, the control system will respond to the average 
beam position signal and control the average of the beam displacements 
just as it would respond to one beam propagating in the guide. 


2.7 The Response of a Beam to Waves of Displacements 


As an example of how the system might perform consider the fol- 
lowing particular type of disturbance. Both lenses and sensors are 
displaced simultaneously. The constant transverse displacements, D,, 
of the lens and sensor assembly relative to the straight reference line 
form an axial wave of spatial frequency o, and are described by the 
following equation: 


D, — Bei**™ (32) 


then 


DE) = 





ane (33) 


where B is the amplitude of the wave, o = 27/d, and dX is the spatial 
wavelength. When n is large enough so that boundary effects can be 
neglected, the Z transform of the displacement of the beam from the 
sensor is obtained by using equation (20) and letting D(z) = N(z) = 
G(z) in equation (12). 


=e 
(1 — Ad/f)z” — 2 cosw.dz + 1 





Ri) — Dz) = — De). (34) 
It follows*® from equation (34) that at any given spatial frequency, 
w, the ratio of the amplitude of the wave of the beam’s displacement 
from the lens sensor assembly, to the amplitude of the wave of the 
lens sensor assembly’s displacement from the axis is given by 


m—-~D,)_ | cosed | 
D, | — | coswd — cosw,d — A(d/2fe'** | 





(35) 


The ratio given in equation (35) is plotted in Figure 3, for various 
real values of A, as a function of the spatial phase, od. It follows from 
the sampling theorem that the largest spatial phase possible between 
lenses is od = x. 

Figure 4 is a plot of the ratio of the amplitude of the beam displace- 
ment to the lens displacement when the sensors remain on axis and the 


BEAM POSITION CONTROL 795 


2.5 





hu 
° 





w 





(T-D) max/ 2 max 
S 





0.5 





2 0.4 06 14.0 2.0 4.0 
SPATIAL PHase (wd = 47d) 


Fig. 8 — Spatial frequency response of the beam position relative to the sensor 
when the lenses and sensors are disturbed simultaneously. 


lens displacements NV, form a wave of axial spatial frequency o. For 
this case, the ratio of the amplitude of the beam displacement wave 
to the amplitude of the lens displacement wave is 








Lh ee.) ed (36) 
N,, | coswd — coswd — A(d/2f)e'*" | 

















SPATIAL PHASE (wd= 57d) 


Fig. 4— Spatial frequency response of the beam position relative to the sensor 
when sensors are undisturbed and the lenses are disturbed. 


796 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


2.8 The Response of a Beam to Random Transverse Displacements 


Again assuming that the entire lens sensor assembly moves together, 
the spatial response of the beam to any general set of displacements, 
D,, is given by the inverse Z transform of equation (34). The dis- 
placement of the beam from the sensor is the following convolution 
summation 


ie ~~ D,, = —D, ae > (fee a Aligoncda) Ds, (37) 


ma) 


where h, is the inverse Z transform of the transfer funetion H(z), 
riven in equation (20). 


d/f __ sin nwed 
(1 — Ad/f)""'” sin wed 
where cos wed = (1 — d/2f)/(1 — Ad/f)/*. We wish to find the ex- 
pected value of the square of the beam displacement, from the sensor, 
((tn—D,)?). For uncorrelated displacements, the expected value of the 
product of any two displacements is given by 





hy = (38) 


Ops ‘ Wi a (39) 


o, if n=m 


where oy, is the rms lens and sensor displacement. The expected value 
of the square of the beam position can be written* 


(r. — Dot) = oi EAM Sa. — ann} ao 


This summation has been evaluated for large n and plotted in Figure 
5 for d/f = 1, 2, and 3 as a function of A. When the magnitude of A is 
large compared with one, the ratio of the rms beam to rms lens and 
sensor displacement approaches 1/|A|. 


2.9 The Ability of a Short Section of Stabilized Guide to Return an Initially 
Off-Axis Beam to the Optical Axis 


The preceding sections are concerned with the effect of lens dis- 
placements on a beam that is injected into the guide with zero slope 


* Equation (40) can be obtained using equation (37) with the following con- 
siderations. (2) The displacements, Dn, are uncorrelated. (#%) The expected value 
of a sum is the sum of the expected values. (%#) The expected value of the prod- 
uct of a deterministic function and a random variable is the product of the 
deterministic function and the expected value of the random variable. (iv) The 
zeroth lens and sensor are on axis. 
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Fig. 5 — The statistical response of beam position relative to the sensor when 
lenses and sensors are disturbed simultaneously. 
and coincident with the guide axis. This is equivalent to saying that 
the boundary conditions of the difference equation are zero. 

When a beam is injected into a straight feedback controlled guide 
off axis and the lens displacements are zero the Z transform of the 
beam position is a function of the boundary conditions and is deter- 
mined using standard techniques.?° 


rete [1 — Ad/fjz” — (2 cosw,d)z]7y5 + (1 — Ad/f)r; 
(1 — Ad/f)z” — (2 cosw,d)z + 1 
where 7) and r; are the beam displacement at the zeroth and first 
lenses, respectively. The inverse Z transform of equation (41) can be 
found in tables and is equal to the following 


r, = I — Ad/f)-’” (B sin wend + C cos wend) (42) 


where cos wed = (1—d/2f)/(1—Ad/f)*”? and the constants B and C 
chosen so that equation (42) will satisfy the boundary conditions. 

Equation (42) shows that when a beam is injected into the feed- 
back controlled guide off axis the maximum possible deviation from 
the axis decreases by a factor of (1—Ad/f)-’/? at each lens. For ex- 
ample, a beam injected parallel to the axis, one centimeter off axis, 
into a stabilized guide four lenses long with Ad/f = —9, will leave the 
guide with a maximum displacement from the axis that is less than 
0.1 mm. 


2.10 The Probability That the Beam Will Not be Contained Within a Given 
Aperture 


An upper limit on the probability of the beam exceeding a threshold 
can be calculated using Chebyshev’s inequality. The probability of 





(41) 
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the magnitude of the beam deviation from the axis of the sensors 
exceeding a value, 7’, at any given lens is less than the mean squared 
deviation divided by T? 


Pim —g 127 3% P On) ), (43) 


If the magnitude of “A” is much larger than one, it follows from equa- 
tion (24) that 








2 

(tm, = gn)") = Fs (44) 
where oz is the rms disturbance of the lens position. Equation (44) is 
quite general and holds for deterministic as well as random lens dis- 
placements. 


III. CONCLUSIONS 


A beam waveguide using the linear beam control system described 
here can be made stable in time and space by proper choice of the 
control mechanism. 

If many beams are propagating in the guide it is only necessary 
that the control system interact with one or more beams propagating 
in the positive direction for the guide to be stabilized for all beams. 

When the magnitude of A, the low frequency value of H,(s) is large, 
the error in the beam position varies inversely with A. 

It follows from equation (24) that, if A is large, the probability that 
the deviation of the beam from the axis of the sensors will exceed a 
value, 7, at any given lens is equal to the probability that the dis- 
placement of the preceding lens will exceed a value of 7 multiplied 
by A. If A = —100, and T = 1, and if the lens displacements have 
a standard deviation of 0.1 and a gaussian amplitude distribution, then 
the probability that the deviation of the beam from the guide axis 
will exceed one is less than the vanishingly small value of 107'°. A 
more conservative number, 107°, for this probability, that holds for 
all amplitude distributions, is obtained using Chebyshev’s inequality. 

There is no growth in the rms deviation of the beam from the guide 
axis as the beam propagates down the stabilized guide. Also, since 
the effects of guide imperfections, for example, transverse lens dis- 
placements, do not build up, the construction tolerances on the guide 
can be relaxed. It is possible to contain the beam within a narrow 
region of the guide axis, which permits minimizing the guide aperture. 
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Communication of Analog Data from a 
Gaussian Source Over a Noisy Channel 


By A. D. WYNER 
(Manuscript received December 20, 1967) 


We consider the problem of the transmission of analog data from a 
Gaussian source over a memoryless channel with capacity C nats per 
second. The source emits R independent zero mean Gaussian variates 
per second with variance o°. These digits are block-coded RN at a time 
into N second channel inputs. The performance criterion 1s the mean 
square error. Let &(N) be the smallest attainable mean square error 
with parameter N(R, C, o° fixed). Shannon has shown that &(N) = 
o exp (—2C/R) = & and ¢(N) > & as N > «. Hence the ideal error 
é, ts attainable in the limit as the coding delay N — ~. We are concerned 
with the rate at which &(N) — & , and our principal result is that &(N) — 
& < Of(log N/N)’']. 


I. INTRODUCTION 


We are interested in the following problem. Suppose we have an 
analog data source which emits a sequence of statistically independent 
Gaussian variates at a rate of R per second. We wish to transmit this 
data through a noisy channel of capacity C nats per second. Our prob- 
lem is the determination of the minimum possible mean-squared-error. 

Specifically we shall study the communication system of Figure 1. 
The output of the analog source is a sequence X,, X., X3 --: of 
statistically independent Gaussian variates with zero mean and variance 
o which appear at the coder input at rate of R per second. After N 
seconds, n = NRF source variates have accumulated at the coder input. 
Let X denote this random n-vector. The channel is a discrete memoryless 
channel* which we assume accepts one input per second, and the coder 
contains a mapping of X to an allowable channel input N-vector S. 
Since it requires N seconds to transmit S, the system can process the 
data continuously without a “backup” at the coder input. 


* Actually our results are valid for a broader class of channels. See the remark 
after Theorem 2 in Section II. 
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The decoder examines the channel output N-vector R and emits a 
Euclidean n-vector X’ which is hopefully ‘‘close”’ to X. The error 
criterion which we adopt is (the ‘‘mean-squared-error’’) 


2 


€ — 


Aes tae fees 
~ Ht \|X — xX’ If (1) 
where “‘|| ||” is the Euclidean norm and £ denotes expectation. 

We shall assume that the parameters o*, R, and C are held fixed 
for this entire paper, and denote by ¢(N), the smallest value of ¢ 
attainable with parameter N (and therefore n = RN). Shannon’ has 
shown that 


(i) &(N) = o° exp (—2C/R) = &, (2a) 
(i) &@(N) —-e a Noo, (2b) 


/ 
ANALOG X Ss R Ss 
SOURCE CODER CHANNEL DECODER 


Fig. 1— Communication System. 


so that & = 6 is attainable in the limit as the delay N > o. We are 
concerned here with the rate at which «(N) approaches the ideal é , 
and our principal result is 


2(N) < ali +g a + o(S5) |. i. Nee GN 


where 6 > 0, a parameter related to the channel, is defined in equation 
(13). 

This result is related to a result of D. Sakrison® which was done 
independently.* In fact we have used one of his ideas (Lemma 1 in 
our paper) to simplify our original proof. 





II. STATEMENT AND DISCUSSION OF RESULTS 


Following Shannon’s technique,’ we separate the coder into two parts 
as shown in Figure 2. The first part, called the source encoder or quantizer, 
contains a fixed set $8 of M Euclidean n-vectors, and associates with 
each possible input n-vector X a member of 8 (say X). Let us denote 


* This paper and Sakrison’s paper were presented at the International Sym- 
posium on Information Theory, San Remo, Italy, September 1967. 
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nN 
Xx SOURCE Xx CHANNEL Ss 
ENCODER ENCODER 


Tig. 2— Decomposition of the Coder. 


the resulting (mean-square) ‘quantization error” by 
1 5 
= 7 \|X— Xf. (4) 


The second part of the encoder, called the channel encoder, associates 
with each X (one of the JJ members of $) a channel input N-vector 
(say S). The decoder, at the receiving end of the channel, associates 
with each received N-vector R one of the members of $ (say X’). Let 
P, = Pr {X’ + &} be the probability of a transmission error, and 
denote the (mean-squared) transmission error by 


Pood... 
G => E||X—X If. (5) 


Clearly 


és P, max||u—v (Ii. (6) 
n u,veS 


Further, the overall error é’ satisfies 
@ = 2B X—X! | < (eo ten). ) 


Thus we want to make both eg and e7 as small as possible. 

Consider the parameter J7, the number of members of the approximat- 
ing set S. In the interest of minimizing eg we want to make J large. 
However, in the interest of minimizing P, and therefore «7, we want 
to make M small. The proper compromise yields our result, equation (3). 

The following theorems indicate just how to choose M. The first is 
proved in Section III, and the second was proved by C. E. Shannon 
(Ref. 3, p. 16). 


Theorem 1: (Source Encoding). Let X be a random n-vector (n = 1, 2, ---) 
whose components are zero mean Gaussian variates with variance o”. 
Let Ry > 0 be given and let {a,}% be a sequence which tends to zero. Then 
there exists a set 8 of M n-vectors, where 


M s exp [n(o — a,)], (8) 
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and a mapping f of X into &, such that asn — © 


1 2 
7B \|X — 7% | 


< ote (1 + 20, + PE") + oat) + of EB"). gy 


Furthermore, & ts such that for all u ¢ §, 
Lal? so. (10) 
i = 


A special case of some interest in itself is that for which a, = 0. 
In this case Theorem 1 asserts the existence of a quantization of X 
with M = exp (nR,) points and mean-squared quantization error no 
more than o” exp (—2R,)(1 + log n/n) + O(log log n/n) asn > &. 
If the channel is noiseless with capacity C nats second, it can transmit 
e°™ messages with no error in the N seconds that it takes for the n- 
vector X to be emitted from the source. Thus if Ry = C/R, M = exp 
(Ron) = e°”, and the members of $ can be reconstructed perfectly by 
the receiver. The overall error is therefore 


pee eee ee 20/R)(1 4; 196 zn) 4 o( a log en), (11) 


Let us turn now to the discrete memoryless channel defined by an 
input set (1, 2, --- , K), an output set (1, 2, --- , J), and a set of transi- 
tion probabilities Pj | kh), 1 S$ 7 S$ J, 1S k S K. Corresponding to 
each input probability distribution p,, 1 < k < K, there is a joint 
distribution p(k, 7) = P(j | k)p, on the product of the input and output 
sets. Define the random variable (called “information’’) 


Ub, ) = log | PUD], iS2b2K 12925. 269 


is De Peis | i’) 


The channel capacity C = max,,,, HU, where the maximization is 
performed with respect to all possible input probability distributions. 
Let (p*)~., be a maximizing input distribution, and let U*(k, 7) be 
the corresponding information, then define 


= (2 var U*)"' = (2EU* — 20771. (13) 


Theorem 2: (Shannon): Let {by}S_-, be a sequence which tends to zero 
from above. Then there exists an N-dimensional code (for the channel 
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described above) with MT members 

M = exp [N(C — by)], (14) 
and (for any a priori distribution on the code words) error probability 

P, S k exp (—6Nby), (15) 
where k is independent of N, and B ts defined by (18). 


Actually we can broaden the class of channels for which our main 
result (3) holds to include that class of channels for which Theorem 2 
holds for some constant 8. This broadened class includes the Gaussian 
channel with signal-to-noise ratio p for which 6 = [(1 + p)’/p(2 + p)] 
(see equation 74 of Ref. 4). 

Theorems 1 and 2 lead us directly to the proper choice of M and 
our main result. Since the channel encoder must encode each of the 
members of $ into channel inputs, we equate the M’s of Theorems 1 
and 2 obtaining (from n = NR) 


Ro, =C/R and by = Ra,. (16) 


_ flog N 
by = | BN ) (17) 


where @ is defined in (13) we have from Theorem 1 a quantization error 


oe Peat 20/R)(1 a aw (eX) + 0 EX) | (18) 


and from (6), (10), and Theorem 2, a transmission error 


If we then choose 


ey < 40°P, < 4o°k ~ (19) 


Thus by combining (7), (18), and (19) we have an over-all mean-squared 
error 


é < oc’ exp (—20/)(1 + ah (ex) aL (x): (20) 


This is our result (8). 


III. PROOF OF THEOREM 1 


We must establish the existence of a mapping f of Euclidean n-space 
into a set 8 of M n-vectors such that E || X — f(X) ||’/n satisfies the 
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upper bound (9). Let 7, = H || X ||. It can be shown that with very 
high probability, as 7 — o, X will be near the surface of the sphere 
of radius r, with center at the origin. It turns out to be convenient 
to establish the existence of a mapping g from the surface S,, of this 
sphere to a set S. 

Accordingly we shall construct the mapping f as follows. Let X, = 
X(E || X{{/]| X |]) = X(r./]| X ||) be the projection of X onto the surface 

r+ Let g be a mapping of S,, to some set § of n-vectors (g has not 
yet been found, of course), then f = g(X,). The following lemma of 
D. J. Sakrisan’, is proved in the Appendix. 


Lemma 1: For any mapping g, and f = g(X,), 
BH \|X — fX) |P = Var || X|]+ 2 || %, — 9%) | 2D) 


Since, as we shall see, Var || X || is relatively small, the principal 
contribution to H || X — f(X) ||’ is Z || X, — g(X,) ||’. 

Our next task is to find the mapping g, and to this end we will establish 
a lemma concerning the covering of S,, by spherical caps. First some 
definitions. 

Let w, z with and without subscripts denote points on S,,, the 
surface of a sphere in n-space of radius r,. Let a(w, z) be the angle* 
between w and z. For 0 S @ S 7, let C(w, 0) = [z : a(w, z) S 94] be 
the spherical cap of half angle @ centered at w. Assign the usual ‘‘area’”’ 
measure to S,,. If A C S,, is measurable, let »(A) be its measure. 
In particular, let 
i Lr} 


( 
C,(9) ae ul C(w, 6)] ae 
n+1 
r( +4) 


be the area (measure) of a cap of half angle @, so that 


F 
/ sin”? » dg (22) 
0 


n/2 n-1 


Cr) = “(tt2) (23) 


Lemma 2:{ Let X, be a random vector which is uniformly distributed 


is the area of S,,. We now state 


* The angle a(w, z) is defined by cos a = (w, z)/|| w || || z ||, where (w, z) is 
the inner product and 0 <a <7. 

7 It is shown in Ref. 5 , that An(r), the surface area of an n-sphere of radius r 
is given by An(r) = nx? 7/T [(n abe 2)/2). Equation (22) follows from the fact 
that Cn(0) = So’ (redq) An- (7. sin ¢). 

+ Lemma 2 is related to a result on the covering of the n-sphere in Ref. 6. 
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on S,,. Let M (a positive integer) and 6 (0 S 6 S zm) be arbitrary. Then 


(for any dimension n) there exists a set of M points {wi ,--- . War} C S,, 
such that 
QW, »7t , War) 


I> 


Pr {Xe U ew, 0} < exp(-MC(O/C@). 24 


Proof: Let us define the function 


i <j<¥-” 
Ce nee re ae ‘ if aw;,z)>6, 1SjK M095) 
0 otherwise. 
Then for a fixed set w,, --: , War, 
M 
Ow, » 7° » War) = Pr i ¢ U C(w; ) as 
J=1 y, 
Tr EF, > °° , War , X,) 
1 ; 
- ze | Pon fee War, 2 du(z). (26) 
Now consider a random experiment in which the JZ points w, ,- ++, Way 


are random vectors chosen independently with uniform distribution 
on S,,. Q is now a random variable given by 


QW, Way 9 Wa) = ais | POW. +, War, 2) ule), 27) 


where upper case W’s represent random vectors. We now compute 
EQ(W,, -:-, War), the average of Q over all choices of W,,--- , War. 
We will show that HQ S exp [—MC,(6)/C,(r)]. Since at least one set 
{w,, °°: , War} must satisfy Q(w,, --:, War) S HQ, the lemma will 
follow. 

We can write 


HQ = Eats | POW, War, 2) due) 
CG) 


= a [ auoBrow, ,---, War, 2), (28) 


the interchange of expectation and integration being justified by the 
fact that F = 0. As indicated in (28), EF(W,, --- , Wa , Z) is computed 
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with z held fixed. Now 


M 
EPW,,++: , War, Z) = Pr {PF = 1} = Pr) (a(W,, z) > 9} 


7=1 


a C,(0)\" 
=(1- tt < exp [—-MC,(0/C.ln)], 


(29) 
so that 


B(Q) < exp (-MC,(0)/Cul@)) | 2 = exp [-MC,(0)/Cule)] (80) 


which concludes the proof. 

We now give the construction of g. Let M = M,, = exp [n(R, — a,)] 
as in Theorem 1, and let 6 = 6, = are sin exp (—R, + 46,), where 
5, > 0 will be specified below. Let {w,, we, --- , War} be a set which 
satisfies (24) for these 6 and M. Let x e S,, and say x e UjL, C(w,, 8). 
Let j, be the smallest index j such that x « C(w; , 6). Then 


g(x) = (cos 8) wy. , (31) 


(see Figure 3), and || g(x) — x |] S 7, sin 6. If x ¢ UjL, C(w;, 6), then 
let g(x) = w,. In this case || x — g(x) || S$ 2r, . Hence 
E || X, — g(%) \|’ < rosin” 6+ 4r3Q(wi , We, +++» War). (32) 


Since the set {w;} satisfies (24), 
E ||X, — 9(X,) ||? S 72 sin’ 6, + 47? exp [—MC,(6)/C,(m)]. (83) 





AN Ut 
ua Seer : 
ee i / 
a tf / 
i i i 


‘ya 


Fig. 3— Construction of g(x). The solid line is the cap C(wy;, 4). 


ANALOG DATA COMMUNICATION 809 


Combining (33), and Lemma 1, and the fact (proved in the Appendix) 
that 


r, <ov/n and * var || X || S 0.9207/n, (34) 
we obtain 
1 _ 2 
1B Ix — 700 || 


< o’ sin’ 6, + 40° exp [—M,C,(6,)/C,(r)] + 0.920°/n, (35) 


where M, = exp [n(R, — a,)] and 6, = arc sin exp (—R, + 6,), where 
5, > 0 is to be specified. Our final step is selection of 6, so that (9) is 
satisfied. For 6, bounded away from 7/2, Shannon (equation 27 of 
Ref. 4) has shown that as n > ©, 


C CalOn) J __ sin" 0, 
C,,(z) eto cos 6, sin 6, Be sin” 6, (36) 


so that (using the definitions of 1/7, and 6,) for n sufficiently large 











C,,(On) exp [n(6, = an) | 
uM, = SS: 37 
Cn) oe en) 
We now define 6, by 
5 <a, 4+ Lieen , loglogn , log V2r, (38) 
2 hh n n 


Then, from (87), for n sufficiently large 
exp [—17,C,(0,)/C,(m)] S 1/n. (39) 
Finally, we have 
exp (26,) = 1 + 26, + O(8:) 
Ai OCs o| (128 ny’ | (40) 
Combining (35), (39), and (40) we obtain 
1 2 
~ B || X — #0) || 


< ote (1 + 20, + 8%) 4 of PBMER) + oe) (an 


which is (9). 
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Inequality (10) follows from (31) and (34) by simply 
[1 #X) |] = |] 9K) |] = 7 cos 8< 1m < Vo". (42) 


APPENDIX 


In this appendix we establish some facts about the random vector 
X = (X,, X2, --- , X,) whose coordinates are independent zero-mean 
Gaussian variates with variance o”. 


Proposition: 
(i) E \| X || S (no*)* 
(ii) var || X || S 0.920” 


Proof: (i) follows from the Schwarz inequality: [Z(1-|| X ||P s 
E\| X |? £V = £ || X |)? = no’. To establish (7), consider the n-fold 


joint probability density for the vector (R, $,, ---, ®,_,), the polar 
coordinate representation of X*, 

g(r, Y1 > f2 Gn—1) = exp (r°/20°)r CH gn a Pn-1)3 (43) 
where h(¢, , *** 5 Gn-1) = (200°)""” cos "~* ¢, Cos” ge ++ * COS Gano. 
The marginal distribution for the random variable R = || X || is 

2 n- 2 2 
10: Se ee Er 20); 
ar) 


so that an integration yields 


BE \|X || = ER) = V2? ret +)/r(3). 


Since # || X ||? = no’, we have 


| 


var || X || Spee 


"L 6) J) 


ww hve (1+ 13 — goat) $00 sow v%R (1 + a), 


* See, for example, M. G. Kendall and A. Stuart, The Advanced Theory of 
Statwstics, vol. 1, London: Griffin, 1963, Section 11.2. 


(44) 


Using the Stirling formula, 
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to underestimate IT'(n + 1/2) and overestimate I'(n/2) we obtain 


Var || X {I < _ -( 1y( = 7 = I) 
ne ONE) Stn 1) SOG ey)? 
Further 


et. +)) 








a 





exp [—1 + n log 1 + 1/n)] 











> exp [—-1 + n(1/n — 1/2n’)] 
1 
= exp (— 1/2n) 2 1 ~~ On? 
and 
] 1 Lote 1 0.21 
Sona 1) M+) =z (ZE+5) SE 
Hence, 
Var || X | 
2 


9 
a Ley GeO S9 0 2 aor) <2 





a (46) 
This is (72). 


We now give a 


Proof of Lemma 1: 
E || X — f(&) | = # ||X-—X,+ X, — g&,) |! 


= E\|X—X, | +2 ||X, — 9) II’ 


+ 2H (X or X, ) X, ri g(X,)), (47) 
where (u, v) is the inner product of the n-vectors 
|X — X, |? = { 


|| X || — # || X |I}", so that 


E || X — X, ||’ = var || X |]. (48) 
Further, the inner product 


u and v. Now 


(x = X, ’ X, = g(X,)) 
= (X, X,) ~~ (X, g(X,)) re (X, ’ X,) She (X, ’ g(X,)) 
= [|X |] ro — ro — || X I |] gCK) || cose, + ro || (Kp) || cos ae , 


(49) 
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where a, is the angle between X and g(X,), and a, is the angle between 
X, and g(X,). Since X and X, are colinear, a, = a, = a(X,) a function 


of X, . Now from (48) we see that the random variable R = || X || 
and the vector (®,, --- ®,_,) are statistically independent. Since X, , 
depends only on ®, --- &,_, and not on R, we conclude that || X || is 


independent of g(X,) and a(X,). Thus from (49) 
E(X — X&, , X, — g(%)) 
= ro || X || — ro — E || X || Ell] g(&,) || cos a(X,)] 
+ roB[|| 9X) |] cos a(X,)] = 0. (50) 
Equations (47), (48) and (50) imply Lemma 1. 
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On Permutation Switching Networks 


By A. E. JOEL, JR. 
(Manuscript received December 1, 1967) 


Rearrangeable switching networks are considered as permutation gen- 
erators. Two-by-two network arrays may be implemented in these networks 
using a B reversing element. Nested tree arrays of these elements may be 
used to synthesize rearrangeable switching networks which appear minimal. 
Multistage network arrays of these elements may be implemented within 
a single coordinate device. 


I. INTRODUCTION 


Two-sided communication switching networks with an equal num- 
ber of input and output terminals (NV) have been looked upon as 
permutation generators of order N.+? Alteration of the permutation 
being produced by any state of the network configuration may be 
- achieved by rearranging or changing the state of the switching net- 
work. By limiting the switches in each stage to binary action devices 
certain economies have been achieved.? This paper describes addi- 
tional savings, particularly in the number of required operating ele- 
ments. Some further savings are achieved by introducing break con- 
tacts into switching networks. 


II. THE 8B ELEMENT 


A small space division switching network is a two-by-two array 
of crosspoints (see Fig. 1). Each of these crosspoints is generally con- 
sidered as either electromechanical make-contacts or electronic two- 
state devices. Confining our attention to make-contacts, they may be 
actuated by a coordinate mechanism, as in a crossbar switch, or as 
individual relays. 

Usually in switching networks the crosspoints are controlled so 
that only one may be operated in any column and row. Such use of 
2 X 2 arrays confine the maximum number of states from the 16 
possible states to those useful in switching networks. The 16 states 
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CROSSPOINT = = 
OPERATIONAL 
ELEMENT 





Fig. 1—2 x 2 crosspoint array. 


include double connections (two operated in the same column or 
row) and the zero state. The zero state is useful if the crosspoint 
devices require holding power. When this state occurs it introduces 
a vacuous partition into the permutation of the network. Actually 
the only useful states, particularly when viewed from the standpoint 
of rearrangeable networks, are the two shown in Fig. 1: (Al, 
B-2), and (A-2, B-1). In these rearrangeable networks the vacu- 
ous states are eliminated. 

It is proposed that the equivalent of 2 x 2 arrays of four make- 
contacts and associated electromagnets be achieved by the simple 
expedient of two transfer contacts connected in the so-called reversing 
configuration (see Fig. 2a). Using this two-state or binary (8) ele- 
ment only two moving contacts and one electromagnet are required. 
The equivalent of 8 elements may be achieved in many other ways 
without break contacts. Figs. 2b and 2c show a @ array made from 
the array of Fig. 1 using only 2 operating elements. 


Ill. THE USE OF B ELEMENTS IN SYMMETRICAL REARRANGEABLE NETWORKS 


Methods have been described by V. E. Bene* for using 2 X 2 ar- 
rays to form rearrangeable link networks. Improvements result when 
applying B elements to these networks. For example, the equivalent 
of a 4 X 4 rearrangeable switching network can be made of three 
stages of 2 X 2 switches or a total of six 8 elements (Fig. 3). Here 
6 B elements replace 16 crosspoints and their operating magnets. An 
8 X 8 switch of 64 crosspoints may be replaced by 20 6 arrays in 5 
stages (see Fig. 4). 

For extending to larger networks, the method taught by Benes 
may be applied, replacing the middle 4 x 4 stage by the three sym- 
metric stages of 8 elements. A rearrangeable network with N terminals 


PERMUTATION SWITCHING NETWORKS 815 








Fig. 2— Binary B element. (a) One operator. (b) Two operators. 


on each side will comprise S stages of these elements with N/2 ele- 
ments per stage where 


S = 2 log, N — 1. 
The total 8 elements per input terminal is 
N 1 
3 @ log, N — 1) = log, N — 3. 
For a network with N = 1024 terminals per side (see Fig. 5), then 
9.5 B elements are used in the network per input terminal. This is 


considerably less than the 82 coordinate array crosspoints per input 
required in the network described by Benes. 


Fig. 3— Three stage 4 & 4 network of 8 elements. 
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| EQUIVALENT ____5 
MIDDLE STAGE 


Fig. 4— Five stage 8 x 8 BeneS type network of 6 elements. 


IV. NESTED TREE TYPE PERMUTATION NETWORKS 


The 4 X 4 symmetric network may be further simplified by omit- 
ting one 8 element to obtain a network (see Fig. 6), where each input 
reaches each output through one or more selecting trees of 6 ele- 
ments. The 8 elements in the first stage is used by inputs 2 and 3 to 
select one of two trees in stages 2 and 3, either of which can reach 
a given output depending upon the state of the B element to which 
it is connected. 

This 4 X 4 network (Fig. 6) is a trivial case of the “nested trees” 
concept. The stages 2 and 3 are two superimposed trees with stage 
1 a one-stage selecting tree. 

The concept is better illustrated with higher order networks. The 
8 X 8 array in Fig. 7 uses 17 B elements with 1, 2, and 3 stage trees 
in a total of 6 stages. This compares very favorably with a single 
stage 8 X 8 array requiring only about one quarter of the number of 
operating magnets, about one half of the number of moving contacts, 
and only 4 more contacts. A 16 xX 16 array with 49 8 elements would 
have 1, 2, 3, and 4 stage trees in a total of 10 stages. : 

In general with nested trees the first stage connected with the out- 
puts (or inputs) has 


Ss, = log, N 


stages of N/2 B elements and is the highest order tree group. Then 
the remaining stages form lower order trees, which might be con- 
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Fig. 5 — Bene’ type 1024 terminal network of 8 elements. 


sidered in different spatial dimensions or disjunctive nested parts of 
the network. 
The total 8 elements in such a network is 


logs N AT N 
> 9° logs ge-1 


This function simplifies to 


N log, N — N +1. 


Derived in another manner this same formula has been obtained by 
Goldstein and Leibholz* 


STAGES-->1 2 3 


Trig. 6 — Three stage “tree type” 4 x 4 network of 8 elements. 
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| STAGE 2 STAGE be re d 
TREE TREES 3 STAGE TREES 


Fig. 7 — Six stage “nested tree” type 8 & 8 network of 8 elements. 


Figure 8 shows the derivation of a nested tree network for N = 32 
where the different order nested trees are shown as boxes to illustrate 
the various parts of the network. Figure 9 illustrates the spatial con- 
cept for the same network. 

The total number of stages in such networks may be defined as 
the total number of stages in all trees. Thus, 





4 STAGES 5 STAGES 
! 
3 STAGES 8's 16B's 
4B's moo000 





32 TOTALS 
STAGES { + 2 + 3 + 4 + 5 = 15 
B 
ELEMENTS 1 + 4 + {2 + 32 + 80 = 129 


Fig. 8 — Schematic of “nested tree” type 32 X 32 network of 8 elements. 


PERMUTATION SWITCHING NETWORKS 819 


S = log, y(lo Xt ), 


The number of 8 elements per input is, for large numbers, loge NV 
— 1. This is an improvement of 0.5 8 element per input compared 
with the symmetric network approach described in the previous sec- 
tion. 

To double the size of a nested network requires the addition of 
new higher order trees and doubling the input capacity of the previous 
network by converting the former output trees to input trees. Increas- 
ing the size by this means is much easier than the reassignment of 
links required when symmetric type networks grow. 


V. SERIAL PERMUTATION NETWORKS 


The nested tree type link network is nonsymmetric. Other non- 
symmetric arrangements of reversing devices have been proposed to 





2 DIRECT 
INPUTS 

k 
_ 32 

| OUTPUTS 

NC 4 STAGES SACS 

‘é 38'S 168 ELEMENTS 

INPUTS 


3 STAGES 











Fig. 9 — Pictorial concept of “nested tree” type 32 X 32 network of 8 elements. 
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form permutation networks.® These networks are scrial in character 
as contrasted to the link type. 

As shown in Fig. 10, assuming a permutation generator of Py — 1, 
then (N — 1) 8 elements will insert the N lead into the permuta- 
tion to produce Py. The total 8 arrays per input for N inputs would 
be 


N-1 
2 


For N = 4 a complete network of 3 stages derived in this manner is 
shown in Fig. 11. For N > 4, serial networks of this type require 
more @ elements than the link type. However, these networks are 
more readily adaptable to networks where N is small, where N is not 
a power of 2, and to obtain network growth in one dimension. How- 
ever, growth is expensive since when JN inputs are present the addi- 
tion of another input requires (NV — 1) £ elements. 

Since the (V — 1) @ elements are operated in only (VN — 1) com- 
binations to achieve the desired output, then savings may be made by 
using combinatorial arrangements. Combinations of log. N — 1 ele- 
ments are achieved by trees of N — 1 stages where the several £ ele- 
ments of a stage may be actuated by the same magnet. Thus, for the 
third stage of Fig. 11, we have Fig. 12, where the two of the £ ele- 
ments are actuated by the same magnet. 

For N = 8, although 33 B elements would be required (almost 
twice the number required for a link type network) one less magnet 
would be required, 17 instead of 18. As N increases the magnet sav- 
ings become large, but generally it is impractical to actuate the 
large number of B elements required for the higher order stages.® 
For N = 1024, for example, the last stage would require 10 relays, 
3 of which would actuate 128, 256, and 512 6 elements, respectively. 

The smaller the number of paths through an array or switch of a 








Fig. 10— Serial permutation network extension. 
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Fig. 11— Three stage serial 4 X 4 network of 8 array. 


switching stage, the fewer existing elements of a permutation are 
disturbed by rearranging. The smaller the switches, the greater the 
number of stages required. However, the number of stages does not 
increase as rapidly as the switch size decreases. The @ array is ideal 
for rearrangeable switching networks, and the £ array is a contribu- 
tion to keeping the number of operating elements (magnets) small. 


S3 
CT 


| 
| | 
| | 
[ead 


Fig. 12 — Stage 3 of serial 4 x 4 network combining 8 elements. 
VI. A NETWORK DEVICE 


The £ elements in the various arrays just described may be thought 
of as crosspoints in a coordinate device. For example, the 8 X 8 
network has 4 rows and 5 columns of @ elements. A 4 X 5 coordinate 
device could be built which has 4 select and 5 hold magnets or operat- 
ing elements, like a crossbar switch. First the select elements would 
be operated in combination to prepare to change the state of £ ele- 
ments in a particular column. Then the column operating element 
would reset the 8 elements in this combination. In five such steps 
the entire 8 X 8 network would be reconfigured. For such a large array 
as a 1024 x 1024 symmetrical network 19 steps would be required. 

Crossbar switches with vertical multiples divided such as in Fig. 
2c, could be used as network devices in this sense. The interstage 
link wiring would be between verticals. A 128 xX 128 network requires 
769 8 elements. A 10 horizontal crossbar switch provides 5 B elements 
per vertical so that 16 10 xX 10 switches would implement this size 
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rearrangeable network. This compares with 30 10-vertical crossbar 
switches for the more familiar three-stage network of 10 x 10 switches 
fora 100 xX 100 rearrangeable network. 

Nonblocking networks without interruptions caused by rearrange- 
ability can be made from rearrangable networks by extension of 
switch sizes to a full Clos network. For a 100 x 100 network, 20 
10 X 19 and 19 10 X 10 switches are required. Nonblocking without 
interruption may also be achieved by alternately using two permuta- 
tion networks connected in parallel. For the 128 x 128 network only 
36 10 X 10 switches with divided vertical multiples would be required 
in this configuration. 


VII. SUMMARY 


The B element reduces the number of devices required to imple- 
ment rearrangeable switching networks. In particular, networks con- 
sisting of nested trees of these elements produce minimal permutation 
switching networks. Crossbar switches or similar coordinate devices 
may be used to implement networks made of @ arrays. Two permuta- 
tion networks in parallel are equivalent of a nonblocking switching 
network without rearrangement. 


REFERENCES 


1. BeneS, V. E., “Permutation Groups, Complexes, and Rearrangeable Connect- 
ing Networks,” B.S.T.J., 43, No. 4, part 2 (July 1965), pp. 1619-1640. 

2. Joel, Jr. A. E, “Relay Permutation Type Switching System,” U S. Patent 
2,625,610, applied for December 5, 1947, issued January 13, 1953. 

3. BeneS, V. E., “Optimal Rearrangeable Multistage Connecting Networks,” 
BS.T.J., 43, No. 4, part 2 (July 1964), pp. 1641-1656. 

4, Goldstein, L. J., and Leibholz, 8. W., “On the Synthesis of Signal Switching 
Networks with Transient Blocking,” IEEE Trans. Elect. Comm. (October 
1967), pp. 637-641. 

5. Moore, E. F., “Relay Selecting Circuit,” U S. Patent 2,864,008, applied for 
December 28, 1955, issued December 9, 1958. 

6. Shannon, Claude E., “Memory Requirements in a Telephone Exchange,” 
BS.TJ., 29, No.3 (July 1950), pp. 343-349. 


Contributors to This Issue 


E. T. Bayuiss, B.S.E.E., 1961, and M.S., 1962, both from Oregon 
State University; Bell Telephone Laboratories, 1962—. Mr. Bayliss 
has worked on methods of antenna pattern synthesis with special ap- 
plication to large array antennas. He is engaged in system engineer- 
ing and design of radars and interceptors for ballistic missile defense. 
Member, Eta Kappa Nu, Pi Mu Epsilon, Sigma Tau, IEEE. 


H. RicHarp Brurrier, Diploma in Radio Engineering, Valparaiso 
Technical Institute, Valparaiso Ind., 1954; Bell Telephone Labora- 
tories, 1954—. Mr. Beurrier has worked on circuits applied to ultrasonic 
transmission systems. More recently he has been concerned with vhf 
and uhf network systems applied to transistor circuits to effect un- 
usually large power capability. 


Ta-Suinc Cuu, B.S., 1955, National Taiwan University; M.S., 
1957, and Ph.D., 1960, Ohio State University ; Bell Telephone Labora- 
tories, 1963—. Mr. Chu has been concerned with problems involving 
electromagnetic wave propagation and microwave antennas. Membcr, 
IEEE, Commission 2 of URSI, American Physical Society, Sigma Xi, 
Pi Mu Epsilon. 


JaMEs C. Daty, B.S., 1960, University of Connecticut; M.E.E., 
1962, and Ph.D., 1967, Rensselaer Polytechnic Institute; faculty 
member of the Electrical Engineering Department at Rensselaer, 
where he taught and did research on bulk semiconductor microwave 
interactions, 1962-1966. Bell Telephone Laboratories, 1966—. He is 
now concerned with optical wave guidance. Member, IEEE, Tau 
Beta Pi, Eta Kappa Nu, Sigma Xi. 


ALLEN N. FriepMAN, B.S., 1956, City College of New York; MS., 
1958, Yale University; Ph.D., 1964, Columbia University; Bell Tele- 
phone Laboratories, 1961—. Mr. Friedman first studied the properties 
of metalic magnetic materials, then became involved in exploratory 
development of devices between optical and subaudio frequencies. Re- 
cently he has been concentrating on vhf and uhf devices. Member, 
American Physical Society and its Solid State section, Phi Beta 
Kappa, Sigma Xi. 


823 


824 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 


Deter Guoce, Dipl. Ing., 1961, D.E.E., 1964, Braunschweig Tech- 
nische Hochschule (Germany); research staff, Braunschweig Tech- 
nische Hochschule 1961-1965; Bell Telephone Laboratories, 1965—. 
In Braunschweig, Mr. Gloge was engaged in research on lasers and 
optical components. At Bell Telephone Laboratories, he has concen- 
trated in the study of optical transmission techniques. Member, VDE, 
IEEE. 


D. C. Hoas, B. Sc., 1949, University of Western Ontario; M. Sc., 
1950, Ph.D., 1958, McGill University; Bell Telephone Laboratories, 
1953—. His work has included studies of artificial dielectrics for 
microwaves, diffraction of microwaves, and over-the-horizon, milli- 
meter wave and optical propagation. Fellow, IEEE; member, Com- 
mission 2 of URSI, Sigma Xi, AAAS. 


Amos E. Jort, Jr., B.S., 1940, and M.S., 1942, Massachusetts Insti- 
tute of Technology; Bell Telephone Laboratories, 1940—. Mr. Joel did 
relay engineering, crossbar system testing, fundamental development 
studies of telephone switching systems, design of circuits for early relay 
digital computers, and preparation of texts to teach switching design. 

Later he was concerned with designing automatic message accounting 
computer circuits, making fundamental engineering studies of new 
switching systems, and was head of a department responsible for the 
development planning of the Bell System’s first electronic telephone 
switching systems. Mr. Joel was Director of Bell Laboratories’ Com- 
mon Systems Switching Laboratory from 1961 to 1967. Mr. Joel is 
now a Switching Consultant at Bell Laboratories. He holds some 50 
patents, among which is the largest U. S. patent ever issued (its 884 
pages cover a machine to automatically compute the cost of each toll 
call). 

He has been active in national and local IEEE affairs, is a member- 
at-large of its Communication Technology Group, chairman of its 
Meetings Committee, and member of its Communication Switching 
Committee. Member, Sigma Xi, Association for Computing Machinery, 
AAAS. 


Prerer Kaiser, Diplom Ingenieur, 19638, Technische Hochschule, 
Munich, Germany; M.S., 1965, and Ph.D., 1966, University of Cali- 
fornia, Berkeley; Bell Telephone Laboratories, 1966—. At Berkeley, 
Mr. Kaiser was working on frequency independent antennas. He now 
is engaged in optical transmission research with emphasis on gas lens 
beam waveguides. Member, IEEE. 


CONTRIBUTORS TO THIS ISSUE 825 


R. E. Parkin, B.Sc. (Eng.), 1962, Ph.D., 1965, Imperial College 
(University of London); Bell Telephone Laboratories, 1965—. Mr. 
Parkin is concerned with computer-aided design techniques applied 
to telephone development, with particular emphasis on hybrid inte- 
grated technology. Associate Member, IEE. 


H. Srwen, B.E.E., 1948, City College of New York; M.E.E., 1947, 
and D.E.E., 1954, Brooklyn Polytechnic Institute; Bell Telephone Lab- 
oratories, 1953-1961, 1962—. He has been concerned with microwave 
ferrite propagation studies and parametric amplification in solid-state 
devices, and active networks in high through microwave frequencies. 


WiuuiaM H. Srerer, B.S.E.E., 1955, Evansville College; M.S.E.E., 
1957, and Ph.D. (E.E.), 1960, University of Illinois; Bell Telephone 
Laboratories, 1962—. Mr. Steier first worked on the millimeter wave 
circular waveguide transmission system. More recently he worked 
on optical transmission lines and gas lenses. Member, IEEE, Sigma Xi. 


Aaron D. Wyner, B. 8., 1960, Queens College; B.S.E.E., 1960, M.S., 
1961, and Ph.D., 1963, all from Columbia University; Bell Telephone 
Laboratories, 1963—. Mr. Wyner has been engaged in research in 
various aspects of information theory. He is also adjunct associate 
professor of electrical engineering at Columbia University and Chair- 
man of the Metropolitan New York Chapter of the IEEE Information 
Theory Group. Member, IEEE, SIAM, Tau Beta Pi, Eta Kappa Nu, 
Sigma Xi. 


